Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,696 G→C 7.4% 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 7.4% 381.8 / 19.2 132Q22E (CAG→GAG) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (50/72);  new base C (3/7);  total (53/79)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-02

Only 100 of 134 total aligned reads displayed.

ATCATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCG  >  NC_000913/4082547‑4082845
                                                                                                                                                     |                                                                                                                                                     
  cATTGGGGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAg                                                                                                                                                     <  2:54918/150‑1 (MQ=255)
        gggcgTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgg                                                                                                                                               <  1:192085/146‑1 (MQ=255)
        ggTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTGCCTGGAAGTACCGCCGCGGAGGGGCgg                                                                                                                                               >  2:1386938/1‑150 (MQ=25)
          ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga                                                                                                                                             <  1:1640844/150‑2 (MQ=255)
          ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga                                                                                                                                             <  1:1834935/150‑2 (MQ=255)
          ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga                                                                                                                                             <  1:1359422/150‑2 (MQ=255)
           tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  1:1034902/150‑1 (MQ=255)
           tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  1:1843172/150‑1 (MQ=255)
           tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  2:1076308/150‑1 (MQ=255)
           tGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  1:1104771/150‑1 (MQ=255)
           tGTCGCACACCCCAATGTTATTGCCGACGGAATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCgggac                                                                                                                                            <  1:690925/150‑3 (MQ=255)
           ggccGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggggg                                                                                                                                            <  1:1245368/147‑1 (MQ=255)
              cGTACACCCCAATGTTATTGCCGACGGTATCGCTAATGTCGTTCCACTCTGACCACGCCACCTGCCAGGCTTTACAGCCGATACATGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCa                                                                                                                                         <  2:1062402/150‑1 (MQ=255)
                aacacCCCAATGTTCTTGCCGACGGTATCGCGACTGTCGTTCCCCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGACTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                                                                                                       <  1:1141874/149‑1 (MQ=37)
                  cacCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCGCCGGCGGGGCGGGGGGCAGAcc                                                                                                                                     >  1:1710539/1‑150 (MQ=25)
                    ccccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                                   >  2:1605448/1‑150 (MQ=255)
                    ccccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGt                                                                                                                                   >  2:284955/1‑150 (MQ=255)
                    ccccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCGCTGGGGCGCGGGGGTCAGACCGt                                                                                                                                   >  1:437216/1‑150 (MQ=255)
                       cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                                                                                                                <  1:1721827/150‑1 (MQ=255)
                       cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAg                                                                                                                                <  2:985935/150‑1 (MQ=255)
                           gTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGCCCGTTAGTcgc                                                                                                                            >  1:418789/1‑150 (MQ=255)
                            ttATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgcg                                                                                                                           <  2:862763/150‑1 (MQ=255)
                            ttATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTcgcg                                                                                                                           <  2:628979/150‑1 (MQ=255)
                                 gCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAacg                                                                                                                      >  2:1312181/1‑150 (MQ=255)
                                    gACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGCCCGTTAGTCGCGGCacgacg                                                                                                                   >  1:1581219/1‑150 (MQ=255)
                                     aCGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCGGGGGTCAGACCGTTAGTCGCGGAAcgacga                                                                                                                  <  1:1642681/150‑1 (MQ=255)
                                        ttATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATg                                                                                                               <  2:120800/149‑1 (MQ=255)
                                        gTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATg                                                                                                               <  1:1095617/150‑1 (MQ=255)
                                         tATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCGTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATgc                                                                                                              >  1:56258/1‑149 (MQ=255)
                                            cgcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGCACGACGAATGagat                                                                                                           >  2:1133735/1‑150 (MQ=255)
                                             gcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGGTAGTCGCGGAACGACGAATGATATc                                                                                                          >  1:154568/1‑150 (MQ=255)
                                              cgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCCCCCCTGAGGCGCGGGGGTCAGCCCGTTTGTCGCGGAACGACGCATCATATCt                                                                                                         >  2:1469358/1‑150 (MQ=25)
                                               taaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTAAAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                        <  2:1569975/149‑1 (MQ=255)
                                                aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTg                                                                                                       <  2:102133/150‑1 (MQ=255)
                                                aaTGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTg                                                                                                       >  2:705173/1‑150 (MQ=255)
                                                       ttCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACCGTTAGTCGCGGCACGCCGAATGATATCTTGAGATTGa                                                                                                >  1:411130/1‑150 (MQ=255)
                                                       ttCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACCGTTAGTCGCGGAACGACGAATGATATGTTGAGATTgt                                                                                                >  1:500973/1‑149 (MQ=255)
                                                       ttCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGa                                                                                                >  2:1820767/1‑150 (MQ=255)
                                                        tCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGAACGACGAATGATATCTTGAGTTTGAt                                                                                               >  2:251448/1‑150 (MQ=255)
                                                         ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              <  2:864456/150‑1 (MQ=255)
                                                         ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              <  2:695918/150‑1 (MQ=255)
                                                         ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATa                                                                                              <  1:1276196/150‑1 (MQ=255)
                                                         ccACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTGCTTCCTGGAAGTCCCGCGCCTGAGGCTCGGGGGTCAGACCGTTAGGCGCGGAACGACGCAGGAAATCTTGAGATTGATa                                                                                              >  1:309925/1‑150 (MQ=37)
                                                            ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGGGGGGGGCAGCCCGGTAGTCGCGGGACGCCGAATGATTTGTTGTGGTTGCTaatg                                                                                           >  1:1570426/1‑148 (MQ=25)
                                                               tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGCGGCGCCGGGGGCAGACCGTTAGTCGCGGCACGACGCATGATATATTGAGGTTGATACGCCat                                                                                        >  1:509377/1‑150 (MQ=37)
                                                               tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCat                                                                                        <  2:114513/150‑1 (MQ=255)
                                                                gAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCata                                                                                       <  1:22567/150‑1 (MQ=255)
                                                                  acacGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATg                                                                                     >  1:1197526/1‑150 (MQ=255)
                                                                   cacGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgt                                                                                    >  2:957432/1‑150 (MQ=255)
                                                                     cGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCGGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtcg                                                                                  >  2:1743214/1‑150 (MQ=255)
                                                                       ccACCTGACAGGCTATACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCCCCTGAAGGCCCGGGGTCAGACCCGTAATAGCGGAACGACGAAAGATATGTTGAGATTGTTAAGCCataattagag                                                                                >  2:1655453/1‑143 (MQ=255)
                                                                         aCCTGACAGGCTTTTCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtctc                                                                              <  2:461654/150‑1 (MQ=255)
                                                                           cTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCg                                                                            <  1:1756869/150‑1 (MQ=255)
                                                                             gACAGGCTTTACAGCCGATACTGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGtt                                                                          <  1:982099/150‑1 (MQ=255)
                                                                             gACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGtt                                                                          <  2:1470819/150‑1 (MQ=255)
                                                                             gACAGGCTATACAGCCGATACAGGTGGTAACGCCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGtt                                                                          <  1:1663109/150‑1 (MQ=255)
                                                                               cAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCAACTTCTTCCTGGAACTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGTATGAAATCTTGAGATTGATAAGCCATATGTCGTCAACGTTaa                                                                        >  2:1129043/1‑149 (MQ=38)
                                                                                aGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTaca                                                                       <  1:21018/150‑1 (MQ=255)
                                                                                 ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTacac                                                                      >  1:1287665/1‑150 (MQ=255)
                                                                                 ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTACGCCCGTTTGTCGCGGAACGACGCATGATAACTTGAGATTGATAAGCCATATGTAGTCGCCCTTacac                                                                      >  2:831375/1‑150 (MQ=38)
                                                                                 ggCTTTACAGCCGATACAGGTGGTAACGTCGATGAGATTCGCCACTTCTTCCTGGAAGACCCCCGCCACACGCGCGGGGGTCCGACCCTTAGTCGCGGAACCACGAAAAATATCTTGACAATGAAAAGCCAAATGTCGTCTcccataccc                                                                      >  2:658480/1‑143 (MQ=255)
                                                                                     ttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGCCCGTTAGTCGCGGAACGACGAATGATATCGTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCttt                                                                  >  2:146324/1‑150 (MQ=255)
                                                                                       aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTc                                                                <  1:1490429/150‑1 (MQ=255)
                                                                                       aCAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTc                                                                <  1:1312181/150‑1 (MQ=255)
                                                                                        cAGCCGATACAGGTGGTAACGTCGATGAGGTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTcc                                                               <  2:1726344/150‑1 (MQ=255)
                                                                                        cAGCCGATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATCAGCCATATGTCGTCTCCGTTACACCTTTTcc                                                               <  2:1200151/150‑1 (MQ=255)
                                                                                         aGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCca                                                              <  2:1430791/150‑1 (MQ=255)
                                                                                           ccGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAGCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCcaca                                                            <  2:1565728/150‑1 (MQ=255)
                                                                                             gATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTATTCCACAtt                                                          >  1:160414/1‑150 (MQ=255)
                                                                                              aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGACGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc                                                         <  1:644387/150‑1 (MQ=255)
                                                                                              aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTc                                                         <  2:1293652/150‑1 (MQ=255)
                                                                                               tACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTca                                                        <  2:1589415/150‑1 (MQ=255)
                                                                                                 cAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTcaca                                                      <  2:614835/150‑1 (MQ=255)
                                                                                                   ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag                                                    <  2:436468/150‑1 (MQ=255)
                                                                                                   ggtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaag                                                    <  1:834093/150‑1 (MQ=255)
                                                                                                    gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGCACGACGAATGATATCTTGACATTGATAAGCCATAAGTCGTCTCCGTTACACCGTTTCCACATTCACaagg                                                   >  2:13792/1‑150 (MQ=255)
                                                                                                       gtAACGTCGATGAGTTTCGCCTCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaaggaag                                                <  1:159548/150‑1 (MQ=255)
                                                                                                            gTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGTTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACtt                                           <  2:1600246/150‑1 (MQ=255)
                                                                                                                   aGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTcc                                    >  2:1087792/1‑150 (MQ=255)
                                                                                                                     tttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCgg                                  <  2:252665/150‑1 (MQ=255)
                                                                                                                       tCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCg                                <  1:1874434/150‑1 (MQ=255)
                                                                                                                        cGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATTTCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCTCATTCACAAGGAAGGACTTAAACTCCGGCGt                               >  2:1227263/1‑150 (MQ=255)
                                                                                                                         gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTc                              <  1:1323836/150‑1 (MQ=255)
                                                                                                                         gCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTc                              <  2:475900/150‑1 (MQ=255)
                                                                                                                            aCTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGc                           >  1:27168/1‑150 (MQ=255)
                                                                                                                             cttcttCCTGGAAGTCCCGCGACTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCg                          <  1:1324482/150‑1 (MQ=255)
                                                                                                                               tcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgtg                        <  1:1005108/150‑1 (MQ=255)
                                                                                                                                   ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg                    <  2:1082758/150‑1 (MQ=255)
                                                                                                                                   ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCg                    <  1:1389574/150‑1 (MQ=255)
                                                                                                                                    cTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGATACACCTTTTCCACATTCACAAGGAAGGACTTAAAATCCGGCGTCTGCGTGTTCGc                   >  1:1288042/1‑150 (MQ=255)
                                                                                                                                    cTGGAAGTCCCGCGCCAGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATGTTGAGATTGAAAAGCCATATGTCGTCTCCGTTAAACATTTACCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGc                   >  2:135177/1‑150 (MQ=255)
                                                                                                                                     tGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCa                  >  1:1050960/1‑150 (MQ=255)
                                                                                                                                     tGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCa                  <  1:1552691/150‑1 (MQ=255)
                                                                                                                                          gTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             <  2:309925/150‑1 (MQ=255)
                                                                                                                                          gTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTACACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCAcc             >  1:826708/1‑150 (MQ=255)
                                                                                                                                           ccccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGCCGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCa            >  1:1325401/2‑150 (MQ=255)
                                                                                                                                             ccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAAc          <  1:284955/150‑1 (MQ=255)
                                                                                                                                                  ccTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCGACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATg     >  2:349168/1‑150 (MQ=255)
                                                                                                                                                    tGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGc   <  2:631456/150‑1 (MQ=255)
                                                                                                                                                     gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACGAATGGCg  <  1:705173/150‑1 (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: