Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,183 T→A 7.8% intergenic (+100/+56) yhbY → / ← greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,328,1830TA 7.8% 319.4 / 12.8 115intergenic (+100/+56)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base T (51/55);  new base A (3/6);  total (54/61)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 116 total aligned reads displayed.

CTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCA  >  NC_000913/3328039‑3328332
                                                                                                                                                |                                                                                                                                                      
cTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCgg                                                                                                                                                   <  2:1475302/150‑1 (MQ=255)
 tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGc                                                                                                                                                  <  2:896583/150‑1 (MQ=255)
  ttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGcc                                                                                                                                                 <  1:1449079/150‑1 (MQ=255)
     tCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCttt                                                                                                                                             <  2:870253/150‑1 (MQ=255)
       gCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCttttt                                                                                                                                           >  2:1029795/1‑150 (MQ=255)
       gCCCAACTAAAGAACGTAAAATCTCGCTGCCACACTAAGATTATCCAAAACTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCttttt                                                                                                                                           <  2:296394/150‑1 (MQ=255)
         ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCttttatc                                                                                                                                         <  1:395141/150‑4 (MQ=255)
         ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCttttatc                                                                                                                                         <  2:1508733/150‑4 (MQ=255)
         ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTcc                                                                                                                                         <  1:734262/150‑1 (MQ=255)
         ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTcc                                                                                                                                         >  1:1647865/1‑150 (MQ=255)
           aagtAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGGGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctt                                                                                                                                       <  2:866575/147‑1 (MQ=255)
             cTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCcttta                                                                                                                                     >  2:573505/1‑149 (MQ=255)
                aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGTCCGCATAGCGGCCTTTTTCctttcttt                                                                                                                                  <  1:262274/150‑1 (MQ=255)
                aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttcttt                                                                                                                                  <  1:121343/150‑1 (MQ=255)
                  taaCGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTAc                                                                                                                                <  2:879249/149‑1 (MQ=255)
                     cGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAAt                                                                                                                             >  1:128397/1‑150 (MQ=255)
                      gTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctt                                                                                                                                       >  1:88292/1‑139 (MQ=255)
                      gTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctt                                                                                                                                       <  2:88292/139‑1 (MQ=255)
                      gTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                            >  2:1187333/1‑150 (MQ=255)
                           aTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCCTAGCGGCCTTTTTCCTTTCTTTACAATACATCa                                                                                                                       <  1:714430/150‑1 (MQ=255)
                               cGCTGCCACGCTAATATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACAt                                                                                                                   <  1:1111194/150‑1 (MQ=255)
                               cGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACAt                                                                                                                   <  1:1105410/150‑1 (MQ=255)
                                     cACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGGGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                             <  2:1685540/150‑1 (MQ=255)
                                     cACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                             <  2:1103700/150‑1 (MQ=255)
                                      aCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGt                                                                                                            >  2:1561587/1‑150 (MQ=255)
                                       cGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                           >  2:1386334/1‑150 (MQ=255)
                                       cGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                           <  1:194519/150‑1 (MQ=255)
                                                  tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTc                                                                                                <  1:326134/150‑1 (MQ=255)
                                                   ccTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGACACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGACATTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                               >  1:270323/1‑150 (MQ=255)
                                                   ccTAAAGTTACACACATTTGCTGTGGAAAACGAGGGGTGTGCCGCAGGCAGGAGAGCAAATTGCCACGCTCTGTTCGTTGATAAAAGTCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                               <  1:262746/150‑1 (MQ=255)
                                                         gTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAgg                                                                                         <  1:1029795/150‑1 (MQ=255)
                                                         gTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAgg                                                                                         >  1:13077/1‑150 (MQ=255)
                                                            acacacATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTACTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTAACAGGTAt                                                                                      >  2:690188/1‑150 (MQ=255)
                                                             cacacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                     >  1:798008/1‑150 (MQ=255)
                                                              acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTc                                                                                    <  1:453801/150‑1 (MQ=255)
                                                               cacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTcc                                                                                   <  2:1237485/150‑1 (MQ=255)
                                                                 caTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAc                                                                                 <  1:190559/150‑1 (MQ=255)
                                                                 caTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAc                                                                                 >  1:618764/1‑150 (MQ=255)
                                                                 caTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTACTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAc                                                                                 >  2:116062/1‑150 (MQ=255)
                                                                  aTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCAcc                                                                                >  1:134667/1‑150 (MQ=255)
                                                                   tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCt                                                                               <  2:610380/150‑1 (MQ=255)
                                                                   tttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCt                                                                               >  2:1169894/1‑150 (MQ=255)
                                                                    ttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCtt                                                                              <  1:1639704/150‑1 (MQ=255)
                                                                     tGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTa                                                                             <  1:1507969/150‑1 (MQ=255)
                                                                        tgtgTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAtt                                                                          <  2:13077/150‑1 (MQ=255)
                                                                        tgtgTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAtt                                                                          >  1:351960/1‑150 (MQ=255)
                                                                         ttgtAAAACGAGGTGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCTGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTa                                                                         <  2:770100/149‑1 (MQ=255)
                                                                         gtgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTa                                                                         >  1:830348/1‑150 (MQ=255)
                                                                          tgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAAGTCTTACAGGTATTACACATTAAGTAc                                                                        >  1:900227/1‑150 (MQ=255)
                                                                           gtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACt                                                                       >  1:1247432/1‑150 (MQ=255)
                                                                           gtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACt                                                                       >  2:629368/1‑150 (MQ=255)
                                                                            tAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACtt                                                                      <  1:1019526/150‑1 (MQ=255)
                                                                                  gatgggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAAtt                                                                <  2:807194/147‑1 (MQ=255)
                                                                                   aGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTc                                                               <  2:1518187/150‑1 (MQ=255)
                                                                                   aGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAAAATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAGTACTTCAAAttg                                                               >  1:1498180/1‑149 (MQ=255)
                                                                                   aGGGGTTTTCCGCAGGCAGGAGAGCAAAATGACACGCTCTGTTAGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACATTAATTACTTCAAATTc                                                               >  2:1639092/1‑150 (MQ=255)
                                                                                    ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCt                                                              >  2:419187/1‑150 (MQ=255)
                                                                                       gTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACt                                                           >  2:331118/1‑150 (MQ=255)
                                                                                       gTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACt                                                           >  2:1491940/1‑150 (MQ=255)
                                                                                           tCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgc                                                       <  1:1187333/150‑1 (MQ=255)
                                                                                           tCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgc                                                       <  2:1647865/150‑1 (MQ=255)
                                                                                               caggcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCgccgcc                                                   >  2:1653444/1‑150 (MQ=255)
                                                                                                 ggcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCg                                                 <  1:373030/150‑1 (MQ=255)
                                                                                                 ggcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCg                                                 >  2:166162/1‑150 (MQ=255)
                                                                                                    aggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCg                                              <  1:103164/150‑1 (MQ=255)
                                                                                                     ggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGt                                             >  1:1383998/1‑150 (MQ=255)
                                                                                                         agCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                         <  1:948812/150‑1 (MQ=255)
                                                                                                         agCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                         <  2:1276254/150‑1 (MQ=255)
                                                                                                         agCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                         <  2:1247432/150‑1 (MQ=255)
                                                                                                           cAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                                       <  2:128397/150‑1 (MQ=255)
                                                                                                            aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatg                                      <  1:1561587/150‑1 (MQ=255)
                                                                                                             aaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatga                                     <  1:1341516/150‑1 (MQ=255)
                                                                                                                tGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAAAACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTAATTCAAATTCTACTTCGACGCCAGGCATTTTTATGACca                                  >  2:260931/1‑150 (MQ=255)
                                                                                                                 gCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcac                                 >  2:1526828/1‑150 (MQ=255)
                                                                                                                 gCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATATTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcaa                                 >  1:161817/1‑149 (MQ=255)
                                                                                                                  ccACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcaca                                >  1:647302/1‑150 (MQ=255)
                                                                                                                    cctctctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCacaac                              <  2:1347589/147‑1 (MQ=255)
                                                                                                                       ctctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatc                           >  2:620214/1‑150 (MQ=255)
                                                                                                                       ctctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGACTCCCGGCGTTTTGATGACCACAAcata                           >  1:250821/1‑149 (MQ=255)
                                                                                                                          tGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACatcatc                        >  2:920914/1‑150 (MQ=255)
                                                                                                                           gTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtct                       <  2:351960/150‑1 (MQ=255)
                                                                                                                           gTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtct                       >  2:507170/1‑150 (MQ=255)
                                                                                                                             tCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc                     >  2:279163/1‑150 (MQ=255)
                                                                                                                                ttGATAAAAGGCCGCATAGCGGCCTTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatcaatgttaa                   >  2:203846/1‑143 (MQ=255)
                                                                                                                                  ggaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATcttcttctt                >  1:69020/4‑150 (MQ=255)
                                                                                                                                  gATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATcttcttctt                <  2:798008/150‑1 (MQ=255)
                                                                                                                                  gATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATcttcttctt                >  1:916810/1‑150 (MQ=255)
                                                                                                                                   gaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt               >  1:1379357/3‑150 (MQ=255)
                                                                                                                                   aTAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt               <  2:852117/150‑1 (MQ=255)
                                                                                                                                   aTAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt               <  1:1635925/150‑1 (MQ=255)
                                                                                                                                    tAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg              >  2:1627387/1‑150 (MQ=255)
                                                                                                                                    tAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg              >  1:1364023/1‑150 (MQ=255)
                                                                                                                                      aaaGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg              >  1:1248538/1‑148 (MQ=255)
                                                                                                                                      aaaGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg              <  2:1248538/148‑1 (MQ=255)
                                                                                                                                        aGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGa          <  2:69020/150‑1 (MQ=255)
                                                                                                                                        aGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACATCAAATTCAACTTCGCCGCCCGGCGTTTACATGACCACAACATCATCTTCTTCTTTGccta          >  1:159326/1‑148 (MQ=255)
                                                                                                                                            cGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAg      >  2:1712623/1‑150 (MQ=255)
                                                                                                                                              cATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGc    <  2:1383998/150‑1 (MQ=255)
                                                                                                                                              cATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGc    <  1:331118/150‑1 (MQ=255)
                                                                                                                                                tAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCa  <  1:573505/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGCCA  >  NC_000913/3328039‑3328332

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