Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 937,771 A→T 7.8% A2618A (GCT→GCA)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947937,7710AT 7.8% 280.5 / 11.1 90A2618A (GCT→GCA) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (47/36);  new base T (3/4);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GCTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATGCTGCT  >  NC_002947/937627‑937920
                                                                                                                                                |                                                                                                                                                     
gCTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCgg                                                                                                                                                  <  1:360633/150‑1 (MQ=255)
 cTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGTCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCggg                                                                                                                                                 <  1:2105734/150‑1 (MQ=255)
 cTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCggg                                                                                                                                                 <  1:2150318/150‑1 (MQ=255)
   gATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtc                                                                                                                                               <  1:1315157/150‑1 (MQ=255)
     tgctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTCCCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcga                                                                                                                                             <  1:958535/150‑1 (MQ=255)
     tgctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGTTCCCGCTGTCACCCGGCAGCGCAGCCGggtcga                                                                                                                                             <  1:64391/150‑1 (MQ=255)
     tgctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggtgaa                                                                                                                                             <  1:423572/150‑4 (MQ=255)
      gctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggtgaac                                                                                                                                            <  2:653467/150‑5 (MQ=255)
       ctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgagg                                                                                                                                           <  2:109538/150‑1 (MQ=255)
        tgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgaggt                                                                                                                                          <  2:622879/150‑1 (MQ=255)
        tgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgaggt                                                                                                                                          >  1:1941591/1‑150 (MQ=255)
                gcGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAgg                                                                                                                                  <  2:1264088/149‑1 (MQ=255)
                   gCAGCGACGACCTGTTCGCTGGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCCGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATGAGGaaa                                                                                                                               <  2:41523/150‑1 (MQ=255)
                   gCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGaaa                                                                                                                               <  1:2236311/150‑1 (MQ=255)
                        gacgacCGGTTCGCTCGGGTTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTAcc                                                                                                                          <  1:568968/150‑1 (MQ=255)
                        gacgacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCCGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTAcc                                                                                                                          >  1:980913/1‑150 (MQ=255)
                        gacgacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTAcc                                                                                                                          >  2:1862364/1‑150 (MQ=255)
                           gacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggtgaa                                                                                                                                             >  1:1904288/1‑125 (MQ=255)
                           gacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggtgaa                                                                                                                                             <  2:1904288/128‑4 (MQ=255)
                           gacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGaaa                                                                                                                               >  2:2328048/1‑142 (MQ=255)
                           gacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGaaa                                                                                                                               <  1:2328048/142‑1 (MQ=255)
                           gacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGtt                                                                                                                       >  1:990507/1‑150 (MQ=255)
                             ccAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTcg                                                                                                                     <  2:1586457/150‑1 (MQ=255)
                             ccAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCGCGGGGCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGACGATTAGGAAAGTACCGTTcg                                                                                                                     >  1:1931298/1‑150 (MQ=255)
                              cAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTcgc                                                                                                                    >  2:1976331/1‑150 (MQ=255)
                              cAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCAACCGGCAGCGCAGCCGTGTCGAGGTCGATCAGGAAAGTACCGTTcgc                                                                                                                    >  2:396112/1‑150 (MQ=255)
                                gTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCgg                                                                                                                  <  2:1081309/150‑1 (MQ=255)
                                     cTCGGGCTGGTCTGTGTGACCACCAGTGGCACTTCGTAGGTCATCGCTTCAGAGGCCTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACt                                                                                                             <  2:990507/150‑1 (MQ=255)
                                         ggCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACtgtg                                                                                                          >  2:644781/1‑149 (MQ=255)
                                         ggCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACtgtg                                                                                                          <  1:644781/149‑1 (MQ=255)
                                             ggTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGt                                                                                                     <  2:1323526/150‑1 (MQ=255)
                                               tCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTAc                                                                                                   >  2:1559737/1‑150 (MQ=255)
                                                cGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGACTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTAcc                                                                                                  <  2:2238191/150‑1 (MQ=255)
                                                  gCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGa                                                                                                >  2:1029579/1‑150 (MQ=255)
                                                     tGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCa                                                                                             <  1:1976331/150‑1 (MQ=255)
                                                       accaccAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGc                                                                                           <  1:177770/150‑1 (MQ=255)
                                                         caccaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGt                                                                                         >  1:1316493/1‑150 (MQ=255)
                                                         caccaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGACGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTATGCCGGTACCGACCAGCGt                                                                                         >  1:859819/1‑150 (MQ=255)
                                                          accaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCgg                                                                                                                                                  >  1:1591943/1‑92 (MQ=255)
                                                          accaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCgg                                                                                                                                                  <  2:1591943/92‑1 (MQ=255)
                                                           ccaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGc                                                                                       <  1:1263848/150‑1 (MQ=255)
                                                              gTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCAt                                                                                    >  1:1512880/1‑150 (MQ=255)
                                                               cggCACGTAGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCAtt                                                                                   <  1:865472/149‑1 (MQ=255)
                                                                 gCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTgg                                                                                 >  2:2341373/1‑150 (MQ=255)
                                                                 gCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTgg                                                                                 <  2:980913/150‑1 (MQ=255)
                                                                  cACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCAtt                                                                                   <  2:257157/147‑1 (MQ=255)
                                                                  cACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCAtt                                                                                   >  1:257157/1‑147 (MQ=255)
                                                                   aCGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCGGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCa                                                                               >  2:1070470/1‑150 (MQ=255)
                                                                     gtcgtAGGTCAGCGCTTCAGAGGCGTTGNNGCTCGGGTCGGTCTNTNCNAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATc                                                                             >  1:1690847/1‑150 (MQ=255)
                                                                        gtAGGCCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATcgcg                                                                          >  2:1252244/1‑150 (MQ=255)
                                                                           ggTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGcc                                                                                                        <  1:1145354/117‑1 (MQ=255)
                                                                           ggTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGcc                                                                                                        >  2:1145354/1‑117 (MQ=255)
                                                                           ggTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCAc                                                                       <  1:2176145/150‑1 (MQ=255)
                                                                             tCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGt                                                                                         >  2:1276223/1‑130 (MQ=255)
                                                                             tCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGt                                                                                         <  1:1276223/130‑1 (MQ=255)
                                                                               aGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCg                                                                   >  2:2203590/1‑150 (MQ=255)
                                                                                gcgcTTCAGGGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGa                                                                  >  2:2216718/1‑150 (MQ=255)
                                                                                gcgcTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGa                                                                  >  1:860218/1‑150 (MQ=255)
                                                                                  gcTTCAGATGCGTTGCCTGTCGGGGCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACg                                                                <  2:1941591/150‑1 (MQ=255)
                                                                                  gcTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACg                                                                >  1:373885/1‑150 (MQ=255)
                                                                                    ttCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgaggt                                                                                                                                          >  1:1200548/1‑74 (MQ=255)
                                                                                    ttCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgaggt                                                                                                                                          <  2:1200548/74‑1 (MQ=255)
                                                                                     tCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTa                                                             <  2:373885/150‑1 (MQ=255)
                                                                                                gCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATc                                                                             <  2:993689/123‑1 (MQ=255)
                                                                                                gCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATc                                                                             >  1:993689/1‑123 (MQ=255)
                                                                                                gCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTc                                                  >  1:1482566/1‑150 (MQ=255)
                                                                                                   gCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTcgcc                                               >  1:764098/1‑150 (MQ=255)
                                                                                                       gggTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGa                                           <  2:1482566/150‑1 (MQ=255)
                                                                                                        ggtcggtcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGAc                                          >  1:507212/1‑150 (MQ=255)
                                                                                                              tcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATAAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGACATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCa                                    >  2:566049/1‑150 (MQ=255)
                                                                                                              tcTGTACCAAGATCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCa                                    >  2:1297605/1‑150 (MQ=255)
                                                                                                                 gTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAgg                                 >  2:1838392/1‑150 (MQ=255)
                                                                                                                     cAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGAACCAACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGCCCGCCCAAGGTCAg                             >  1:1020308/1‑150 (MQ=255)
                                                                                                                      aGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCAACTAGACGGAAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAgg                            <  1:1862364/150‑1 (MQ=255)
                                                                                                                          tCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGt                        >  2:631628/1‑150 (MQ=255)
                                                                                                                           cAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTc                       >  1:1760219/1‑150 (MQ=255)
                                                                                                                           cAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTc                       >  2:461234/1‑150 (MQ=255)
                                                                                                                            aGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGTCGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCg                      >  2:1738394/1‑150 (MQ=255)
                                                                                                                                  tctctctgcagcgcagcCGTGTCTATGTCTATCAGGAAATTACCTTTCGCGGTCACTGTGCCGGTACCGACCAGCATGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcg                <  2:464004/143‑1 (MQ=38)
                                                                                                                                   gaCCCGGCTGCGCTGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcga               >  2:299600/2‑150 (MQ=255)
                                                                                                                                    aCCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgac              >  1:1909360/1‑150 (MQ=255)
                                                                                                                                    aCCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgac              >  1:1345128/1‑150 (MQ=255)
                                                                                                                                    aCCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGCCAg                             >  2:1020359/1‑135 (MQ=255)
                                                                                                                                    aCCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGCCAg                             <  1:1020359/135‑1 (MQ=255)
                                                                                                                                     cccGGCTGCGCTGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacg             >  1:933862/1‑150 (MQ=255)
                                                                                                                                     cccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacg             >  1:2407878/1‑150 (MQ=255)
                                                                                                                                      ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACTGACAAGCGTGCCATTGGGATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGACCAAGGTCCGGCCGTCGGCGCcgaaga            >  1:464004/1‑150 (MQ=255)
                                                                                                                                      ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAg                                  <  2:432090/128‑1 (MQ=255)
                                                                                                                                      ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAg                                  >  1:432090/1‑128 (MQ=255)
                                                                                                                                      ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacga            >  2:2073330/1‑150 (MQ=255)
                                                                                                                                      ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacga            >  2:1276279/1‑150 (MQ=255)
                                                                                                                                             cgcagcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAAtgct     <  1:1297605/150‑1 (MQ=255)
                                                                                                                                               cagcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCGGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAAtgctgc   <  2:5180/150‑1 (MQ=255)
                                                                                                                                                agcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATgctgct  <  2:507212/150‑1 (MQ=255)
                                                                                                                                                agcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATgctgct  <  2:764098/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATGCTGCT  >  NC_002947/937627‑937920

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