Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 205,404 T→A 5.4% L3637Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947205,4040TA 5.4% 207.4 / 4.4 75L3637Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base A (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACT  >  NC_002947/205255‑205549
                                                                                                                                                     |                                                                                                                                                 
cACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGCTCACCATCACCGCCGACCCGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCt                                                                                                                                                   <  1:628084/150‑1 (MQ=17)
  ccGACGAGCCTGGTACCCCCGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTgg                                                                                                                                                 <  2:1275332/150‑1 (MQ=17)
   cgacgaGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGt                                                                                                                                                >  2:1581122/1‑150 (MQ=32)
        aGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                           <  2:614868/150‑1 (MQ=32)
        aGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                           <  1:866906/150‑1 (MQ=32)
         gCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAc                                                                                                                                          <  1:222664/150‑1 (MQ=32)
           cTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCCCCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAccc                                                                                                                                        >  2:957310/1‑150 (MQ=17)
            tGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGAgcca                                                                                                                                       <  2:350901/150‑2 (MQ=12)
            tGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCCGTAGGGGCAAGGGCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                       <  1:511700/150‑1 (MQ=11)
              gTACCCCAGGCAACCCTGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                     <  2:1642557/150‑1 (MQ=17)
              gTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                     <  2:374301/150‑1 (MQ=32)
              gTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCAAGACCTGGTAGTGACCCTGa                                                                                                                                     <  2:762409/150‑1 (MQ=17)
                 ccccAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGca                                                                                                                                  <  2:1356394/150‑1 (MQ=32)
                   ccAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                >  1:1126333/1‑150 (MQ=32)
                      ggCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                             <  2:964135/150‑1 (MQ=33)
                      ggCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                             <  2:811406/150‑1 (MQ=33)
                       gCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACg                                                                                                                            >  1:302410/1‑150 (MQ=34)
                         aaCCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                          <  2:1872832/150‑1 (MQ=35)
                           cccGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGCCCTGGTCGTGACCCTGAGCAACAACGCCCa                                                                                                                        <  2:453910/150‑1 (MQ=25)
                            ccGGGCGTCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                       <  1:734804/150‑1 (MQ=25)
                             cGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAgg                                                                                                                      >  1:8127/1‑150 (MQ=35)
                                     agcaACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATc                                                                                                              >  1:223339/1‑150 (MQ=33)
                                                    cTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAAcc                                                                                               >  2:303032/1‑150 (MQ=38)
                                                       gtctaGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                            <  1:1226357/146‑1 (MQ=255)
                                                          aaGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                         >  2:1891162/1‑150 (MQ=38)
                                                            ggTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAATGCCGGCGAAACCAGCGCGcc                                                                                       <  1:481317/150‑1 (MQ=38)
                                                             gTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCg                                                                                      >  1:1233854/1‑150 (MQ=38)
                                                                      aCCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgc                                                                                          <  2:1345867/137‑1 (MQ=39)
                                                                      aCCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgc                                                                                          >  1:1345867/1‑137 (MQ=39)
                                                                      aCCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCAc                                                                             <  1:957310/150‑1 (MQ=38)
                                                                      aCCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCATGAGCAACAACGCCCAGGTCACCATCAAGGCCTGCGAAACCAGCGCGCCGTACACCCAc                                                                             <  2:334820/150‑1 (MQ=255)
                                                                             aCCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTgcgc                                                                      >  1:1588384/1‑150 (MQ=35)
                                                                              ccAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGTCGAAACCAGCGCGCCGTACACCCACGCTGCGCa                                                                     <  1:1423357/150‑1 (MQ=25)
                                                                               cAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACATCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAg                                                                    >  1:1575630/1‑150 (MQ=25)
                                                                                    ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGa                                                               >  2:851401/1‑150 (MQ=38)
                                                                                    ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCACGCCGTACGCCCACGCTGCGCAGGGCGa                                                               >  2:1276693/1‑150 (MQ=255)
                                                                                     tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                              <  1:851401/150‑1 (MQ=38)
                                                                                     tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                              <  1:1891162/150‑1 (MQ=38)
                                                                                       ggTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGa                                                            >  1:605332/1‑150 (MQ=38)
                                                                                        gTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACTACCGGTTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGAc                                                           <  2:1232409/150‑1 (MQ=255)
                                                                                             cGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTTAGCAACAACGCCCAGGTCACCATCAAGGGCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTa                                                      >  2:955563/1‑150 (MQ=255)
                                                                                              gagaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTAc                                                     >  1:196824/1‑150 (MQ=39)
                                                                                               agaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACa                                                    >  2:481812/1‑150 (MQ=39)
                                                                                                    gTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGAt                                               >  2:1346182/1‑150 (MQ=39)
                                                                                                     tCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACCACCTCCAGTTTACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATg                                              <  2:1126333/150‑1 (MQ=11)
                                                                                                     tCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATg                                              >  2:345168/1‑150 (MQ=38)
                                                                                                      cAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                              >  2:907594/1‑133 (MQ=34)
                                                                                                      cAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGAt                                                              <  1:907594/133‑1 (MQ=34)
                                                                                                        aaCCGACCTTCACTGTTCACGTCAACCAGCCACTGGCTCACGACCTGGTGGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTg                                           >  2:1191785/1‑150 (MQ=255)
                                                                                                           cGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGcc                                        >  1:1658246/1‑150 (MQ=35)
                                                                                                            gACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCt                                                       <  1:1221215/134‑1 (MQ=37)
                                                                                                            gACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCt                                                       >  2:1221215/1‑134 (MQ=37)
                                                                                                             aCCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAATGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAg                                      >  2:742828/1‑150 (MQ=25)
                                                                                                                    cTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGcc                               >  1:515103/1‑150 (MQ=33)
                                                                                                                     tGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCACGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCt                              <  1:1276693/150‑1 (MQ=17)
                                                                                                                       ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                            <  2:196824/150‑1 (MQ=31)
                                                                                                                       ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                            <  2:322953/150‑1 (MQ=31)
                                                                                                                         cACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGt                          <  2:8127/150‑1 (MQ=11)
                                                                                                                          aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTa                         >  1:1453377/1‑150 (MQ=31)
                                                                                                                          aCGTCAACCAGCCGCTGGCTCAAGACATGGTCGTGACCCTGAGCAACAACGCCCAGGACACCATCAAGGCCGGCGAAACCAGCGAGCAGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTa                         >  2:402471/1‑150 (MQ=11)
                                                                                                                           cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTAt                        >  2:1715844/1‑150 (MQ=32)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATc                       <  1:1715844/150‑1 (MQ=32)
                                                                                                                             tCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAATGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCa                      <  1:742828/150‑1 (MQ=17)
                                                                                                                               aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAAc                    >  1:995781/1‑150 (MQ=32)
                                                                                                                                     ccGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTcggcg              <  1:303032/150‑1 (MQ=32)
                                                                                                                                       tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCa                      >  2:1114627/7‑140 (MQ=11)
                                                                                                                                       tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCa                      <  1:1114627/134‑1 (MQ=255)
                                                                                                                                       gCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt            <  1:481812/150‑1 (MQ=32)
                                                                                                                                          gggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGa         >  1:293312/4‑150 (MQ=255)
                                                                                                                                           gCTCACGACCTGGTCGTGACCCTGAGCAACACCGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGAt        <  1:412316/150‑1 (MQ=17)
                                                                                                                                           gCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGAt        <  2:1575630/150‑1 (MQ=33)
                                                                                                                                           gCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGAt        >  2:1240540/1‑150 (MQ=33)
                                                                                                                                             tCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGc      >  2:1095291/1‑150 (MQ=34)
                                                                                                                                               aCGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTGGACAGCCATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCa    <  2:1524178/150‑1 (MQ=255)
                                                                                                                                               aCGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTACCAACTCGGCGGTGGATGCCa    <  1:1346182/150‑1 (MQ=25)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACt  >  2:1200080/1‑150 (MQ=35)
                                                                                                                                                     |                                                                                                                                                 
CACCGACGAGCCTGGTACCCCAGGCAACCCGGGCGGCAGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGATGCCACT  >  NC_002947/205255‑205549

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: