Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,184 A→T 6.9% intergenic (+101/+55) 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,1840AT 6.9% 220.5 / 4.6 81intergenic (+101/+55)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base A (42/33);  new base T (4/2);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGC  >  NC_000913/3328039‑3328330
                                                                                                                                                 |                                                                                                                                                  
cTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCgg                                                                                                                                                <  2:133033/150‑1 (MQ=255)
 tttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCTAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGc                                                                                                                                               >  2:221044/1‑150 (MQ=255)
  ttATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGcc                                                                                                                                              <  1:222473/150‑1 (MQ=255)
   tATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTTAGCGGCCt                                                                                                                                             <  2:816761/150‑1 (MQ=255)
       gCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGCGGGGTGTTCCGCAGGCAGGAGAGCAAAATGCACCGCTCTCTTCGTTGATAAAAGGCCGCATAGCGGCCttttt                                                                                                                                         >  1:209892/1‑150 (MQ=255)
       gCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCttttt                                                                                                                                         >  2:1243819/1‑150 (MQ=255)
        cccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTATGCGGCCttttat                                                                                                                                        <  1:1111928/150‑3 (MQ=255)
        cccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTc                                                                                                                                        >  1:1285001/1‑150 (MQ=255)
          cAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCATGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGCCAGCATAGCGGCATTTTTCct                                                                                                                                      >  2:1118624/1‑150 (MQ=255)
          cAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCct                                                                                                                                      >  2:128574/1‑150 (MQ=255)
           aaCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAACGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTGCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctt                                                                                                                                     <  1:990607/150‑1 (MQ=255)
           aaCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctt                                                                                                                                     <  1:726719/150‑1 (MQ=255)
           aaCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACAAATTTGCTGTGTAAAACGAGGTGTTTTCCGCAGGCAGGAGAGCAAAATTCCACGCGCTGTTCGTTGATAAAAGGCCGCTATGCGGCCttttatcaa                                                                                                                                     <  1:658844/150‑6 (MQ=255)
             cTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttc                                                                                                                                   >  1:662317/1‑150 (MQ=255)
                aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttcttt                                                                                                                                <  2:945754/150‑1 (MQ=255)
                aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCctttcttt                                                                                                                                <  2:142289/150‑1 (MQ=255)
                aaGAACGTAAAATCTCGCTGCCACGCTAAGAGTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAATGGCCGCATAGCGGCCGTTTTCctttcttt                                                                                                                                <  1:839811/150‑1 (MQ=255)
                  gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTAc                                                                                                                              <  2:332147/150‑1 (MQ=255)
                      gTAAAATCTCGCTGCCCCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTACCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                          >  2:865593/1‑150 (MQ=255)
                       tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATAc                                                                                                                         >  1:4445/1‑150 (MQ=255)
                       tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATAc                                                                                                                         >  1:327333/1‑150 (MQ=255)
                         aaaTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACAt                                                                                                                       >  1:1025366/1‑150 (MQ=255)
                               cGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACAt                                                                                                                 <  2:300034/150‑1 (MQ=255)
                                 cTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCt                                                                                                               >  2:435628/1‑150 (MQ=255)
                                     cACGCTAAGATTATCCTAAAGTTACACACATTTGCTATGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGACTTTTTACTTTCTTTACAATACATCAACATCTTGAg                                                                                                           >  2:62497/1‑150 (MQ=255)
                                      aCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCTCCTAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGt                                                                                                          >  1:793954/1‑150 (MQ=255)
                                       cGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                         >  2:1209022/1‑150 (MQ=255)
                                       cGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGATTTCCGCAGGCAGGAGAGCAAAATGCCAAGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                         >  1:1222239/1‑150 (MQ=255)
                                         cTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                       >  1:1221796/1‑150 (MQ=255)
                                             taTTATCCTAAAGTTACACACATTTGCTGTGAAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCTAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                                   <  2:1116005/149‑1 (MQ=255)
                                             gATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                                   <  2:127342/150‑1 (MQ=255)
                                                  tCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTc                                                                                              <  1:1203784/150‑1 (MQ=255)
                                                     tAAAGTTACAAACATTTGCTGTGTAAAACTAGGGGTTTTCCGCAGGCTGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCATTTCTTTACAATACATCAACATATTGAGTATTGGGGAATTCttt                                                                                           >  1:626014/1‑149 (MQ=255)
                                                      aaaGTTACAAACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTAc                                                                                          <  2:828775/150‑1 (MQ=255)
                                                       aaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACa                                                                                         <  2:956784/150‑1 (MQ=255)
                                                       aaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGAAAAAATGCCGCATAGCGGCCTTTTTCCTTTCTTTACAAGACCTCAACATCTTCAGTATTGGGTAATTCTTACa                                                                                         >  1:130146/1‑150 (MQ=255)
                                                        aggtaaccaccTTTGCTGTGGAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCGCTGTTCGTTGATAAACGGCCGCATAGCGGCCTTTTTCCTTTCTTTCCAATACATCAACATCTTGAGTATTGGGTAATTCTTACAg                                                                                        <  2:818929/143‑1 (MQ=255)
                                                         gTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAgg                                                                                       <  2:1074354/150‑1 (MQ=255)
                                                             cacacaTTTGCTGTGTAAAACTAGGGGGTTTAAGCAGGCAGTAGAGCAAAATACCACGCTCTGTTCGTTGATAAAAGGCCGCATAGAGGCAATATTCCTATCTTGACAATACATCAACATCTTGAGTAATGCTTAATTCTTACAGGTAtt                                                                                   >  1:653572/1‑150 (MQ=255)
                                                             cacacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                   >  2:93102/1‑150 (MQ=255)
                                                                  aTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGTGTAATTCTTACAGGTATTCAAcc                                                                              >  2:41282/1‑150 (MQ=255)
                                                                         gtgtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTa                                                                       >  2:120171/1‑150 (MQ=255)
                                                                           gtAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAAGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACt                                                                     <  2:362358/150‑1 (MQ=255)
                                                                             aaaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTc                                                                   >  1:811902/1‑150 (MQ=255)
                                                                             aaaaCGAGGGGTTTTCAGCAGTAAGGAGAGCAAAATGCCACGCTCTTTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTTGGTAATTCTTACAGGTATTCCACCTTAATTACTTc                                                                   <  1:769404/150‑1 (MQ=255)
                                                                              ataCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACTATTCATCAACATCTTGAGTTTTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCa                                                                  <  1:886403/148‑1 (MQ=255)
                                                                              aaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCa                                                                  <  2:344184/150‑1 (MQ=255)
                                                                                   aGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTAGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTc                                                             <  1:203612/150‑1 (MQ=255)
                                                                                        ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACtt                                                        <  2:1102300/150‑1 (MQ=255)
                                                                                        ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTACCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACtt                                                        <  2:338757/150‑1 (MQ=255)
                                                                                           tCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgc                                                     <  1:748525/150‑1 (MQ=255)
                                                                                            ccGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTcgcc                                                    >  2:203274/1‑150 (MQ=255)
                                                                                               caggcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAATGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCgccgcc                                                 <  1:88536/150‑1 (MQ=255)
                                                                                                      gagaGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGtt                                          >  1:1185163/1‑150 (MQ=255)
                                                                                                       agagCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGttt                                         >  1:839965/1‑150 (MQ=255)
                                                                                                         agCAAAATGCCAATCTATGTTATTTGATAAAAGTCCGCATAGCGTCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCTCCCGGCGTTTtg                                       <  1:184061/150‑1 (MQ=255)
                                                                                                          gCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtga                                      >  1:170834/1‑150 (MQ=255)
                                                                                                          gCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTGTCTTTACACTACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtga                                      >  1:632536/1‑150 (MQ=255)
                                                                                                           cAAAATGCCACGCTCTGTTCCTTGATAAAAGGCCTCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgat                                     <  1:917948/150‑1 (MQ=255)
                                                                                                             aaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatga                                   >  2:306193/1‑150 (MQ=255)
                                                                                                               aTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGAcc                                 <  2:1229255/150‑1 (MQ=255)
                                                                                                                  ccACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcaca                              >  2:1112510/1‑150 (MQ=255)
                                                                                                                   cACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCCCCTGAATTACTTCAAATACTAATTCGCCGCCCGGCGTTTTGATGACcacca                             >  1:563776/1‑148 (MQ=255)
                                                                                                                    aCGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGGTGACCacaaa                            >  1:345368/1‑149 (MQ=255)
                                                                                                                       ctctGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatc                         >  1:640188/1‑150 (MQ=255)
                                                                                                                          tGTTCGTTGATAAAAGGCCGCATAGCGGCGTTTTACCGTTCTTTACAGTACATCAACATCTTGAGTATTGCGTAATTCTTAAAGGTATTCCACCTTAATTAATACAAATACTACTTCGCCGCCCGGCGTTTTGATGACCACAACatcgtc                      >  1:58103/1‑150 (MQ=255)
                                                                                                                           gTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtct                     <  2:209892/150‑1 (MQ=255)
                                                                                                                           gTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtct                     <  2:1261175/150‑1 (MQ=255)
                                                                                                                             tCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc                   >  2:920012/1‑150 (MQ=255)
                                                                                                                              cGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttct                  <  2:1005368/150‑1 (MQ=255)
                                                                                                                                aaggaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttctta                >  1:45674/6‑149 (MQ=255)
                                                                                                                                aaggaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttattc                >  2:1085531/6‑150 (MQ=255)
                                                                                                                                 tGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttcttct               >  1:837171/1‑150 (MQ=255)
                                                                                                                                   gaaaaaGGCCGCTATGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt             >  1:691846/3‑150 (MQ=255)
                                                                                                                                   gaaaaaGGCCGCTATGCGGCCTTTATCCATTCTTTAAAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCGACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTGTTGAGGCCGACAACATCAGATTCTTGttt             >  2:203562/3‑150 (MQ=255)
                                                                                                                                     aaaaGGCCGCATAGTGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAATACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc           <  1:653308/150‑1 (MQ=255)
                                                                                                                                         ggCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAATATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGAt       >  2:753607/1‑150 (MQ=255)
                                                                                                                                          gCCGCATATCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTACACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATc      >  2:1297090/1‑150 (MQ=255)
                                                                                                                                          gCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATc      >  1:970650/1‑150 (MQ=255)
                                                                                                                                           ccGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCa     >  1:184558/1‑150 (MQ=255)
                                                                                                                                              cATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGc  >  1:343758/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGC  >  NC_000913/3328039‑3328330

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