Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,860,752 A→G 6.6% L76P (CTA→CCA)  purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,860,7520AG 6.6% 671.4 / 23.9 212L76P (CTA→CCA) purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base A (94/103);  new base G (8/6);  total (102/110)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 217 total aligned reads displayed.

GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGG  >  NC_002947/1860603‑1860901
                                                                                                                                                      |                                                                                                                                                     
gCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTa                                                                                                                                                       >  2:1281305/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                             <  2:152041/150‑1 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                             >  1:2245731/1‑150 (MQ=255)
          gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                             <  2:1766154/150‑1 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          <  1:3974924/150‑1 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          <  1:239812/150‑1 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          >  2:1226085/1‑150 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          <  2:2828483/150‑1 (MQ=255)
             cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                          >  2:4310129/1‑150 (MQ=255)
               tGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACgg                                                                                                                                        <  1:84132/150‑1 (MQ=38)
                  aCTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCtt                                                                                                                                     <  2:944826/150‑1 (MQ=38)
                  aCTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGGCGATGCGCGCCATTAGCGCGGCATCGAAGGCtt                                                                                                                                     >  2:3789651/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   >  2:1412519/1‑150 (MQ=255)
                    ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                   >  1:2307487/1‑150 (MQ=255)
                      gggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACg                                                                                                                                 >  2:3960373/1‑150 (MQ=255)
                       ggCAGCAGCGACGGGTGGATGTTGTGCAGGCGGCCCTGTTAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACgg                                                                                                                                <  1:2895399/150‑1 (MQ=255)
                        gcagcaGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGc                                                                                                                               <  1:2571707/150‑1 (MQ=255)
                              gcGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGa                                                                                                                         <  2:2626969/150‑1 (MQ=255)
                               cGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGaa                                                                                                                        >  2:45390/1‑150 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    >  1:4384162/1‑150 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  1:961680/150‑1 (MQ=255)
                                   gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                    <  1:2692011/150‑1 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgggg                                                                                                                  >  2:1806337/1‑149 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                  <  1:905473/150‑1 (MQ=255)
                                     gTGGATGTTGAGCAGGCGGCCCTGATAGGGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                  <  1:1158587/150‑1 (MQ=255)
                                          tGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGt                                                                                                             <  1:2142677/150‑1 (MQ=255)
                                           gTTGAGCAGGCGGCCCTTATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCTGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                            <  1:3127173/150‑1 (MQ=255)
                                            ttGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCg                                                                                                           >  1:2581712/1‑150 (MQ=255)
                                                      ggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGc                                                                                                 <  1:1190365/150‑1 (MQ=255)
                                                       gcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCt                                                                                                >  1:4443545/1‑150 (MQ=255)
                                                         ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                              <  2:2854357/150‑1 (MQ=255)
                                                          cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                             >  2:3529215/1‑150 (MQ=255)
                                                           tGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCg                                                                                            >  1:3320187/1‑150 (MQ=255)
                                                               aGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgc                                                                                        >  1:3885615/1‑150 (MQ=255)
                                                                 tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCCGTAGCGCGGCATCGCAGGCTTCCCGGCCGTCGAACTGGGGGGGGTCGCGCCCGACGCTGTCGATgccg                                                                                      >  2:49851/1‑150 (MQ=255)
                                                                 tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                      >  1:3753884/1‑150 (MQ=255)
                                                                    cgcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggcc                                                                                   <  2:4205863/150‑1 (MQ=255)
                                                                    cgcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggcc                                                                                   <  2:3479724/150‑1 (MQ=255)
                                                                     gcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggccg                                                                                  <  1:1803168/150‑1 (MQ=255)
                                                                       aCGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCt                                                                                >  1:493319/1‑150 (MQ=255)
                                                                        cGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                                               <  2:3254796/150‑1 (MQ=255)
                                                                            gccgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgcggc                                                                           >  1:1695277/1‑150 (MQ=255)
                                                                                ccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGt                                                                       >  2:2663940/1‑150 (MQ=255)
                                                                                  aCTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTg                                                                     >  2:742948/1‑150 (MQ=255)
                                                                                    tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGATGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                                   <  1:553291/150‑1 (MQ=255)
                                                                                    tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                                   >  2:4120439/1‑150 (MQ=255)
                                                                                    tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                                   <  2:2352995/150‑1 (MQ=255)
                                                                                      aGGATTCGCATGAAAACTGACAGCACCACCAGGTCCGGGGAAAATCATTCGATGATCTACATTAGCTCGTAATCGAAGGCTTCACGTCCTTATAACTGTTTGTGGTCGAGCACGACGCTGTCGATGCCGTCCGCTGAGTCGCGTTGCAgg                                                                 <  2:3287246/150‑1 (MQ=255)
                                                                                      aGGATGCGCATGACACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 <  2:3169283/150‑1 (MQ=255)
                                                                                      aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 <  1:2854493/150‑1 (MQ=255)
                                                                                      aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 <  1:1743675/150‑1 (MQ=255)
                                                                                      aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                 <  1:1496227/150‑1 (MQ=255)
                                                                                       ggATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAAGTGGGGGGGGTCGAGCACGACGCTGTCGATGCCGGCCGCGGCGGCGCGTTGCTGGc                                                                >  2:2752213/1‑150 (MQ=255)
                                                                                           gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                            <  2:1284737/150‑1 (MQ=255)
                                                                                           gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                            <  2:43640/150‑1 (MQ=255)
                                                                                           gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                            <  1:2268628/150‑1 (MQ=255)
                                                                                            cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATTGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                           <  1:1257411/150‑1 (MQ=255)
                                                                                            cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                           <  2:3124570/150‑1 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                          >  2:1667790/1‑150 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                          >  2:2032824/1‑150 (MQ=255)
                                                                                             gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                          >  2:930830/1‑150 (MQ=255)
                                                                                              cATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGc                                                         >  1:2320754/1‑150 (MQ=255)
                                                                                              cATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGc                                                         >  1:2652014/1‑150 (MQ=255)
                                                                                               aTGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCa                                                        >  2:1931309/1‑150 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                       >  1:4239914/1‑150 (MQ=255)
                                                                                                tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                       >  1:2300062/1‑150 (MQ=255)
                                                                                                 gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                      >  2:2229959/1‑150 (MQ=255)
                                                                                                   aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCgg                                                    >  1:1968140/1‑150 (MQ=255)
                                                                                                    aCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgc                                                   >  2:816449/1‑150 (MQ=255)
                                                                                                      cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGGCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCCTCGgcgc                                                 >  1:2861616/1‑150 (MQ=255)
                                                                                                       ggCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGcgcg                                                <  2:3393043/150‑1 (MQ=255)
                                                                                                       ggCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGcgcg                                                >  2:1162631/1‑150 (MQ=255)
                                                                                                          cAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCTCGGtt                                             <  1:1442191/150‑1 (MQ=255)
                                                                                                           aGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                                            >  1:3365493/1‑150 (MQ=255)
                                                                                                              accaccAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGTCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaa                                         <  2:2521495/150‑1 (MQ=255)
                                                                                                               ccaccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaaa                                        <  2:4251804/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                       <  1:2334361/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                       <  2:1111519/150‑1 (MQ=255)
                                                                                                                caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                       >  2:2115699/1‑150 (MQ=255)
                                                                                                                  ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                     >  1:995547/1‑150 (MQ=255)
                                                                                                                  ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                     <  1:4088037/150‑1 (MQ=255)
                                                                                                                        ccGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgc                               <  1:2941737/150‑1 (MQ=255)
                                                                                                                          ggggCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGACCTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCgg                             >  2:4050784/1‑150 (MQ=255)
                                                                                                                                  gCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc                     <  1:1281464/150‑1 (MQ=255)
                                                                                                                                    cttCGATGCCGCGCTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                   >  1:277197/3‑150 (MQ=38)
                                                                                                                                    cGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                   <  2:2471926/150‑1 (MQ=255)
                                                                                                                                     ttCGATGCCGCGCTAATGGCGCGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg                  >  2:4036275/2‑150 (MQ=255)
                                                                                                                                      tCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                 >  1:3578025/1‑150 (MQ=255)
                                                                                                                                        gatgccgcgcTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCt               >  1:3431868/6‑150 (MQ=255)
                                                                                                                                        gATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCt               >  1:3659032/1‑150 (MQ=255)
                                                                                                                                         aTGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTg              >  2:139236/1‑150 (MQ=255)
                                                                                                                                          tgccgcgcTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGt             >  2:2741818/4‑150 (MQ=255)
                                                                                                                                             cgcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTcc           <  1:3126406/150‑1 (MQ=255)
                                                                                                                                              gcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCt          <  1:45390/150‑1 (MQ=255)
                                                                                                                                                gcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGc        >  2:781760/1‑150 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc       <  2:2580844/150‑1 (MQ=255)
                                                                                                                                                 ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc       <  2:2354631/150‑1 (MQ=255)
                                                                                                                                                  cATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGccc      <  2:3699396/150‑1 (MQ=255)
                                                                                                                                                      aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg  <  1:4043587/150‑1 (MQ=255)
                                                                                                                                                      aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg  <  1:1925667/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCGTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGG  >  NC_002947/1860603‑1860901

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