Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,912,492 T→A 7.0% D341V (GAC→GTC)  PP_1712 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,912,4920TA 7.0% 239.3 / 4.2 78D341V (GAC→GTC) PP_1712hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base A (3/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGATTAGCCCCCCCTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCGAGAGACCAC  >  NC_002947/1912353‑1912639
                                                                                                                                           |                                                                                                                                                   
cccGATTAGCCCCCCCTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGtcgatcga                                                                                                                                           <  2:192971/150‑1 (MQ=255)
 ccGATTAGCCCCCCCTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGag                                                                                                                                          <  2:712899/150‑1 (MQ=255)
 ccGATTAGCCCCCCCTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGag                                                                                                                                          >  1:2342213/1‑150 (MQ=255)
   gATTAGCCCCCCCTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGagta                                                                                                                                        <  1:1888689/150‑1 (MQ=255)
         cccccccTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTg                                                                                                                                  >  2:2140046/1‑150 (MQ=255)
           cccccTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGAtcgtcgt                                                                                                                                                  <  2:6941/132‑1 (MQ=255)
           cccccTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGAtcgtcgt                                                                                                                                                  >  1:6941/1‑132 (MQ=255)
           cccccTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGc                                                                                                                                >  1:1113381/1‑150 (MQ=255)
           cccccTCGGTTTCTTCAACATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGc                                                                                                                                >  1:1900275/1‑150 (MQ=255)
            ccccTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCAGTAGAGGGTCGGTTTGCTGCGATCGTCGTCGATCGAGTAGTGCTGGcc                                                                                                                               >  2:1308105/1‑150 (MQ=255)
             cccTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCg                                                                                                                              >  2:352316/1‑150 (MQ=255)
                      tcttcAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTg                                                                                                                     <  2:1296603/150‑1 (MQ=255)
                      tcttcAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTg                                                                                                                     <  2:1086908/150‑1 (MQ=255)
                        ttcAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGAt                                                                                                                   >  1:1278809/1‑150 (MQ=255)
                        ttcAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGAt                                                                                                                   >  2:1583937/1‑150 (MQ=255)
                          cAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGagta                                                                                                                                        <  1:1934309/127‑1 (MQ=255)
                          cAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGagta                                                                                                                                        >  2:1934309/1‑127 (MQ=255)
                          cAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGc                                                                                                                 <  1:1583937/150‑1 (MQ=255)
                          cAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGc                                                                                                                 <  2:2342213/150‑1 (MQ=255)
                             tatGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGATTCAATTCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGtt                                                                                                              <  1:747242/150‑1 (MQ=255)
                              atGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTg                                                                                                             >  2:2451655/1‑150 (MQ=255)
                                 ggACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGcc                                                                                                                               >  1:346278/1‑129 (MQ=255)
                                 ggACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGcc                                                                                                                               <  2:346278/129‑1 (MQ=255)
                                        ccgccgTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGtct                                                                                                   <  1:2140046/150‑1 (MQ=255)
                                        ccgccgTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGtct                                                                                                   >  1:2346570/1‑150 (MQ=255)
                                                gATGCCGAAACGCACTACCTGCCCCTGGTTGCCATGGTGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGt                                                                                           >  1:1774881/1‑150 (MQ=255)
                                                 aTGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTg                                                                                          >  2:1852754/1‑150 (MQ=255)
                                                           gCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGcc                                                                                                                               <  2:10018/103‑1 (MQ=255)
                                                           gCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGcc                                                                                                                               >  1:10018/1‑103 (MQ=255)
                                                           gCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTAt                                                                                <  1:496102/150‑1 (MQ=255)
                                                              cTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTAGTCTGTTTATTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATtga                                                                             <  2:1774881/150‑1 (MQ=255)
                                                               tACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTAAAGAAGCTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATtgat                                                                            <  2:665648/150‑1 (MQ=255)
                                                                  cTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTgatgat                                                                         >  1:1806525/1‑150 (MQ=255)
                                                                    gCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGGGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATcc                                                                       >  1:590505/1‑150 (MQ=255)
                                                                        cTGGTTGAGCTCGCGGGCAGGCTCGCAAGCATAGTAAAACCGTGAGAGTGTCGGGTGGGTGCGATCGTCGTTGATTGGGTGGGGGAAGCCCAACAACGGGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGAGGAGCCCTTc                                                                   >  1:2388497/1‑150 (MQ=255)
                                                                        cTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTc                                                                   <  2:2346570/150‑1 (MQ=255)
                                                                           gTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTc                                                                <  2:1421808/150‑1 (MQ=255)
                                                                           gTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTc                                                                >  2:548030/1‑150 (MQ=255)
                                                                           gTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTc                                                                <  2:1113381/150‑1 (MQ=255)
                                                                            ttGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGGCGTTTTGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCt                                                               >  1:1233418/1‑150 (MQ=255)
                                                                               aaCTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGAc                                                            <  1:458212/150‑1 (MQ=255)
                                                                                  tCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATtga                                                                             <  2:137862/130‑1 (MQ=255)
                                                                                  tCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATtga                                                                             >  1:137862/1‑130 (MQ=255)
                                                                                  tCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGAc                                                            <  2:910243/147‑1 (MQ=255)
                                                                                  tCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGAc                                                            >  1:910243/1‑147 (MQ=255)
                                                                                   cgcgGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAgt                                                        <  1:1852754/150‑1 (MQ=255)
                                                                                    gcgGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGTTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAgtg                                                       <  1:1308105/150‑1 (MQ=255)
                                                                                      gggCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAgtggt                                                     <  1:1885861/150‑1 (MQ=255)
                                                                                       ggCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAgtggtg                                                    <  1:2451655/150‑1 (MQ=255)
                                                                                           ggCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCAGCGATCGGCGTCGATCGAGTAGTGCTGGCTGTACAACGCCACGAGTTGGTGGCGGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCGGACCAGTGGTTGccc                                                >  1:1751252/1‑150 (MQ=255)
                                                                                            gCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGt                                                                                                               >  2:322658/1‑86 (MQ=255)
                                                                                            gCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGt                                                                                                               <  1:322658/86‑1 (MQ=255)
                                                                                            gCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCt                                                                                  <  1:787308/115‑1 (MQ=255)
                                                                                            gCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCt                                                                                  >  2:787308/1‑115 (MQ=255)
                                                                                              tCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCAcc                                             >  1:1705161/1‑150 (MQ=255)
                                                                                                gCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGa                                                             <  1:859693/132‑1 (MQ=255)
                                                                                                gCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGa                                                             >  2:859693/1‑132 (MQ=255)
                                                                                                 ccAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCg                                          >  1:1070293/1‑150 (MQ=255)
                                                                                                      aTAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGaa                                     <  2:1705161/150‑1 (MQ=255)
                                                                                                          tCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTg                                 >  1:1188869/1‑150 (MQ=255)
                                                                                                            cAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCACTATTGATGATCCCTTCGTCTGACCAGtggtgg                                                   <  1:1123442/130‑1 (MQ=255)
                                                                                                            cAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCACTATTGATGATCCCTTCGTCTGACCAGtggtgg                                                   >  2:1123442/1‑130 (MQ=255)
                                                                                                              gCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGccc                                                >  1:477597/1‑131 (MQ=255)
                                                                                                              gCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGccc                                                <  2:477597/131‑1 (MQ=255)
                                                                                                              gCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCAt                             <  2:1806525/150‑1 (MQ=255)
                                                                                                                   tAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCgaga                        >  1:639686/1‑150 (MQ=255)
                                                                                                                   tAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCgaga                        <  1:352316/150‑1 (MQ=255)
                                                                                                                              ttGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGGCCCGTCAg             >  2:1314489/1‑150 (MQ=255)
                                                                                                                              ttGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAg             >  1:1070563/1‑150 (MQ=255)
                                                                                                                                ctactcGAACGACGACGATCGAGTAGTGCTGGGCGTACAAAGTGATGCCGTGGCTGCTGTCTTCTACCGTGGAGCCCCTATTTATGATCCCTTCCTCTGACCCGTGGTGGCCCACCCCCCCGAATTTGTCACCGAGGCGCCCCCTCcgcg           >  2:294322/6‑150 (MQ=255)
                                                                                                                                 tactcGATCGACGACGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCga          >  2:10436/5‑150 (MQ=255)
                                                                                                                                 gCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGACGCCGCTTTTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCga          <  2:1900275/150‑1 (MQ=255)
                                                                                                                                  cTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGc                                         <  1:630660/118‑1 (MQ=255)
                                                                                                                                  cTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGc                                         >  2:630660/1‑118 (MQ=255)
                                                                                                                                   tGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCgaga        >  1:2240514/1‑150 (MQ=255)
                                                                                                                                   ctcGATCGACGACGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCgaga        <  1:294322/148‑1 (MQ=255)
                                                                                                                                   ctcGATCGACGACGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCgaga        >  1:693371/3‑150 (MQ=255)
                                                                                                                                   ctcGATCGACGACGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCgaga        <  1:1489076/148‑1 (MQ=255)
                                                                                                                                   ctcGATCGACGACGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCgaga        <  1:1489649/148‑1 (MQ=255)
                                                                                                                                     cGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCgagaga      <  2:2388497/150‑1 (MQ=255)
                                                                                                                                        tcgtcgtcgATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCGAGAGACCa   <  2:639686/150‑1 (MQ=255)
                                                                                                                                         cgtcgtcgATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCGAGAGACCAc  <  1:1314489/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CCCGATTAGCCCCCCCTCGGTTTCTTCAATATGGGACCACCCGCCGTCGATGCCGAAACGCACTACCTGCCCCTGGTTGAACTCGCGGGCAGGCTCGCCAGCATAGTCCAGCCGGTAGAGGGTCGGTTGGCTGCGATCGTCGTCGATCGAGTAGTGCTGGCCGTACAACGTGATGCGTTGGCTGCTGTCTTCTACCGTGAAGCCGCTATTGATGATCCCTTCGTCTGACCAGTGGTGGCCCACCTCGCCGAATTTGTCATCGAGACGCCCCGTCAGCGAGAGACCAC  >  NC_002947/1912353‑1912639

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