Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,648 A→T 8.5% intergenic (+46/‑104) ybgE → / → ybgC cyd operon protein YbgE/esterase/thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,6480AT 8.5% 206.8 / 5.6 69intergenic (+46/‑104)ybgE/ybgCcyd operon protein YbgE/esterase/thioesterase
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base T (3/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAG  >  NC_000913/774502‑774797
                                                                                                                                                  |                                                                                                                                                     
gCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAAttt                                                                                                                                                    >  2:50458/1‑150 (MQ=255)
    gggCTTTCGTCCGCCGACGGTTCTTTGGCCAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAACTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTAc                                                                                                                                                <  2:576720/150‑1 (MQ=255)
     ggCTTTCGTCCGCAGATGGTTCTTTGGCAAGCGATTTTTTGCGCTTTGCTTTGCGATATTGTTCTCATTTTCGGGCTGATTTTCTTCTTCTTTTAAATGAGAATTCTCTTTATAAATTTATGTGCCTGGCGCGGCCCATAATTTTTTaaa                                                                                                                                               <  2:551118/150‑3 (MQ=255)
     ggCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACg                                                                                                                                               <  2:50593/150‑1 (MQ=255)
      gCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTaaag                                                                                                                                              <  2:1383133/150‑4 (MQ=38)
       cTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCt                                                                                                                                             >  1:454854/1‑150 (MQ=255)
       cTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCt                                                                                                                                             >  2:1225437/1‑150 (MQ=255)
       cTTGCGGCCGCAGAAGGCTCTTTGGCAAGGGATGGTTTGCCCATTGCGTGCCGATATTGTGCTCATTGTCGGGCTGATTTTCTTCTGCTTTTAAATCAGAAGTCTCTTTAAAAATTTATGGGCCGGGAGCGGCCCATAATTTTTTaaaga                                                                                                                                             <  1:1280820/150‑5 (MQ=25)
         ttCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTcc                                                                                                                                           >  1:230472/1‑150 (MQ=255)
                  aGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTg                                                                                                                                  <  1:936542/150‑1 (MQ=255)
                   gAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGATCCCTTAACTTGc                                                                                                                                 >  1:391542/1‑150 (MQ=255)
                   gAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCGGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGc                                                                                                                                 <  1:927600/150‑1 (MQ=255)
                      ggTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCt                                                                                                                              >  2:410710/1‑150 (MQ=255)
                          cTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCAtt                                                                                                                          >  1:145162/1‑150 (MQ=255)
                                 aaGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAc                                                                                                                   <  1:410710/150‑1 (MQ=255)
                                  aGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAAcc                                                                                                                  <  2:602969/150‑1 (MQ=255)
                                   gggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCt                                                                                                                 >  2:987386/1‑150 (MQ=255)
                                   gggAATTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTATTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTACCTTGCCCTCATTCCCACACCt                                                                                                                 >  2:180483/1‑150 (MQ=255)
                                    ggATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTc                                                                                                                >  2:290754/1‑150 (MQ=255)
                                                 aTTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGta                                                                                                   <  1:1336441/150‑1 (MQ=255)
                                                    gcttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtag                                                                                                >  1:399434/1‑150 (MQ=255)
                                                     cttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGGTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagt                                                                                               >  1:471236/1‑150 (MQ=255)
                                                     cttgcCGATACTGTTCTCATTGTCGCGCTGATCTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCACGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagt                                                                                               <  2:182736/150‑1 (MQ=255)
                                                      ttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagta                                                                                              >  2:79591/1‑150 (MQ=255)
                                                      ttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagta                                                                                              <  1:369566/150‑1 (MQ=255)
                                                      ttgcCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTAtagta                                                                                              <  1:658837/150‑1 (MQ=255)
                                                               ttGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGttt                                                                                     <  1:409638/150‑1 (MQ=255)
                                                                 gTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTaa                                                                                   >  2:1095354/1‑150 (MQ=255)
                                                                    ctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAg                                                                                >  1:1136454/1‑150 (MQ=255)
                                                                    ctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAg                                                                                <  1:290754/150‑1 (MQ=255)
                                                                    ctcATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAAACGCGCGCGTATAGTAGCAGCGTTTAAAAg                                                                                >  1:1141448/1‑150 (MQ=255)
                                                                      cATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTa                                                                                                                                       >  1:560699/1‑93 (MQ=255)
                                                                      cATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTa                                                                                                                                       <  2:560699/93‑1 (MQ=255)
                                                                      cATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgcg                                                                                                     <  1:1144985/127‑1 (MQ=255)
                                                                      cATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATcgcgcgcg                                                                                                     >  2:1144985/1‑127 (MQ=255)
                                                                         tGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTa                                                                           >  1:247220/1‑150 (MQ=255)
                                                                          gTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTaa                                                                          >  2:178736/1‑150 (MQ=255)
                                                                                  gATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTATCTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgtt                                                                  <  1:998350/150‑1 (MQ=255)
                                                                                   aTTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTttgttg                                                                 <  1:1101191/150‑1 (MQ=255)
                                                                                    ttttCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGc                                                                >  1:1077196/1‑150 (MQ=255)
                                                                                        cttcttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACGTTTGTTGCATTa                                                            >  1:322623/1‑150 (MQ=255)
                                                                                            ttcttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCgg                                                        >  2:1066369/1‑150 (MQ=255)
                                                                                              cttTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGa                                                      <  1:1095354/150‑1 (MQ=255)
                                                                                               ttttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGAt                                                     >  2:67070/1‑150 (MQ=255)
                                                                                               ttttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGAt                                                     <  2:655182/150‑1 (MQ=255)
                                                                                               ttttAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGAt                                                     >  1:455234/1‑150 (MQ=255)
                                                                                                   aaaTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTaa                                                 >  2:531276/1‑150 (MQ=255)
                                                                                                    aaTCAGAATACTCTTTAAAAAAATATTAGCCGATACAAGCACAAAAATTTATACGAAACATTAAAATAAACTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTaaa                                                >  2:611878/1‑150 (MQ=255)
                                                                                                     aTCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAg                                               <  1:67070/150‑1 (MQ=255)
                                                                                                      tCAGAATTATCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTGTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAgc                                              <  1:427514/150‑2 (MQ=255)
                                                                                                       cAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTg                                             >  1:640301/1‑150 (MQ=255)
                                                                                                         gAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGaa                                           <  2:251597/150‑1 (MQ=255)
                                                                                                          aaTTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAAt                                          <  2:1077196/150‑1 (MQ=255)
                                                                                                                 cgtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACgg                                   >  1:339622/3‑149 (MQ=38)
                                                                                                                 cgtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGc                                   >  1:1197850/3‑150 (MQ=38)
                                                                                                                 cgtAAAAATTTATGGGCCTGGAGCGGCCCATAATTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGc                                   >  1:1205566/3‑150 (MQ=38)
                                                                                                                  ttAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCt                                  <  1:34893/150‑1 (MQ=255)
                                                                                                                     aaaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGtt                               >  2:527767/1‑150 (MQ=255)
                                                                                                                      aaaaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGttt                              <  1:987386/150‑1 (MQ=255)
                                                                                                                        aaTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCg                            >  2:612284/1‑150 (MQ=255)
                                                                                                                         aTTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGa                           >  2:411307/1‑150 (MQ=255)
                                                                                                                           tATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATg                         <  1:160249/150‑1 (MQ=255)
                                                                                                                              gggCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGTcc                      <  2:399434/150‑1 (MQ=255)
                                                                                                                              gggCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGcc                      <  2:230472/150‑1 (MQ=255)
                                                                                                                                   gCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTc                 >  2:399592/1‑150 (MQ=255)
                                                                                                                                        acgCCCATAAATATTTACGCTCCCTAAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACATTTGTTCCATTACCCGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc            <  2:156337/148‑1 (MQ=255)
                                                                                                                                        aGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTc            <  1:1225437/150‑1 (MQ=255)
                                                                                                                                          gCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCTCGTcta          >  1:169155/1‑150 (MQ=255)
                                                                                                                                             cATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTctacta       <  1:178736/150‑1 (MQ=255)
                                                                                                                                                 aaTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTGGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGaa   >  2:914429/1‑150 (MQ=255)
                                                                                                                                                  aTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAg  <  2:391542/150‑1 (MQ=255)
                                                                                                                                                  aTTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAg  >  1:514533/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCGTGGGCTTTCGTCCGCAGAAGGTTCTTTGGCAAGGGATTTTTTGCCCATTGCTTGCCGATATTGTTCTCATTGTCGGGCTGATTTTCTTCTTCTTTTAAATCAGAATTCTCTTTAAAAAATTATGGGCCGCTCCAGGCCCATAAATTTTTACGCTCCCTTAACTTGCCCTCATTCCCAAACCTCAATCGCGCGCGTATAGTAGCAGCGTTTAAAAGCTCTAACTTTTGTTGCATTACCGGGATGTAAAGTGAATACAACGCTGTTTCGATGGCCGGTTCGCGTCTACTATGAAG  >  NC_000913/774502‑774797

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