Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,696 G→C 5.0% 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 5.0% 483.2 / 13.1 154Q22E (CAG→GAG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (60/84);  new base C (4/4);  total (65/89)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 157 total aligned reads displayed.

ATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  NC_000913/4082547‑4082845
                                                                                                                                                     |                                                                                                                                                     
aTCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATAGCGAATGTCGTTCCACTCTTAACAAGCCACCTGACAGGCTTTACAGCCGATACAGGTAGTAACGTCGATTAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCTg                                                                                                                                                       >  2:1855245/1‑150 (MQ=37)
 tCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGa                                                                                                                                                      <  2:1336565/150‑1 (MQ=255)
   aTTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAgg                                                                                                                                                    >  2:658137/1‑150 (MQ=255)
      ggggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCCCTTCTTCCTGGAAGTCCCGCGCCTGAGgcgc                                                                                                                                                 >  2:279218/1‑150 (MQ=255)
        ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                                                                                                               <  1:2084133/150‑1 (MQ=255)
        ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                                                                                                               >  1:1933795/1‑150 (MQ=255)
         tttGTAGTACACCCAAATGTTATTGCCGACGTTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                              <  2:699360/149‑1 (MQ=255)
         gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                              <  2:2059648/150‑1 (MQ=255)
         gTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                              <  1:437867/150‑1 (MQ=255)
           tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  2:601607/150‑1 (MQ=255)
           tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  2:1465485/150‑1 (MQ=255)
              cGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCAGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGCGGTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGGGcggcgggcc                                                                                                                                         >  2:1978248/1‑147 (MQ=37)
               gTACACCCCAATGTAATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAg                                                                                                                                        <  2:1892921/150‑1 (MQ=255)
                tACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAgc                                                                                                                                       >  1:19385/1‑149 (MQ=255)
                   ccccaAATGTTATTGCCGACGGTATCGCGAATTGCGTTCCACTCTCAACACGCCCCCTGACAGGCTTTACAGCCGATACAGGTGGAAACGTCTATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCg                                                                                                                                    <  1:1879347/149‑1 (MQ=37)
                    cccaAATGTTATTGCCGACGGTTTCGCGAATGTCGTGCCACTCTCAACACGCCCCCTGACAGGCTTTACCGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGt                                                                                                                                   <  2:1821434/150‑1 (MQ=255)
                     cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCTCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGtt                                                                                                                                  <  1:45792/150‑1 (MQ=255)
                     cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCAGCAAGTCCCCCGCCTCACGCGCGGCGGTCAGaaccgt                                                                                                                                  >  1:485942/1‑145 (MQ=25)
                       cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                                                                                                                <  2:1960552/150‑1 (MQ=255)
                        aaTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                                                                                                               >  1:19020/1‑150 (MQ=255)
                          tGTTATTGCCGACTGTATCGAGAATGTCGTTCAACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcg                                                                                                                             <  2:908978/150‑1 (MQ=255)
                           gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgc                                                                                                                            >  1:321566/1‑150 (MQ=255)
                           gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCTGAGGCCCGGGGGTCAGACCCTTAGTcgc                                                                                                                            >  2:640888/1‑150 (MQ=255)
                            ttATTGCCGACGGTATCGCGAAGGTCGTTCCACCCTCAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgcg                                                                                                                           <  1:231037/150‑1 (MQ=255)
                              aTTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGa                                                                                                                         <  2:886918/150‑1 (MQ=255)
                               ttGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGGCAGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGaa                                                                                                                        >  1:1174920/1‑150 (MQ=255)
                                tGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAac                                                                                                                       >  1:869682/1‑150 (MQ=255)
                                tGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCACACCCGTAGGCGCGGCac                                                                                                                       >  2:1127942/1‑150 (MQ=25)
                                   cGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACTGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgac                                                                                                                    <  2:951868/150‑1 (MQ=255)
                                    gACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgacg                                                                                                                   <  2:1627993/150‑1 (MQ=255)
                                    gACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacgacg                                                                                                                   <  1:1515032/150‑1 (MQ=255)
                                    gACGGGATCGCGAATGTCGTTCCACTCTGAGCACGCCACCTGCCAGGCTTTACAGCCGATACAGGTGGTCACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCGGGAGGCGCGCGGGTCAGCCCGGTACTAGCGGAccgccc                                                                                                                   >  2:2145618/1‑147 (MQ=21)
                                       ggTATCGCGAATGGCGTTCCCCTCTGACCACGCCACCTGACCGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTGCGCCACTGCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAAt                                                                                                                <  1:2028918/150‑1 (MQ=38)
                                         tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                                                                                                              <  1:1521673/150‑1 (MQ=255)
                                         tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                                                                                                              >  1:1295411/1‑150 (MQ=255)
                                         tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGa                                                                                                              >  2:1233351/1‑150 (MQ=255)
                                           tCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGata                                                                                                            <  2:248887/150‑1 (MQ=255)
                                            cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                                                                                                           >  2:1789103/1‑150 (MQ=255)
                                            cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGAAGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGatat                                                                                                           >  2:14974/1‑150 (MQ=255)
                                                aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCAATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTg                                                                                                       <  2:1135313/150‑1 (MQ=255)
                                                 aTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTga                                                                                                      >  2:65476/1‑150 (MQ=255)
                                                  tGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgag                                                                                                     >  1:1251280/1‑150 (MQ=255)
                                                    gcGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAt                                                                                                   <  1:2094489/149‑1 (MQ=255)
                                                     cGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGGGGCGCGGGGGTCAGACCGGTAGTAGCGGGACGACGCATGCTATCTTGAGAtt                                                                                                  >  1:2171801/1‑150 (MQ=37)
                                                     cGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                                  >  1:1570880/1‑150 (MQ=255)
                                                      gTTCCACTCTGAACACGCCAGCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg                                                                                                 <  2:104057/150‑1 (MQ=255)
                                                         ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              >  2:839678/1‑150 (MQ=255)
                                                         ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              <  1:818661/150‑1 (MQ=255)
                                                            ctctGAACACGCCACCTGACAGGCTGTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCTCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGc                                                                                           <  1:1985243/150‑1 (MQ=255)
                                                             tctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGcc                                                                                          <  1:691727/150‑1 (MQ=255)
                                                              ctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGCTGTCACGCGACTGAGGCGCTGGGGTCCGCCCGGTTGTCGCGGGGCGCCGAGGGATTTGTTGCGATTGATAGGcct                                                                                         >  2:363752/1‑149 (MQ=18)
                                                               tGAACACGCCACCTGCCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCat                                                                                        <  2:1404074/150‑1 (MQ=255)
                                                                gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata                                                                                       <  1:1974072/150‑1 (MQ=255)
                                                                 aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCatat                                                                                      >  2:348374/1‑150 (MQ=255)
                                                                 aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGAGGCGCGGGGGGGAGACCGTTTGTGGCGGGACGCCGAATGATATGTTGAGATTGATAAGCCatat                                                                                      >  1:485627/1‑150 (MQ=37)
                                                                   cacGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGGGGCGCGGGGGTGAGACCGGTTGTGGCGGGCCGCCGGATGATTTGTGGTGTTTGATAAGCCCTATgt                                                                                    >  1:788342/1‑150 (MQ=21)
                                                                     cGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcg                                                                                  >  2:125763/1‑150 (MQ=255)
                                                                     cGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTGAGGCCGGTTGTCGCGGGACGCCGAATGATATCTTGAGATTGATAAGCCATATgtcg                                                                                  >  1:2201693/1‑150 (MQ=38)
                                                                          ccTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTcc                                                                             <  2:627278/150‑1 (MQ=255)
                                                                           cTGAGAGGCTTTACAGCCTATTCCGGTGGGAACGTCGATGAGTTTCTCCACTTCTTCCTGGAAGTCCCGCGCCTGATGCGCGGGGGTCAGACCGTTAGGCGCGGAACGACGAATGATATCTTGAGCTTGATAAGCCATATGTCGTCTCCg                                                                            <  2:2112386/150‑1 (MQ=255)
                                                                             gACAGGCTTTAGAGCCGATACAGGTGGTACCGGCGATGCGTTTCGCCAATTCTTCCTGGAAGTCCAGCGCCTGAGGCGCGGGGGTCAGACCGGTAGGCGCGGAACGACGAATGATATCTTGAGAGTGATAAGCCATATGTCGTCTCCGtt                                                                          <  1:437805/150‑1 (MQ=255)
                                                                              aCAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTa                                                                         <  2:79298/150‑1 (MQ=255)
                                                                               cAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGAGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTATCCGTTaa                                                                        >  2:1302025/1‑149 (MQ=255)
                                                                                 ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTAGCCACTTCTTCCTGGAAGTCCCGCGCCAGCGGCGCGGGGGTCAGACCGGTAGTCGCGGAACGTCGCATGATATCTTGAGATTGATTAGCCATATGTAGTCTCCGTTagac                                                                      >  2:1790507/1‑150 (MQ=38)
                                                                                   cTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCt                                                                    <  1:1299851/150‑1 (MQ=255)
                                                                                     ttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTTGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCttt                                                                  <  1:260764/150‑1 (MQ=255)
                                                                                      tACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt                                                                 >  2:1542484/1‑150 (MQ=255)
                                                                                      tACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtttt                                                                 >  2:434290/1‑150 (MQ=255)
                                                                                        cAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGCGGCGCGGGGGGGAGACCGGTTGTCGCGGGACGACGAATGGTTTATTGTGATTGGTAAGCCATATTTCGTCTCCCGTAGACCttttga                                                               >  2:763568/1‑148 (MQ=255)
                                                                                          gCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcac                                                             <  2:823648/150‑1 (MQ=255)
                                                                                            cGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCACAtt                                                          <  2:758959/150‑1 (MQ=255)
                                                                                             gATAGAGGTGGTAAGGTCGATGAGTTTCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACAACGAAAGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt                                                          <  1:57405/150‑1 (MQ=255)
                                                                                              aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc                                                         <  1:1717663/150‑1 (MQ=255)
                                                                                               tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTca                                                        <  2:775466/150‑1 (MQ=255)
                                                                                                   ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag                                                    <  1:1132403/150‑1 (MQ=255)
                                                                                                        tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaagg                                               <  1:640888/150‑1 (MQ=255)
                                                                                                         aaCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaaggc                                              >  1:1138540/1‑149 (MQ=255)
                                                                                                          aCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGGCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGAc                                             <  1:1089967/150‑1 (MQ=255)
                                                                                                           cGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACt                                            >  2:1223000/1‑150 (MQ=255)
                                                                                                           cGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTATCCCTTACACCTTTTCCACATTAACAAGGAAGGAca                                            >  2:562178/1‑149 (MQ=255)
                                                                                                                aTGAGTTTCGACACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACAGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAAc                                       >  1:351323/1‑150 (MQ=255)
                                                                                                                    gTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCg                                   >  2:1213430/1‑150 (MQ=255)
                                                                                                                         gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTc                              >  2:5340/1‑150 (MQ=255)
                                                                                                                           cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATTTCGTCTCCGTTACACCTTTTACACATTCACAAGGAAGGACGTAAACTACGGCGTCTg                            >  1:636120/1‑150 (MQ=255)
                                                                                                                           cACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTg                            <  1:754854/150‑1 (MQ=255)
                                                                                                                             cttcttCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                          <  1:1026376/150‑1 (MQ=255)
                                                                                                                              ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                         <  1:830292/150‑1 (MQ=255)
                                                                                                                               tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGATTTAAACTCCGGCGTCTGCgtg                        <  2:321566/150‑1 (MQ=255)
                                                                                                                                 ttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGtt                      <  1:125763/150‑1 (MQ=255)
                                                                                                                                 ttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGtt                      <  2:796085/150‑1 (MQ=255)
                                                                                                                                   ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg                    <  1:1789103/150‑1 (MQ=255)
                                                                                                                                    cTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGc                   >  1:1524772/1‑150 (MQ=255)
                                                                                                                                          gacccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGGCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             <  2:130375/148‑1 (MQ=255)
                                                                                                                                           tCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCa            >  1:1030384/1‑150 (MQ=255)
                                                                                                                                           ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGCCGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCa            >  1:1808146/2‑150 (MQ=255)
                                                                                                                                            cccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCaa           >  2:445620/1‑150 (MQ=255)
                                                                                                                                             ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCCCCAGGAAGGCCTTAAACTCCGGCGTCTGCGTGTTCGCATCACCacc          >  1:207346/1‑148 (MQ=255)
                                                                                                                                               gcgcCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGa        <  2:376998/150‑1 (MQ=255)
                                                                                                                                                     gAGGCGCGGGGGTCAGACCGTTAGTCTCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGTACTTAAAATCATGCGTCTGAGTGTTCGCATAAACAAAGAAtttct  >  1:1559278/1‑146 (MQ=255)
                                                                                                                                                     gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCg  >  2:708761/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  NC_000913/4082547‑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: