Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,651,399 C→T 5.2% P80S (CCG→TCG)  PP_3217 → aliphatic sulfonate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,651,3990CT 5.2% 223.2 / 2.2 77P80S (CCG→TCG) PP_3217aliphatic sulfonate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (44/29);  new base T (2/2);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

TACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATCT  >  NC_002947/3651255‑3651543
                                                                                                                                                |                                                                                                                                                  
tACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGAATGCTTGTTGACTCGCCGccc                                                                                                                                               <  2:1911676/150‑1 (MQ=255)
 aCAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCCGGCTCCCCGGAGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGcccc                                                                                                                                              <  2:1606630/150‑1 (MQ=255)
   aaCGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTg                                                                                                                                            >  2:256133/1‑150 (MQ=255)
      gCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGGCGCGTCGcccagcagg                                                                                                                                         <  2:46071/150‑7 (MQ=255)
        tCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCt                                                                                                                                       >  1:719519/1‑150 (MQ=255)
          gTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCg                                                                                                                                     <  2:719519/150‑1 (MQ=255)
           tGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCCCGCCGTGGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCgg                                                                                                                                    >  1:1669979/1‑150 (MQ=255)
                       aGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGgc                                                                                                                        <  2:1727970/150‑1 (MQ=255)
                       aGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCCATCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGgc                                                                                                                        >  1:1771428/1‑150 (MQ=255)
                          cgccgGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGcc                                                                                                                     >  2:1219017/1‑150 (MQ=255)
                          cgccgGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGcc                                                                                                                     >  1:1777705/1‑150 (MQ=255)
                                 gTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCgg                                                                                                                                    >  1:653637/1‑128 (MQ=255)
                                 gTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCgg                                                                                                                                    <  2:653637/128‑1 (MQ=255)
                                  tGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAgg                                                                                                             >  2:785777/1‑150 (MQ=255)
                                  tGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGACCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAgg                                                                                                             >  1:1636825/1‑150 (MQ=255)
                                   gatAGAACAACTGCCTTACCCCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGt                                                                                                            <  2:204598/148‑1 (MQ=255)
                                    cTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTg                                                                                                           >  2:306896/1‑150 (MQ=255)
                                     tAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGa                                                                                                          >  2:937017/1‑150 (MQ=255)
                                       gAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGAtt                                                                                                        >  2:1328420/1‑150 (MQ=255)
                                          caacTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                     <  1:256133/150‑1 (MQ=255)
                                              tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGcgccgc                                                                                                                                                 >  2:1518385/1‑102 (MQ=255)
                                              tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGcgccgc                                                                                                                                                 <  1:1518385/102‑1 (MQ=255)
                                              tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                     <  2:1711180/146‑1 (MQ=255)
                                              tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                     >  1:1711180/1‑146 (MQ=255)
                                              tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCa                                                                                                 >  2:1276048/1‑150 (MQ=255)
                                              tGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCa                                                                                                 <  2:1777705/150‑1 (MQ=255)
                                                ccTTACACCGTGGAGTGGAAGCAGTTCACCGCAGGCCCATGGACGGCCGAATCGCTCAACGTCGGCAGTCTCGGGATTTGCCTGCTGGGCGACGCGCCGCCCCTGTGCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGc                                                                                               <  1:424406/150‑1 (MQ=255)
                                                  ttACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCg                                                                                             >  2:1867457/1‑150 (MQ=255)
                                                  ttACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCg                                                                                             >  1:1099021/1‑150 (MQ=255)
                                                  ttACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCg                                                                                             <  1:537300/150‑1 (MQ=255)
                                                      acCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAg                                                                                         >  1:573914/1‑150 (MQ=255)
                                                        ccTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTg                                                                                                       <  1:684615/134‑1 (MQ=255)
                                                        ccTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTg                                                                                                       >  2:684615/1‑134 (MQ=255)
                                                        ccTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGaa                                                                                       <  2:573914/150‑1 (MQ=255)
                                                        ccTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGaa                                                                                       <  1:785777/150‑1 (MQ=255)
                                                        ccTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGaa                                                                                       >  2:928631/1‑150 (MQ=255)
                                                         cTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAAc                                                                                      >  2:1869714/1‑150 (MQ=255)
                                                                ggAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAg                                                                               <  1:1219017/150‑1 (MQ=255)
                                                                 gAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAg‑‑g                                                                            >  2:546494/1‑150 (MQ=255)
                                                                    gCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGt                                                                         >  1:1234944/1‑150 (MQ=255)
                                                                      aGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTgg                                                                       <  1:367385/150‑1 (MQ=255)
                                                                        ttCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGcc                                                                     >  2:1306006/1‑150 (MQ=255)
                                                                             cGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGGTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCt                                                                >  2:833362/1‑150 (MQ=255)
                                                                              gCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTg                                                                                                                             >  1:1384838/1‑90 (MQ=255)
                                                                              gCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTg                                                                                                                             <  2:1384838/90‑1 (MQ=255)
                                                                              gCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTg                                                               <  2:855911/150‑1 (MQ=255)
                                                                               cAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTgg                                                              <  2:1568170/150‑1 (MQ=255)
                                                                                      cccGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCaa                                                       >  2:1497863/1‑150 (MQ=255)
                                                                                       ccGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAg                                                      <  1:1328420/150‑1 (MQ=255)
                                                                                       ccGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAg                                                      >  2:1673318/1‑150 (MQ=255)
                                                                                            ggCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATtc                                                 <  1:1276048/150‑1 (MQ=255)
                                                                                             gCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGccc                                                                                                                                               <  1:745713/57‑1 (MQ=255)
                                                                                             gCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGccc                                                                                                                                               >  2:745713/1‑57 (MQ=255)
                                                                                             gCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGGCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATtct                                                >  1:972452/1‑150 (MQ=255)
                                                                                               cGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTcc                                              >  2:1908211/1‑150 (MQ=255)
                                                                                                 aaGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACTTCCTAG‑‑TGTTGGCCATCCTGGCGCGCAAGGATTCTCCGa                                            <  1:833362/150‑1 (MQ=255)
                                                                                                  aGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGTCGGCAGAACCTCGAAG‑‑GCGTg                                                                        <  2:606117/121‑1 (MQ=255)
                                                                                                  aGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGTCGGCAGAACCTCGAAG‑‑GCGTg                                                                        >  1:606117/1‑121 (MQ=255)
                                                                                                       tCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCcaca                                      >  2:1501306/1‑150 (MQ=255)
                                                                                                       tCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCcaca                                      >  2:1215650/1‑150 (MQ=255)
                                                                                                       tCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCcaca                                      >  2:1097868/1‑150 (MQ=255)
                                                                                                              cGGCAGTCTCGATATTGGCCTGCTGGGCAACGAGCACCAGGTGGGCGGCGTTGCCGGGGGGGCGCCGATCAAATGCATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCa                               <  1:1867457/150‑1 (MQ=255)
                                                                                                               ggCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAg                              <  1:546494/150‑1 (MQ=255)
                                                                                                                    tctcGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCt                         >  2:32280/1‑150 (MQ=255)
                                                                                                                      tcGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTgc                       <  1:306896/150‑1 (MQ=255)
                                                                                                                      tcGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTgc                       >  1:918120/1‑150 (MQ=255)
                                                                                                                              ggCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCaaa               >  2:248980/1‑150 (MQ=255)
                                                                                                                               gCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                     >  1:1171204/1‑65 (MQ=255)
                                                                                                                               gCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCg                                                                                                     <  2:1171204/65‑1 (MQ=255)
                                                                                                                               gCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGCCGTCAAGTTGCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGCCCTGCGCGGCAAAc              >  1:594016/1‑150 (MQ=255)
                                                                                                                                 ctgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGaa                                                                                       >  2:1626163/1‑77 (MQ=255)
                                                                                                                                 ctgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGaa                                                                                       <  1:1626163/77‑1 (MQ=255)
                                                                                                                                  tgctgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAAcgcg           >  2:1330003/1‑150 (MQ=255)
                                                                                                                                     tgGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCgg        >  2:1765864/1‑150 (MQ=255)
                                                                                                                                     ggggCGGCGCGTCGCCCCTGTTCCTCGGTGCCTTGTGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTgg                                                                       >  2:934242/2‑87 (MQ=255)
                                                                                                                                     ggggCGGCGCGTCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTgg                                                                       <  1:934242/86‑1 (MQ=255)
                                                                                                                                     ggggCGGCGCGTCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCgg        >  1:276533/2‑150 (MQ=255)
                                                                                                                                          gACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATc   >  2:1236164/1‑150 (MQ=255)
                                                                                                                                           aCGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATCt  >  1:1323526/1‑150 (MQ=255)
                                                                                                                                             gcgcCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAg                              <  2:336398/120‑1 (MQ=255)
                                                                                                                                             gcgcCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAg                              >  1:336398/1‑120 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TACAACGCTCGTGCCGTGATGGAAGCCGCCGGCGTGCTAGAACAACTGCCTTACACCCTGGAGTGGAAGCAGTTCACCGCAGGCTCCCCGGTGGCCGAAGCGCTCAACGTCGGCAGTCTCGATATTGGCCTGCTGGGCGACGCGCCGCCCCTGTTCCTCGGTGCCTTGGGGGCGCCGATCAAGGTGATTGCGGTCAGCCGGCAGAACCTCGAAG‑‑GCGTGGCCATCCTGGCGCGCAAGGATTCTCCGATCCACAGCCTGCAGGACCTGCGCGGCAAACGCGCGGCCATCT  >  NC_002947/3651255‑3651543

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