Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,696 G→C 6.9% Q22E (CAG→GAG)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,6960GC 6.9% 396.5 / 15.5 128Q22E (CAG→GAG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (50/69);  new base C (5/4);  total (55/73)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-02

Only 100 of 132 total aligned reads displayed.

CATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  NC_000913/4082549‑4082845
                                                                                                                                                   |                                                                                                                                                      
ctttGGGGTTGTCGTACACCCAAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                                                                      <  2:793757/148‑1 (MQ=255)
cATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGGCTGAg                                                                                                                                                      >  1:150015/1‑150 (MQ=255)
  ttGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                                                                                                    >  2:822534/1‑150 (MQ=255)
  ttGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                                                                                                    <  2:146614/150‑1 (MQ=255)
  ttGGGGTTGTAGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGCCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgc                                                                                                                                                    <  1:228903/150‑1 (MQ=255)
   tGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGgcg                                                                                                                                                   >  1:1133031/1‑150 (MQ=255)
   tGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCAGAGgcg                                                                                                                                                   >  1:554220/1‑150 (MQ=255)
     gggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGcgcg                                                                                                                                                 <  2:21340/150‑1 (MQ=255)
      ggTTGTCGTACACCCCATTGTTTTTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                                                                                                                <  2:631068/150‑1 (MQ=255)
      ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGACCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgg                                                                                                                                                <  1:828991/150‑1 (MQ=255)
      ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTGAGGAGCgg                                                                                                                                                >  1:170089/1‑150 (MQ=255)
       gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggg                                                                                                                                               <  2:758737/150‑1 (MQ=255)
       gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                               <  1:1396146/150‑1 (MQ=255)
         tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                             >  1:263064/1‑150 (MQ=255)
         tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                             >  2:67426/1‑150 (MQ=255)
          gTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGt                                                                                                                                            >  2:418473/1‑150 (MQ=255)
          gTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                               >  1:615201/1‑147 (MQ=255)
          gTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                               <  2:615201/147‑1 (MQ=255)
                  ccccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                                    <  1:55298/150‑1 (MQ=255)
                  ccccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCGGAGGCGCGGGGGTCAGAcccg                                                                                                                                    >  2:369380/1‑148 (MQ=255)
                    ccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTa                                                                                                                                  >  1:204493/1‑150 (MQ=255)
                     cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGtttt                                                                                                                                 >  2:359139/1‑148 (MQ=255)
                     aattGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                                                                                                                 <  2:699348/147‑1 (MQ=255)
                       aTGTTATTGCCGACGGTATCGCGAATGTCGGTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTGCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑c                                                                                                                              <  1:97817/150‑1 (MQ=255)
                         gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTATTACTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGA‑cgc                                                                                                                            >  2:1033135/1‑150 (MQ=255)
                           tATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCgg                                                                                                                          >  2:578304/1‑150 (MQ=255)
                           tATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCG‑CACTTCTTCCTGGAAGTCCCGCGCCGGCGGCGCGGGGGTCAGAACGTTAGT‑CGCggc                                                                                                                         >  1:654054/1‑149 (MQ=37)
                               gCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAacg                                                                                                                      >  2:919338/1‑150 (MQ=255)
                                ccGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAacgc                                                                                                                     >  2:465885/1‑149 (MQ=255)
                                  gACGGCATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCGGGGGCGGGGGGGGCAGGCCGGTTGT‑CGCGGCacgccg                                                                                                                   >  1:1100230/1‑150 (MQ=21)
                                      gTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGT‑CGCGGCACGACGAATg                                                                                                               >  2:496200/1‑150 (MQ=255)
                                       tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGa                                                                                                              <  1:962638/150‑1 (MQ=255)
                                         tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑TGCGGAACGACGAATGata                                                                                                            <  1:629132/150‑1 (MQ=255)
                                         tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGata                                                                                                            <  2:744549/150‑1 (MQ=255)
                                          cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGatat                                                                                                           >  2:783917/1‑150 (MQ=255)
                                           gcgAATGTCGTTCCACTCTGAACACGCCACCTGACTGGCTTTACAGCCGATTCAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATTATATc                                                                                                          <  2:1353599/150‑1 (MQ=255)
                                              aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTTACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGTGGCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTg                                                                                                       >  1:1089907/1‑150 (MQ=255)
                                                 gTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGGGCGGGGGTCAGACCGGTAGT‑CGCGGAACGACGAATGATATCTTgaga                                                                                                    >  1:1362057/1‑150 (MQ=255)
                                                  tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGT‑CGCGGAACGACGAATGATATATTgagaa                                                                                                   >  1:1041630/1‑149 (MQ=255)
                                                    gTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTg                                                                                                 <  1:919338/150‑1 (MQ=255)
                                                    gTTCCACTCTCAACACCCCACCTGACAGGCTTTACACCCGATACAGGTGGTAACGTCGATGAGTATCGACACTTCTTCCTGTAAGTCCCGCGCCTGAGGCGCGGGGGTCGGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTg                                                                                                 <  1:1045810/150‑1 (MQ=255)
                                                     ttCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGa                                                                                                <  2:1133031/150‑1 (MQ=255)
                                                       ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              <  2:230739/150‑1 (MQ=255)
                                                         aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGCCCGGTAGT‑CGCGGAACCACGAATGATATCTTGAGCTTGATAAg                                                                                            >  1:402379/1‑150 (MQ=37)
                                                          ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGc                                                                                           <  1:57723/150‑1 (MQ=255)
                                                              gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata                                                                                       >  2:1030141/1‑150 (MQ=255)
                                                                acacGCCACCTGACAGGCTTTACATCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg                                                                                     <  2:1636756/150‑1 (MQ=255)
                                                                    gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt                                                                                 >  2:70599/1‑150 (MQ=255)
                                                                    gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt                                                                                 >  2:400639/1‑150 (MQ=255)
                                                                     ccACCTGTCAGGCTTTACAGCCGATACAGGTGGTTACGTCGATGAGTTTCGCCAGTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGTTTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTTAGATTGATAATCCATATgtcgtc                                                                                <  2:304012/150‑1 (MQ=38)
                                                                      ccccTGACACGCTTTCCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtct                                                                               <  1:316145/148‑1 (MQ=255)
                                                                        ccTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                                                                             >  1:331809/1‑150 (MQ=255)
                                                                         cTGACAGGCGTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCg                                                                            <  2:452923/150‑1 (MQ=255)
                                                                             cAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTac                                                                        <  1:999490/150‑1 (MQ=255)
                                                                               ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGGGGGGGTCAGACCGTTAGT‑CGCGTAACGACGAATTATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTacac                                                                      <  2:1708990/150‑1 (MQ=255)
                                                                                gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACAcc                                                                     <  2:509276/150‑1 (MQ=255)
                                                                                 gtttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  2:1249489/149‑1 (MQ=255)
                                                                                 cTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  1:418473/150‑1 (MQ=255)
                                                                                 cTTTACAGCCGATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGTGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  2:81887/150‑1 (MQ=255)
                                                                                    tccAGCCGATACAGGTGGTAACGTCGATGAGGTTCGCCACATCTACCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt                                                                 <  2:628523/148‑1 (MQ=255)
                                                                                      cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCGGAGGCGCGGGGGGCAGACCGTTAGT‑AGCGGGACCACGAAATATATATTGAGATAGATAAGCCAAATGTCGTATACGTTAAACCatttaa                                                               >  2:602078/1‑148 (MQ=255)
                                                                                         ccGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcaca                                                            <  2:49549/150‑1 (MQ=255)
                                                                                          gatctcAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAt                                                           <  1:438386/145‑1 (MQ=255)
                                                                                          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAt                                                           >  2:1234733/1‑150 (MQ=255)
                                                                                          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTCCCACAt                                                           <  2:331809/150‑1 (MQ=255)
                                                                                               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAAAGATATATTGAGATTGATAAGCCATATGTAGTCTACGTTACACCTTTTAcacaaagata                                                      >  2:438386/1‑144 (MQ=255)
                                                                                                 ggtggtAACGTCGATGAGTTTCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag                                                    <  2:554220/150‑1 (MQ=255)
                                                                                                  gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaagg                                                   <  1:822534/150‑1 (MQ=255)
                                                                                                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaagg                                               <  1:304098/150‑1 (MQ=255)
                                                                                                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaagg                                               <  1:359139/150‑1 (MQ=255)
                                                                                                           tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTTCCCCTTTTCCACATTGACAAGGAAGGACTTa                                          >  1:716651/1‑150 (MQ=255)
                                                                                                           tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTa                                          <  2:1700846/150‑1 (MQ=255)
                                                                                                            ggttGAGTTTCGGCACTTCTTCGTGGAAGTCCGGCGCCTGATGCGCGGGGGTCAGACCGTTAGT‑CGCGGATCGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                         <  2:1509823/147‑1 (MQ=255)
                                                                                                            cGATGAGTTTCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                         <  2:277832/150‑1 (MQ=255)
                                                                                                            cGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAAGGCTATCTTGAGAGTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                         >  1:1033417/1‑150 (MQ=255)
                                                                                                              aTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAAc                                       >  1:896518/1‑150 (MQ=255)
                                                                                                                gAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTc                                     <  2:231736/150‑1 (MQ=255)
                                                                                                                   tttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGTGGCGCGGGGGTCAGCCCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGGCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCgg                                  <  1:1251819/150‑1 (MQ=255)
                                                                                                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGc                                 <  1:793631/150‑1 (MQ=255)
                                                                                                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGc                                 <  2:1100230/150‑1 (MQ=255)
                                                                                                                     tCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCg                                <  2:1693136/150‑1 (MQ=255)
                                                                                                                      cGCCACTTCTTCCTGGAAGTCCCCCCCCTGAGGCGCGGGGGTCAGACCGTTAGT‑AGCGGAACGACCAATGATATCTTTAGATTGATAAGCCCAAAGTCGTCACCGTTACACCTTTTACACAATAACAAAGAAGGACTTTAACTACCGCGt                               >  2:1387804/1‑150 (MQ=255)
                                                                                                                       gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTc                              <  2:1719684/150‑1 (MQ=255)
                                                                                                                        ccACTTCTTCCTGGAAGTCCCGCGCCGGGGGCGCGGGGGTGAGAACGGTTGTAGGGGGAACG‑GGAATTATATTGTTTTATTTATTAAGCATATGTTGTATCCCTTTAACCGTTTGCACATTAACAAGGAAGGGCGTTAACTcccgggggt                             >  2:1303417/1‑143 (MQ=255)
                                                                                                                           cttcttCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                          <  1:1420601/150‑1 (MQ=255)
                                                                                                                           cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                          <  1:67426/150‑1 (MQ=255)
                                                                                                                            ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                         <  1:1535974/150‑1 (MQ=255)
                                                                                                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgtg                        <  1:1129615/150‑1 (MQ=255)
                                                                                                                               ttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCCGACCGTTATT‑CGCGGAACGCCGAATGATATCTTGAGATTGATAAGCCATATTTCGTCTCCGTTACACCGTTTCCACATTAACAAGGAAGGACTTAAAAAACGGCGTAAGCGTGtt                      >  2:841243/1‑150 (MQ=255)
                                                                                                                                 ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg                    <  2:1421560/150‑1 (MQ=255)
                                                                                                                                      gacccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCa               >  2:1131394/5‑150 (MQ=255)
                                                                                                                                        cccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             >  2:803221/3‑150 (MQ=255)
                                                                                                                                        cccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             >  1:1313609/3‑150 (MQ=255)
                                                                                                                                          cccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCaa           >  1:401531/1‑150 (MQ=255)
                                                                                                                                              cgcCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGaa       <  1:1387804/150‑1 (MQ=255)
                                                                                                                                                ccggAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     <  1:1342051/147‑1 (MQ=255)
                                                                                                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     >  1:332437/1‑150 (MQ=255)
                                                                                                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     >  1:853056/1‑150 (MQ=255)
                                                                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCg  <  1:70599/150‑1 (MQ=255)
                                                                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCg  <  2:1712877/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGT‑CGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  NC_000913/4082549‑4082845

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