Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,939,164 G→T 8.5% intergenic (‑66/‑113) PP_1738 ← / → PP_1739 hypothetical protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,939,1640GT 8.5% 223.2 / 6.8 70intergenic (‑66/‑113)PP_1738/PP_1739hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base T (3/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCTGCCG  >  NC_002947/1939016‑1939306
                                                                                                                                                    |                                                                                                                                              
gCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGt                                                                                                                                               >  1:450375/1‑150 (MQ=255)
 cAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCGAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGtt                                                                                                                                              <  1:1216453/150‑1 (MQ=255)
             tcgacgcTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGt                                                                                                                                  <  1:715570/146‑1 (MQ=255)
              ccGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTg                                                                                                                                 >  1:1385915/1‑150 (MQ=255)
                gcgcTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGc                                                                                                                               <  2:376434/150‑1 (MQ=255)
                      aaGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTgg                                                                                                                         <  2:527967/150‑1 (MQ=255)
                      aaGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTgg                                                                                                                         >  1:972591/1‑150 (MQ=255)
                             cacaGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACtt                                                                                                                  <  2:450375/150‑1 (MQ=255)
                               caGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGc                                                                                                                <  2:1385915/150‑1 (MQ=255)
                                    aTGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCt                                                                                                           <  2:972591/150‑1 (MQ=255)
                                       aGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCa                                                                                                        >  2:1493531/1‑150 (MQ=255)
                                              tGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAACTTTGCCCGcctggt                                                                                                                                   >  2:1608924/1‑111 (MQ=255)
                                              tGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAACTTTGCCCGcctggt                                                                                                                                   <  1:1608924/116‑6 (MQ=255)
                                              tGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGc                                                                                                 >  1:830820/1‑150 (MQ=255)
                                                    gCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCgg                                                                                           <  2:208816/150‑1 (MQ=255)
                                                              gAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATGTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGAttt                                                                                 >  2:1418658/1‑150 (MQ=255)
                                                                 gcgcAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGc                                                                              >  2:1437190/1‑150 (MQ=255)
                                                                     aaaGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTa                                                                          >  1:426491/1‑150 (MQ=255)
                                                                         aTCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGaaa                                                                      >  2:339558/1‑150 (MQ=255)
                                                                         aTCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGaaa                                                                      >  1:708331/1‑150 (MQ=255)
                                                                          tCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAAc                                                                     >  2:701776/1‑150 (MQ=255)
                                                                            gAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   <  2:139255/150‑1 (MQ=255)
                                                                                cATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCAt                                                               <  2:1554357/150‑1 (MQ=255)
                                                                                cATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGCGGGCAACTTTGCCCGcctg                                                                                                                                     <  1:1250995/79‑4 (MQ=37)
                                                                                cATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGCGGGCAACTTTGCCCGcctg                                                                                                                                     >  2:1250995/1‑76 (MQ=37)
                                                                                    tAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTACGATATGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGc                                                           >  1:515605/1‑150 (MQ=255)
                                                                                     aaTCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTg                                                                                                                          <  1:245681/86‑1 (MQ=255)
                                                                                     aaTCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTg                                                                                                                          >  2:245681/1‑86 (MQ=255)
                                                                                       tCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGAt                                                        <  2:426491/150‑1 (MQ=255)
                                                                                         ccccAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTAc                                                                                                                    >  1:1203252/1‑88 (MQ=255)
                                                                                         ccccAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTAc                                                                                                                    <  2:1203252/88‑1 (MQ=255)
                                                                                          cccAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCAt                                                     >  2:1558260/1‑150 (MQ=255)
                                                                                             aaaaGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCAt                                                                                                                             >  1:647796/1‑75 (MQ=255)
                                                                                             aaaaGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCAt                                                                                                                             <  2:647796/75‑1 (MQ=255)
                                                                                              aaaGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGAt                                                 <  1:1493531/150‑1 (MQ=255)
                                                                                                   aCGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCttt                                            <  1:701776/150‑1 (MQ=255)
                                                                                                     gCAAGGTCGAGCAGTGTGCAGGCTATGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCg                                          >  1:703960/1‑150 (MQ=255)
                                                                                                      cagggTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGa                                         >  1:1643078/4‑150 (MQ=255)
                                                                                                       aaGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAATTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGAt                                        >  2:1028279/1‑150 (MQ=255)
                                                                                                       aaGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGAt                                        <  2:1643078/150‑1 (MQ=255)
                                                                                                            cGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAATTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCg                                   <  1:1028279/150‑1 (MQ=255)
                                                                                                                cAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTAt                               <  1:1437190/150‑1 (MQ=255)
                                                                                                                 aGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTAtt                              <  1:1418658/150‑1 (MQ=255)
                                                                                                                  gtgtGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTg                             >  1:202184/1‑150 (MQ=255)
                                                                                                                       cAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGaa                        >  2:527170/1‑150 (MQ=255)
                                                                                                                       cAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGaa                        <  1:1558260/150‑1 (MQ=255)
                                                                                                                        aGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   <  1:1467416/106‑1 (MQ=255)
                                                                                                                        aGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   >  2:1467416/1‑106 (MQ=255)
                                                                                                                        aGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAAc                       >  1:536704/1‑150 (MQ=255)
                                                                                                                          gCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   >  2:1487268/1‑104 (MQ=255)
                                                                                                                          gCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   <  1:1487268/104‑1 (MQ=255)
                                                                                                                          gCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAg                     >  1:152958/1‑150 (MQ=255)
                                                                                                                           cTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    >  1:651479/1‑150 (MQ=255)
                                                                                                                            tAGGGCCACAACCAGGCGGGGAAAGTTGCCTGCTGGGTTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGa                   <  2:158018/150‑1 (MQ=255)
                                                                                                                             aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATc                                                                                                            <  2:1419027/60‑1 (MQ=255)
                                                                                                                             aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATc                                                                                                            >  1:1419027/1‑60 (MQ=255)
                                                                                                                             aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGgtgt                                                                                      <  1:174628/82‑1 (MQ=255)
                                                                                                                             aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGgtgt                                                                                      >  2:174628/1‑82 (MQ=255)
                                                                                                                             aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGAt                  <  2:515605/150‑1 (MQ=255)
                                                                                                                               gACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGt                >  2:1144185/1‑150 (MQ=255)
                                                                                                                                  cacaACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACGTAGCTGCCGGGGTGTGAtttt                                                                                <  1:1147538/83‑1 (MQ=255)
                                                                                                                                  cacaACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGAtttt                                                                                >  2:1147538/1‑83 (MQ=255)
                                                                                                                                  cacaACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGtctc             >  2:787524/1‑150 (MQ=255)
                                                                                                                                      aCCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTc         >  2:628331/1‑150 (MQ=255)
                                                                                                                                         aGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCt      <  1:527170/150‑1 (MQ=255)
                                                                                                                                             gggCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAAAAGGGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCTGCCg  >  1:40771/1‑150 (MQ=255)
                                                                                                                                              ggCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    <  1:599220/131‑1 (MQ=255)
                                                                                                                                              ggCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    >  2:599220/1‑131 (MQ=255)
                                                                                                                                                 aaaGTTGCCCGCTGCGGTGGCATTTGACGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    <  1:1218742/128‑1 (MQ=255)
                                                                                                                                                 aaaGTTGCCCGCTGCGGTGGCATTTGACGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    >  2:1218742/1‑128 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCTGCCG  >  NC_002947/1939016‑1939306

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