Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,183 T→A 6.9% 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 6.9% 223.7 / 4.9 81intergenic (+100/+56)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base T (42/33);  new base A (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 = 9.88e-01

TGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGC  >  NC_000913/3328034‑3328330
                                                                                                                                                     |                                                                                                                                                   
tGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCAt                                                                                                                                                     >  2:567311/1‑150 (MQ=255)
     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:1261175/150‑1 (MQ=255)
                                                                                                                                gTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtct                     <  2:209892/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)
                                                                                                                                                     |                                                                                                                                                   
TGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATAGCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGC  >  NC_000913/3328034‑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: