Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,860,757 Δ1 bp 6.1% coding (222/654 nt) purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,860,7570G. 6.1% 731.3 / 40.0 214coding (222/654 nt)purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base G (95/106);  new base . (8/5);  total (103/111)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 226 total aligned reads displayed.

TGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGGCAGCT  >  NC_002947/1860608‑1860906
                                                                                                                                                      |                                                                                                                                                     
     gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCCGCGCTAATGGcgcgcatc                                                                                                                                                   <  1:385193/150‑5 (MQ=255)
       ccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCCGCGCTAATGGCGCGCATCGa                                                                                                                                                <  2:1868082/150‑1 (MQ=255)
        cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                               <  1:3974924/150‑1 (MQ=255)
        cTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGaa                                                                                                                                               >  2:4310129/1‑150 (MQ=255)
          tGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAgg                                                                                                                                             >  2:4043587/1‑150 (MQ=255)
           gTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGc                                                                                                                                            <  2:1263354/150‑1 (MQ=38)
            tACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCt                                                                                                                                           <  1:1950364/150‑1 (MQ=255)
             aCTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCtt                                                                                                                                          <  2:944826/150‑1 (MQ=38)
             aCTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGGCGATGCGCGCCATTAGCGCGGCATCGAAGGCtt                                                                                                                                          >  2:3789651/1‑150 (MQ=255)
              cTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCttt                                                                                                                                         <  1:2932073/150‑2 (MQ=38)
              cTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTc                                                                                                                                         <  2:537173/150‑1 (MQ=255)
               ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                                        >  1:2307486/1‑150 (MQ=255)
                  ggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACgg                                                                                                                                     <  1:4426914/150‑1 (MQ=255)
                  ggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACgg                                                                                                                                     <  1:4332187/150‑1 (MQ=255)
                   gcagcaGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGc                                                                                                                                    <  1:2571707/150‑1 (MQ=255)
                              gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGCACTgg                                                                                                                         >  1:3019779/1‑150 (MQ=255)
                              gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                                         <  1:2692011/150‑1 (MQ=255)
                                gTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                       >  2:1806315/1‑150 (MQ=255)
                                gTGGATGTTGAGCAGGCGGCCCTGATAGGGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGgt                                                                                                                       <  1:1158587/150‑1 (MQ=255)
                                       ttGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCg                                                                                                                >  1:2581712/1‑150 (MQ=255)
                                         gAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAg                                                                                                              >  1:2281351/1‑150 (MQ=255)
                                               gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                                        >  1:4049907/1‑150 (MQ=255)
                                                 ggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGc                                                                                                      <  1:1190365/150‑1 (MQ=255)
                                                    ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCCCGACGCTGt                                                                                                   >  1:771411/1‑150 (MQ=255)
                                                     cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                                  <  1:3472589/150‑1 (MQ=255)
                                                      tGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCg                                                                                                 <  2:3499702/150‑1 (MQ=255)
                                                      tGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCg                                                                                                 >  1:3320362/1‑150 (MQ=255)
                                                            tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                           >  1:3753884/1‑150 (MQ=255)
                                                             ggCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccgg                                                                                          <  1:448093/150‑1 (MQ=255)
                                                               cgcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggcc                                                                                        <  2:4205863/150‑1 (MQ=255)
                                                                gcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggccg                                                                                       <  1:1803168/150‑1 (MQ=255)
                                                                   cGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTg                                                                                    <  2:3254796/150‑1 (MQ=255)
                                                                          gccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGcgcg                                                                             >  1:1654882/1‑150 (MQ=255)
                                                                           ccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGt                                                                            >  2:2663940/1‑150 (MQ=255)
                                                                               tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                                        <  1:2236723/150‑1 (MQ=255)
                                                                               tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                                        >  2:4120439/1‑150 (MQ=255)
                                                                                 aGGATTCGCATGAAAACTGACAGCACCACCAGGTCCGGGGAAAATCATTCGATGATCTACATTAGCTCGTAATCGAAGGCTTCACGTCCTTATAACTGTTTGTGGTCGAGCACGACGCTGTCGATGCCGTCCGCTGAGTCGCGTTGCAgg                                                                      <  2:3287246/150‑1 (MQ=255)
                                                                                 aGGATGCGCATGACACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                      <  2:3169283/150‑1 (MQ=255)
                                                                                 aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                      >  2:1211706/1‑150 (MQ=255)
                                                                                 aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                      <  1:1496227/150‑1 (MQ=255)
                                                                                 aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                                      <  1:2854493/150‑1 (MQ=255)
                                                                                  ggATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAAGTGGGGGGGGTCGAGCACGACGCTGTCGATGCCGGCCGCGGCGGCGCGTTGCTGGc                                                                     >  2:2752213/1‑150 (MQ=255)
                                                                                   gATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGcc                                                                    <  1:2253673/150‑1 (MQ=255)
                                                                                   gATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGcc                                                                    >  1:2355457/1‑150 (MQ=255)
                                                                                      gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                                 <  1:514447/150‑1 (MQ=255)
                                                                                      gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                                 <  1:867653/150‑1 (MQ=255)
                                                                                      gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                                 <  2:43640/150‑1 (MQ=255)
                                                                                      gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                                 <  2:1061766/150‑1 (MQ=255)
                                                                                       cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATTGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                                <  1:1257411/150‑1 (MQ=255)
                                                                                       cgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAg                                                                <  2:3124570/150‑1 (MQ=255)
                                                                                        gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                               >  2:930830/1‑150 (MQ=255)
                                                                                        gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                               >  2:2032824/1‑150 (MQ=255)
                                                                                        gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                               >  2:1667790/1‑150 (MQ=255)
                                                                                         cATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGc                                                              >  1:2652014/1‑150 (MQ=255)
                                                                                         cATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGc                                                              >  1:2320754/1‑150 (MQ=255)
                                                                                          aTGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCa                                                             >  1:963375/1‑150 (MQ=255)
                                                                                           tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                            >  1:4239914/1‑150 (MQ=255)
                                                                                            gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                           >  2:2229959/1‑150 (MQ=255)
                                                                                              aaCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCAGAAGCGCGGCAGCGAAGGCATCACCGCCGGCGAACGGGGTGTGGTCGAGCACGACGCGGGCGATGCCGGCCGCAGCGGCGCGTTGCAGGCCCTAGGCCTCgg                                                         >  1:909324/1‑150 (MQ=255)
                                                                                               aCCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgc                                                        >  2:816449/1‑150 (MQ=255)
                                                                                                 cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgcgc                                                      <  2:1678906/150‑1 (MQ=255)
                                                                                                 cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGGCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCCTCGgcgc                                                      >  1:2861616/1‑150 (MQ=255)
                                                                                                  ggCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGcgcg                                                     >  2:133167/1‑150 (MQ=255)
                                                                                                  ggCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGcgcg                                                     <  2:3393043/150‑1 (MQ=255)
                                                                                                   gCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCgg                                                    >  2:2844907/1‑150 (MQ=255)
                                                                                                    ccAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                                   <  1:2072230/150‑1 (MQ=255)
                                                                                                     cAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCTCGGtt                                                  <  1:1442191/150‑1 (MQ=255)
                                                                                                       gCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTgg                                                >  2:4023709/1‑150 (MQ=255)
                                                                                                       gCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTgg                                                >  2:471572/1‑150 (MQ=255)
                                                                                                        caccacCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGa                                               <  2:723110/150‑1 (MQ=255)
                                                                                                         accaccAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGTCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaa                                              <  2:2521495/150‑1 (MQ=255)
                                                                                                         accaccAGGTCCGGCGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACCGGGTGTGGACGAGCACGCCGCTGTCGATGCCGGTCGCTGCGGCGCGTTGCACGCCGGAGGCAACGGCGCGGTTGGaa                                              <  2:263828/150‑1 (MQ=255)
                                                                                                          ccaccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaaa                                             <  2:4251804/150‑1 (MQ=255)
                                                                                                           caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                            >  2:4114157/1‑150 (MQ=255)
                                                                                                           caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                            <  1:3391968/150‑1 (MQ=255)
                                                                                                           caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                            <  1:2334361/150‑1 (MQ=255)
                                                                                                           caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                            <  1:3917644/150‑1 (MQ=255)
                                                                                                           caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                            >  2:2115699/1‑150 (MQ=255)
                                                                                                             ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                          <  1:4088102/150‑1 (MQ=255)
                                                                                                             ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                          >  2:1113875/1‑150 (MQ=255)
                                                                                                               aGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAcc                                        >  2:770592/1‑150 (MQ=255)
                                                                                                                   ccGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgc                                    >  2:1418380/1‑150 (MQ=255)
                                                                                                                   ccGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgc                                    <  2:2341868/150‑1 (MQ=255)
                                                                                                                             gCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc                          <  1:1281464/150‑1 (MQ=255)
                                                                                                                               cttCGATGCCGCGCTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                        >  1:277252/3‑150 (MQ=38)
                                                                                                                               cGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg                        <  2:2471926/150‑1 (MQ=255)
                                                                                                                                gTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg                       <  1:3143691/150‑1 (MQ=255)
                                                                                                                                 tCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCggg                      <  1:1724687/150‑1 (MQ=255)
                                                                                                                                   gATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCt                    >  1:3658968/1‑150 (MQ=255)
                                                                                                                                    aTGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCCTCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTg                   >  2:139301/1‑150 (MQ=255)
                                                                                                                                     tgccgcgcTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGt                  >  2:2741818/4‑150 (MQ=255)
                                                                                                                                        cgcgcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTcc                <  1:3126470/150‑1 (MQ=255)
                                                                                                                                           gcCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGc             >  2:3050460/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)
                                                                                                                                               ttAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCtt         <  2:3578984/150‑1 (MQ=255)
                                                                                                                                                 aGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTgg       <  1:471421/150‑1 (MQ=255)
                                                                                                                                                   cgcgGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGGCa     <  1:1915903/150‑1 (MQ=255)
                                                                                                                                                   cgcgGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGGCa     <  1:2386566/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCCCTTGGCAGCT  >  NC_002947/1860608‑1860906

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