Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,183,802 A→G D853G (GAC→GGC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,8020AG98.5% 250.2 / ‑3.8 68D853G (GAC→GGC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (35/32);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGC  >  NC_000913/4183656‑4183948
                                                                                                                                                  |                                                                                                                                                  
cATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATc                                                                                                                                                <  1:250049/151‑1 (MQ=255)
     cGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGa                                                                                                                                            >  2:80854/1‑150 (MQ=255)
      gTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATccc                                                                                                                                              <  1:568571/147‑1 (MQ=255)
         gAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTg                                                                                                                                       >  1:155946/1‑151 (MQ=255)
          aaTGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTgg                                                                                                                                      >  1:156783/1‑151 (MQ=255)
               ttACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGa                                                                                                                                  >  2:341916/1‑150 (MQ=255)
                tccAACTTCGATGTCTCCATCCTCGTATCCGAGCGCATTGTTCAGTAAAACCGTTTCACCACCTTCCACACTCAGGCACTGGCGCGTGTGTCAAGTGCCACCTAGCTGAATCCTGAAGAGATCAGCGCGGGCAGCCCGAACGTGGGTGaa                                                                                                                                 <  2:170497/148‑1 (MQ=255)
                tACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTg                                                                                                                                   >  2:512615/1‑148 (MQ=255)
                 aCAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGc                                                                                                                               >  1:259403/1‑151 (MQ=255)
                  cAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCt                                                                                                                              <  1:472549/151‑1 (MQ=255)
                    aCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTgc                                                                                                                            <  1:57249/151‑1 (MQ=255)
                    aCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTgc                                                                                                                            <  1:164826/151‑1 (MQ=255)
                    aCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATGCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGc                                                                                                                               <  2:870687/148‑1 (MQ=255)
                      ttCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCCGGCATCCCGAACGTGGGTGAAGCTgcg                                                                                                                           <  2:18846/150‑1 (MQ=255)
                       tCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTgcgc                                                                                                                          >  2:971548/1‑150 (MQ=255)
                        cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGGACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCTCTGGCATCCCGAACGTGGGTGAAGCTgcg                                                                                                                           >  2:61585/1‑148 (MQ=255)
                        cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGctc                                                                                                                        >  2:564523/1‑151 (MQ=255)
                        cGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGctc                                                                                                                        >  2:635021/1‑151 (MQ=255)
                                 cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                               <  1:41156/151‑1 (MQ=255)
                                 cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                               <  1:549789/151‑1 (MQ=255)
                                 cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTg                                                                                                               >  2:847210/1‑151 (MQ=255)
                                  aTCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGAGCTCTCCAAACTgg                                                                                                              >  2:823062/1‑151 (MQ=255)
                                    ccTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAt                                                                                                            >  1:24262/1‑151 (MQ=255)
                                            ccGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCgg                                                                                                    >  1:791226/1‑151 (MQ=255)
                                              gAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                                                                                   >  1:224519/1‑150 (MQ=255)
                                              gAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTa                                                                                                  <  2:618290/151‑1 (MQ=255)
                                               aGCGTTTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt                                                                                                 >  2:840884/1‑151 (MQ=255)
                                                  gtgtTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGt                                                                                              >  1:290139/1‑151 (MQ=255)
                                                    gttgttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGttt                                                                                            >  2:333265/1‑151 (MQ=255)
                                                       gttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGCCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTAc                                                                                          >  1:237241/1‑150 (MQ=255)
                                                             gAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTg                                                                                   >  2:250506/1‑151 (MQ=255)
                                                                   cGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAg                                                                             <  1:336410/151‑1 (MQ=255)
                                                                    gTTTCACCACCATCCACATTCAGGAGCTGGCGTGGGTGTCCCGCGACACCAAGCTGGGTCCGGAAGATATCACCGCTGGCATCCCGAACGTGGGTGATGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGt                                                                            <  2:47779/151‑1 (MQ=255)
                                                                             ccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTgg                                                                   >  1:851823/1‑151 (MQ=255)
                                                                             ccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTg                                                                    <  2:509195/150‑1 (MQ=255)
                                                                               aTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg                                                                 <  2:197203/151‑1 (MQ=255)
                                                                                  cacaTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACa                                                              <  1:221533/151‑1 (MQ=255)
                                                                                   acaTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGTTCACCGCTGGCATCCAGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCgtcgt                                                             >  2:724103/1‑147 (MQ=255)
                                                                                   acaTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt                                                             >  1:509004/1‑151 (MQ=255)
                                                                                    caTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt                                                            >  1:245795/1‑151 (MQ=255)
                                                                                      ttCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCt                                                          >  1:336944/1‑151 (MQ=255)
                                                                                      ctCAGGAACTGGCGTGTGTGTCCCGGGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTGTCGTTTAGATTTGTTCGGAAGTGACCGGTTGCGACa                                                              >  2:467052/2‑147 (MQ=255)
                                                                                       tCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCt                                                          <  1:97676/150‑1 (MQ=255)
                                                                                        cAGGAACTGGCGTGTGTGTCCCGTGACCCCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCt                                                          <  2:440857/149‑1 (MQ=255)
                                                                                          ggAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggtt                                                      <  2:155924/151‑1 (MQ=255)
                                                                                                ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaag                                                <  1:446500/151‑1 (MQ=255)
                                                                                                ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaag                                                <  2:71956/151‑1 (MQ=255)
                                                                                                ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                  <  2:449761/149‑1 (MQ=255)
                                                                                                   gtgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt                                              >  1:75025/1‑150 (MQ=255)
                                                                                                    tgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggta                                             >  2:291219/1‑150 (MQ=255)
                                                                                                     gtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggtaa                                            <  2:187647/150‑1 (MQ=255)
                                                                                                     gtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAAc                                           <  1:318219/151‑1 (MQ=255)
                                                                                                       gtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGc                                         >  2:645449/1‑151 (MQ=255)
                                                                                                         gtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGcc                                        <  1:566003/150‑1 (MQ=255)
                                                                                                           cccGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaa                                     <  1:791985/151‑1 (MQ=255)
                                                                                                           cccGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaa                                     <  1:461648/151‑1 (MQ=255)
                                                                                                              gTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAgg                                  <  1:282077/151‑1 (MQ=255)
                                                                                                                gACACCAAGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGt                                 <  2:824691/150‑1 (MQ=255)
                                                                                                                      aaGCTGGGTCCGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTc                          >  1:76830/1‑151 (MQ=255)
                                                                                                                               ccGGAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAccc                 <  1:375398/151‑1 (MQ=255)
                                                                                                                                 ggAAGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc                >  2:817562/1‑150 (MQ=255)
                                                                                                                                   aacagatcccatcTGGCATCCCGAACGTGGGGCCAGCTGCGACCTCCAAACTGGATGAATCCGGTATCCTTTTCATTGGTGCGGACGTGCCCGGTGGCGCCATTCTGGTTGGCAAGGTAACGCCGCAAGGTCAAACTCAGCTGAcccc                >  2:509317/13‑148 (MQ=255)
                                                                                                                                   aaGAGATCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc                >  1:587727/1‑148 (MQ=255)
                                                                                                                                        aTCACCGCTGGCATCCCGAACGTGGGTGCAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAgaagaa         >  2:196225/1‑150 (MQ=255)
                                                                                                                                        aTCACCGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGGTGACCCCAGAAGaaa        <  2:502941/151‑1 (MQ=255)
                                                                                                                                            ccGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAct    >  2:748974/1‑151 (MQ=255)
                                                                                                                                             cGCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActg   <  1:143670/151‑1 (MQ=255)
                                                                                                                                              gCTGGCATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActgc  >  2:628872/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CATGCCGTGGAATGGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGC  >  NC_000913/4183656‑4183948

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑