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.

CATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  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)
                       aTGTTATTGCCGACGGTATCGCGAATGTCGGTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTGCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTc                                                                                                                              <  1:97817/150‑1 (MQ=255)
                         gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTATTACTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGAcgc                                                                                                                            >  2:1033135/1‑150 (MQ=255)
                           tATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCgg                                                                                                                          >  2:578304/1‑150 (MQ=255)
                           tATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCACTTCTTCCTGGAAGTCCCGCGCCGGCGGCGCGGGGGTCAGAACGTTAGTCGCggc                                                                                                                         >  1:654054/1‑149 (MQ=37)
                               gCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacg                                                                                                                      >  2:919338/1‑150 (MQ=255)
                                ccGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgc                                                                                                                     >  2:465885/1‑149 (MQ=255)
                                  gACGGCATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCGGGGGCGGGGGGGGCAGGCCGGTTGTCGCGGCacgccg                                                                                                                   >  1:1100230/1‑150 (MQ=21)
                                      gTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGCACGACGAATg                                                                                                               >  2:496200/1‑150 (MQ=255)
                                       tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                                                                                                              <  1:962638/150‑1 (MQ=255)
                                         tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTTGCGGAACGACGAATGata                                                                                                            <  1:629132/150‑1 (MQ=255)
                                         tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata                                                                                                            <  2:744549/150‑1 (MQ=255)
                                          cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                                                                                                           >  2:783917/1‑150 (MQ=255)
                                           gcgAATGTCGTTCCACTCTGAACACGCCACCTGACTGGCTTTACAGCCGATTCAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATTATATc                                                                                                          <  2:1353599/150‑1 (MQ=255)
                                              aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTTACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGTGGCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTg                                                                                                       >  1:1089907/1‑150 (MQ=255)
                                                 gTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGGGCGGGGGTCAGACCGGTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                                    >  1:1362057/1‑150 (MQ=255)
                                                  tCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGAACGACGAATGATATATTgagaa                                                                                                   >  1:1041630/1‑149 (MQ=255)
                                                    gTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg                                                                                                 <  1:919338/150‑1 (MQ=255)
                                                    gTTCCACTCTCAACACCCCACCTGACAGGCTTTACACCCGATACAGGTGGTAACGTCGATGAGTATCGACACTTCTTCCTGTAAGTCCCGCGCCTGAGGCGCGGGGGTCGGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg                                                                                                 <  1:1045810/150‑1 (MQ=255)
                                                     ttCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa                                                                                                <  2:1133031/150‑1 (MQ=255)
                                                       ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              <  2:230739/150‑1 (MQ=255)
                                                         aCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGCCCGGTAGTCGCGGAACCACGAATGATATCTTGAGCTTGATAAg                                                                                            >  1:402379/1‑150 (MQ=37)
                                                          ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc                                                                                           <  1:57723/150‑1 (MQ=255)
                                                              gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata                                                                                       >  2:1030141/1‑150 (MQ=255)
                                                                acacGCCACCTGACAGGCTTTACATCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg                                                                                     <  2:1636756/150‑1 (MQ=255)
                                                                    gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt                                                                                 >  2:70599/1‑150 (MQ=255)
                                                                    gCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcgt                                                                                 >  2:400639/1‑150 (MQ=255)
                                                                     ccACCTGTCAGGCTTTACAGCCGATACAGGTGGTTACGTCGATGAGTTTCGCCAGTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGTTTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTTAGATTGATAATCCATATgtcgtc                                                                                <  2:304012/150‑1 (MQ=38)
                                                                      ccccTGACACGCTTTCCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtct                                                                               <  1:316145/148‑1 (MQ=255)
                                                                        ccTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                                                                             >  1:331809/1‑150 (MQ=255)
                                                                         cTGACAGGCGTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCg                                                                            <  2:452923/150‑1 (MQ=255)
                                                                             cAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTac                                                                        <  1:999490/150‑1 (MQ=255)
                                                                               ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGGGGGGGTCAGACCGTTAGTCGCGTAACGACGAATTATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTacac                                                                      <  2:1708990/150‑1 (MQ=255)
                                                                                gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACAcc                                                                     <  2:509276/150‑1 (MQ=255)
                                                                                 gtttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  2:1249489/149‑1 (MQ=255)
                                                                                 cTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  1:418473/150‑1 (MQ=255)
                                                                                 cTTTACAGCCGATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGTGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  2:81887/150‑1 (MQ=255)
                                                                                    tccAGCCGATACAGGTGGTAACGTCGATGAGGTTCGCCACATCTACCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt                                                                 <  2:628523/148‑1 (MQ=255)
                                                                                      cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCGGAGGCGCGGGGGGCAGACCGTTAGTAGCGGGACCACGAAATATATATTGAGATAGATAAGCCAAATGTCGTATACGTTAAACCatttaa                                                               >  2:602078/1‑148 (MQ=255)
                                                                                         ccGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcaca                                                            <  2:49549/150‑1 (MQ=255)
                                                                                          gatctcAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAt                                                           <  1:438386/145‑1 (MQ=255)
                                                                                          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAt                                                           >  2:1234733/1‑150 (MQ=255)
                                                                                          cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTCCCACAt                                                           <  2:331809/150‑1 (MQ=255)
                                                                                               cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAAGATATATTGAGATTGATAAGCCATATGTAGTCTACGTTACACCTTTTAcacaaagata                                                      >  2:438386/1‑144 (MQ=255)
                                                                                                 ggtggtAACGTCGATGAGTTTCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag                                                    <  2:554220/150‑1 (MQ=255)
                                                                                                  gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaagg                                                   <  1:822534/150‑1 (MQ=255)
                                                                                                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaagg                                               <  1:304098/150‑1 (MQ=255)
                                                                                                      tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaagg                                               <  1:359139/150‑1 (MQ=255)
                                                                                                           tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTTCCCCTTTTCCACATTGACAAGGAAGGACTTa                                          >  1:716651/1‑150 (MQ=255)
                                                                                                           tCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTa                                          <  2:1700846/150‑1 (MQ=255)
                                                                                                            ggttGAGTTTCGGCACTTCTTCGTGGAAGTCCGGCGCCTGATGCGCGGGGGTCAGACCGTTAGTCGCGGATCGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                         <  2:1509823/147‑1 (MQ=255)
                                                                                                            cGATGAGTTTCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                         <  2:277832/150‑1 (MQ=255)
                                                                                                            cGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAGGCTATCTTGAGAGTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTaa                                         >  1:1033417/1‑150 (MQ=255)
                                                                                                              aTGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAAc                                       >  1:896518/1‑150 (MQ=255)
                                                                                                                gAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTc                                     <  2:231736/150‑1 (MQ=255)
                                                                                                                   tttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGTGGCGCGGGGGTCAGCCCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGGCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCgg                                  <  1:1251819/150‑1 (MQ=255)
                                                                                                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGc                                 <  1:793631/150‑1 (MQ=255)
                                                                                                                    ttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGc                                 <  2:1100230/150‑1 (MQ=255)
                                                                                                                     tCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCg                                <  2:1693136/150‑1 (MQ=255)
                                                                                                                      cGCCACTTCTTCCTGGAAGTCCCCCCCCTGAGGCGCGGGGGTCAGACCGTTAGTAGCGGAACGACCAATGATATCTTTAGATTGATAAGCCCAAAGTCGTCACCGTTACACCTTTTACACAATAACAAAGAAGGACTTTAACTACCGCGt                               >  2:1387804/1‑150 (MQ=255)
                                                                                                                       gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTc                              <  2:1719684/150‑1 (MQ=255)
                                                                                                                        ccACTTCTTCCTGGAAGTCCCGCGCCGGGGGCGCGGGGGTGAGAACGGTTGTAGGGGGAACGGGAATTATATTGTTTTATTTATTAAGCATATGTTGTATCCCTTTAACCGTTTGCACATTAACAAGGAAGGGCGTTAACTcccgggggt                             >  2:1303417/1‑143 (MQ=255)
                                                                                                                           cttcttCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                          <  1:1420601/150‑1 (MQ=255)
                                                                                                                           cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                          <  1:67426/150‑1 (MQ=255)
                                                                                                                            ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                         <  1:1535974/150‑1 (MQ=255)
                                                                                                                             tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgtg                        <  1:1129615/150‑1 (MQ=255)
                                                                                                                               ttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCCGACCGTTATTCGCGGAACGCCGAATGATATCTTGAGATTGATAAGCCATATTTCGTCTCCGTTACACCGTTTCCACATTAACAAGGAAGGACTTAAAAAACGGCGTAAGCGTGtt                      >  2:841243/1‑150 (MQ=255)
                                                                                                                                 ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg                    <  2:1421560/150‑1 (MQ=255)
                                                                                                                                      gacccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCa               >  2:1131394/5‑150 (MQ=255)
                                                                                                                                        cccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             >  2:803221/3‑150 (MQ=255)
                                                                                                                                        cccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             >  1:1313609/3‑150 (MQ=255)
                                                                                                                                          cccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCaa           >  1:401531/1‑150 (MQ=255)
                                                                                                                                              cgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGaa       <  1:1387804/150‑1 (MQ=255)
                                                                                                                                                ccggAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     <  1:1342051/147‑1 (MQ=255)
                                                                                                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     >  1:332437/1‑150 (MQ=255)
                                                                                                                                                ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     >  1:853056/1‑150 (MQ=255)
                                                                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCg  <  1:70599/150‑1 (MQ=255)
                                                                                                                                                   gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCg  <  2:1712877/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  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: