Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 200,748 T→A 5.4% L2085Q (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_002947200,7480TA 5.4% 285.9 / 6.2 93L2085Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base A (3/2);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

GTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGT  >  NC_002947/200600‑200896
                                                                                                                                                    |                                                                                                                                                    
gTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTg                                                                                                                                                     <  2:41182/150‑1 (MQ=255)
   aGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                                                                  <  1:427466/150‑1 (MQ=255)
       ccGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGCTATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGTCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              <  2:71452/150‑1 (MQ=255)
       ccGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              <  2:786540/150‑1 (MQ=255)
       ccGAGGTCAAAGACGAGCCAGGCACACCAGGCAACGAAGGCGATATCGTGAAGGTCACCAGCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              >  2:2202698/1‑150 (MQ=11)
         gAGGTGAACGCCGAGCCAGGCCCACCAAGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGAACAAGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCTCCGTGCACGACAACCAGCCGCTGGCCCACGACCTGGTCGTGAcc                                                                                                                                            <  2:409709/150‑1 (MQ=14)
             tGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAGGACCTGGTCGTGACCCTGa                                                                                                                                        <  1:1883506/150‑1 (MQ=255)
              gACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                       <  1:920609/150‑1 (MQ=255)
                        ccAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCt                                                                                                                             >  1:1405709/1‑150 (MQ=255)
                         ccggCACACCAGGTAACGAAGGCCATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCACCCCGACCTTCACCGTGCACGTAAACCATCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTc                                                                                                                            <  1:2224282/148‑1 (MQ=14)
                         cAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGCGAAAAACTTACGGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCCGAGCAACAACGCTc                                                                                                                            >  1:900340/1‑150 (MQ=14)
                          aGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                           <  2:1383606/150‑1 (MQ=255)
                           ggCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                          <  2:1334815/150‑1 (MQ=255)
                           ggCACACCAGGTAACGAAGGCGAGATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCCGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                          <  2:355257/150‑1 (MQ=255)
                              acacCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACAACGTCTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTc                                                                                                                       >  1:761144/1‑150 (MQ=17)
                                acCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                     <  2:1986462/150‑1 (MQ=255)
                                acCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                    >  2:2105190/1‑119 (MQ=255)
                                acCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                    <  1:2105190/119‑1 (MQ=255)
                                 ccAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                                                                  <  1:1309910/120‑1 (MQ=255)
                                 ccAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                                                                  >  2:1309910/1‑120 (MQ=255)
                                   aGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                                  >  2:123346/1‑150 (MQ=255)
                                   aGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCAAGTGATAGTGGCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                                  >  1:1233207/1‑150 (MQ=14)
                                    ggTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATc                                                                                                                 <  1:935290/150‑1 (MQ=255)
                                       aaCGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                              <  2:2034371/150‑1 (MQ=255)
                                          gAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAAAGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTTACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCt                                                                                                           >  1:261697/1‑150 (MQ=255)
                                           attgggATAACTTCAATGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTg                                                                                                          <  2:1714647/145‑1 (MQ=255)
                                           aaGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTg                                                                                                          <  2:239504/150‑1 (MQ=255)
                                             ggCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCATGAGCAACAACGCTCAGGTCACCATCAAGGCTGGc                                                                                                        <  2:422435/150‑1 (MQ=255)
                                                gATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                    >  2:2336261/1‑103 (MQ=255)
                                                gATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTgg                                                                                                                                                    <  1:2336261/103‑1 (MQ=255)
                                                 atatCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGaaa                                                                                                    <  1:251245/150‑1 (MQ=40)
                                                    tCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCa                                                                                                 >  1:382611/1‑150 (MQ=39)
                                                               aCCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTac                                                                                      >  2:152042/1‑150 (MQ=38)
                                                                 cATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGGCCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTa                                                                                       <  2:1051551/147‑1 (MQ=38)
                                                                 cATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTa                                                                                       >  1:1051551/1‑147 (MQ=38)
                                                                  aTCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAcc                                                                                   >  2:337127/1‑150 (MQ=38)
                                                                   tCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACAccc                                                                                  >  1:1713710/1‑150 (MQ=39)
                                                                      ccgccgACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGCGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCacg                                                                               >  2:1172925/1‑150 (MQ=37)
                                                                      ccgccgACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                               <  1:1413841/134‑1 (MQ=37)
                                                                      ccgccgACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAgc                                                                                               >  2:1413841/1‑134 (MQ=38)
                                                                        gccgACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCCGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCacgac                                                                             <  2:900340/150‑1 (MQ=37)
                                                                        gccgACCAGGCTTCGGTGGCCGAGAAAGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCacgac                                                                             <  2:261697/150‑1 (MQ=25)
                                                                          cgACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAcgc                                                                           <  2:2273331/150‑1 (MQ=38)
                                                                           gACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACTACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGAcgcg                                                                          <  1:152042/150‑1 (MQ=38)
                                                                              cAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCaa                                                                       <  2:761144/150‑1 (MQ=38)
                                                                              cAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCaa                                                                       <  1:2202698/150‑1 (MQ=39)
                                                                                 gCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCaa                                                                                                               >  2:2157274/1‑107 (MQ=35)
                                                                                 gCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCaa                                                                                                               <  1:2157274/107‑1 (MQ=35)
                                                                                   ttCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                >  2:2122615/1‑120 (MQ=33)
                                                                                   ttCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAg                                                                                                <  1:2122615/120‑1 (MQ=33)
                                                                                   ttCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGa                                                                  >  1:166090/1‑150 (MQ=35)
                                                                                     cGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGGACACCCACGACGCGCAAGGCGATg                                                                <  1:337127/150‑1 (MQ=25)
                                                                                         ggCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                            <  1:123346/150‑1 (MQ=35)
                                                                                         ggCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                            <  1:2185483/150‑1 (MQ=35)
                                                                                                    tCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATg                                                 >  2:2127226/1‑150 (MQ=38)
                                                                                                     cAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGc                                                <  1:1299520/150‑1 (MQ=38)
                                                                                                      aaaCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCg                                               >  2:2216388/1‑150 (MQ=39)
                                                                                                      aaaCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCCATCAGGATGCg                                               >  2:1317323/1‑150 (MQ=255)
                                                                                                        aCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCggg                                             >  1:2389392/1‑150 (MQ=40)
                                                                                                              ccagaCCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGAt                                       >  1:1892630/5‑150 (MQ=255)
                                                                                                                   ccGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGcc                                  >  2:227263/1‑150 (MQ=255)
                                                                                                                   ccGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGcc                                  >  2:2248978/1‑150 (MQ=40)
                                                                                                                       gCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCgg                              <  2:1405709/150‑1 (MQ=255)
                                                                                                                           gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATc                          >  1:1817878/1‑150 (MQ=255)
                                                                                                                            tCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCa                                     >  2:1864353/1‑138 (MQ=255)
                                                                                                                            tCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCa                                     <  1:1864353/138‑1 (MQ=255)
                                                                                                                              aaCCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACg                       >  2:65959/1‑150 (MQ=255)
                                                                                                                                 cAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTcc                    >  1:1076228/1‑150 (MQ=255)
                                                                                                                                  aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCaa                                                                       <  1:1469161/98‑1 (MQ=39)
                                                                                                                                  aGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCaa                                                                       >  2:1469161/1‑98 (MQ=39)
                                                                                                                                   gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCCTGAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGc                  >  1:1500252/1‑150 (MQ=255)
                                                                                                                                   gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                                                                      >  2:2139199/1‑98 (MQ=39)
                                                                                                                                   gCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                                                                      <  1:2139199/98‑1 (MQ=39)
                                                                                                                                       cTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCtgtg              <  1:890138/150‑1 (MQ=255)
                                                                                                                                        agggtcACGACCAGGTCGTGGCCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCa                            >  1:2174292/6‑135 (MQ=255)
                                                                                                                                        agggtcACGACCAGGTCGTGGCCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCa                            <  2:2174292/130‑1 (MQ=255)
                                                                                                                                         gggtcACGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGa            >  2:2060825/5‑150 (MQ=255)
                                                                                                                                         ggCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGa            >  1:168898/1‑150 (MQ=255)
                                                                                                                                            ccACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                            <  1:556840/99‑1 (MQ=38)
                                                                                                                                            ccACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                            >  2:556840/1‑99 (MQ=38)
                                                                                                                                            ccACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTc                                >  1:702656/1‑127 (MQ=255)
                                                                                                                                            ccACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTc                                <  2:702656/127‑1 (MQ=255)
                                                                                                                                            ccACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGt         >  1:1629760/1‑150 (MQ=255)
                                                                                                                                            ccACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGt         >  1:1331229/1‑150 (MQ=255)
                                                                                                                                            ccACAACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGt         >  2:479138/1‑150 (MQ=255)
                                                                                                                                             cACGACCTGGTCGTGACCCTGAGCAACAACGCTTATTTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTa        >  1:1770509/1‑150 (MQ=255)
                                                                                                                                               cGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATc                                      >  1:1798041/1‑118 (MQ=255)
                                                                                                                                               cGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATc                                      <  2:1798041/118‑1 (MQ=255)
                                                                                                                                                gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTGACCAGCAAGTCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGAc     <  2:1770509/150‑1 (MQ=255)
                                                                                                                                                  ccTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCg                                               >  2:2045212/1‑106 (MQ=255)
                                                                                                                                                  ccTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCg                                               <  1:2045212/106‑1 (MQ=255)
                                                                                                                                                  ccTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACgg   >  2:363003/1‑150 (MQ=255)
                                                                                                                                                   cTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGt  <  2:166090/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCAGCACCGAGGTGACCGACGAGCCAGGCACACCAGGTAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGGCTTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTCGGCATCACGTCCGCTGTGGATGTAGACGGT  >  NC_002947/200600‑200896

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: