Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,976,120 C→T 5.7% noncoding (70/88 nt) omrA ← small regulatory RNA OmrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,976,1200CT 5.7% 326.9 / 4.3 104noncoding (70/88 nt)omrAsmall regulatory RNA OmrA
Reads supporting (aligned to +/- strand):  ref base C (47/50);  new base T (3/3);  total (51/53)
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

Only 100 of 108 total aligned reads displayed.

CGATACAAAACACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTA  >  NC_000913/2975971‑2976269
                                                                                                                                                     |                                                                                                                                                     
cgaTACAAAACACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATcc                                                                                                                                                       >  1:1533296/1‑150 (MQ=255)
 gaTACAAAACACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCcg                                                                                                                                                      <  2:974465/150‑1 (MQ=255)
     cAAAACACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCa                                                                                                                                                  >  1:164382/1‑150 (MQ=255)
     cAAAACACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCa                                                                                                                                                  <  1:511350/150‑1 (MQ=255)
       aaaCACGCCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAgg                                                                                                                                                <  1:840687/150‑1 (MQ=255)
        aaCACTTCACAAATTTCGAAAAAATCTAAAAAGCATTACTAGTCGCTTTGCTCTGATTTATTAAATGCGATTGGAGGGAGATTACACGAGATATAGAACTCGAGCTACAGTAAATTAGGTGCGAAAAAAAACCTTGGCATCCGCGGAGGt                                                                                                                                               <  2:255348/150‑1 (MQ=255)
        aaCACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGt                                                                                                                                               <  2:491565/150‑1 (MQ=255)
         acacGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCGGATGCGCAGGtt                                                                                                                                              <  1:1189749/150‑1 (MQ=255)
         acacGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGtt                                                                                                                                              <  1:214277/150‑1 (MQ=255)
           acGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCGGATGCGCAGGtttt                                                                                                                                            <  2:469063/150‑3 (MQ=255)
            cGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGt                                                                                                                                           <  1:655230/150‑1 (MQ=255)
            cGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGt                                                                                                                                           <  1:903871/150‑1 (MQ=255)
             ggCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCGGATGCGCAGGtttttt                                                                                                                                          <  1:920584/150‑5 (MQ=255)
              gCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGc                                                                                                                                         <  2:1229058/150‑1 (MQ=255)
                 caAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAag                                                                                                                                      <  2:182364/150‑1 (MQ=255)
                 aaaaTTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAag                                                                                                                                      <  1:606809/149‑1 (MQ=255)
                  aaaTTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGGGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAaga                                                                                                                                     <  1:19124/150‑1 (MQ=255)
                   aaTTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAagag                                                                                                                                    >  1:637283/1‑150 (MQ=255)
                      ttCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACa                                                                                                                                 >  1:52398/1‑150 (MQ=255)
                      ttCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACa                                                                                                                                 >  1:1070232/1‑150 (MQ=255)
                        cGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGATCTACATTAAATTAGGTGCGAAAAAAAACCTTCTCATCCGCGCAGGTTGGTGCAAGAGACAgg                                                                                                                               <  2:636661/150‑1 (MQ=255)
                        cGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAgg                                                                                                                               <  1:1269034/150‑1 (MQ=255)
                          aaaaaaGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGt                                                                                                                             <  2:709801/150‑1 (MQ=255)
                           aaaaaGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTa                                                                                                                            >  2:715218/1‑150 (MQ=255)
                            aaaaGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGATATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTAc                                                                                                                           <  2:687437/150‑1 (MQ=255)
                              aaGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGa                                                                                                                         <  2:332079/150‑1 (MQ=255)
                                gCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAag                                                                                                                       >  2:801830/1‑150 (MQ=255)
                                 cTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAaga                                                                                                                      <  1:1415285/150‑1 (MQ=255)
                                 cTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAaga                                                                                                                      >  1:742026/1‑150 (MQ=255)
                                  tAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACgctctt                                                                                                                     >  1:626094/1‑145 (MQ=255)
                                  tAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAagag                                                                                                                     <  2:164382/150‑1 (MQ=255)
                                   aaaaaGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGc                                                                                                                    >  2:1170971/1‑150 (MQ=255)
                                     aaaGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGt                                                                                                                  <  2:235086/150‑1 (MQ=255)
                                        gCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTTGTGCAAGAGACAGGGTACGAAGAGCGTAcc                                                                                                               <  1:1408000/150‑1 (MQ=255)
                                          aTTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCGAGAGACGGGGTACGAAGcgcggggcga                                                                                                             >  1:287890/1‑144 (MQ=255)
                                            tCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGATACGAAGAGCGTACCGaat                                                                                                           >  1:1280917/1‑150 (MQ=255)
                                              cAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGTGAGATTACACGAGATAAAGAAAGCTAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGATGGTGCAAGAGACAGGCTACGAAGAGCGTACCGaataa                                                                                                         >  2:741872/1‑150 (MQ=255)
                                                   cccTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCa                                                                                                    >  2:105397/1‑150 (MQ=255)
                                                    ccTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCAc                                                                                                   <  2:579040/150‑1 (MQ=255)
                                                      tttGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACca                                                                                                 <  2:734137/150‑1 (MQ=255)
                                                          ctctGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACcaatca                                                                                             >  1:92132/1‑150 (MQ=255)
                                                          ctctGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGATACGAAGAGCGTACCGAATAATCTCACcaatca                                                                                             <  2:1280917/150‑1 (MQ=255)
                                                            ctGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACcaatcaat                                                                                           <  2:1070232/150‑1 (MQ=255)
                                                             tGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATa                                                                                          <  2:1010098/150‑1 (MQ=255)
                                                             tGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATa                                                                                          <  2:1004968/150‑1 (MQ=255)
                                                              gATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATAc                                                                                         <  2:304685/150‑1 (MQ=255)
                                                               aTTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATAcc                                                                                        <  1:681852/150‑1 (MQ=255)
                                                               aTTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATAcc                                                                                        >  1:1537665/1‑150 (MQ=255)
                                                                  cAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACctct                                                                                     >  2:1410859/1‑150 (MQ=255)
                                                                   aGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTg                                                                                    >  1:723210/1‑150 (MQ=255)
                                                                      aaaaGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGa                                                                                 >  2:1182683/1‑150 (MQ=255)
                                                                       aaaGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGAt                                                                                >  2:446563/1‑150 (MQ=255)
                                                                         aGCGAATGGACGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAAATACGTGCGAAAAAAAACCTGCGCATCCCCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCt                                                                              >  2:1404745/1‑150 (MQ=255)
                                                                          gCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCtt                                                                             <  2:742026/150‑1 (MQ=255)
                                                                                ggagggagATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATtgt                                                                       <  2:1073730/150‑1 (MQ=255)
                                                                                 gagggagATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATtgtg                                                                      <  1:543552/150‑1 (MQ=255)
                                                                                  agggagATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCGGGGATCTTGATTGTgg                                                                     <  2:1017546/150‑1 (MQ=255)
                                                                                   gggagATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGt                                                                    <  2:1533296/150‑1 (MQ=255)
                                                                                    ggagATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTc                                                                   <  1:264876/150‑1 (MQ=255)
                                                                                       gaTTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGc                                                                >  1:191565/1‑150 (MQ=255)
                                                                                        aTTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCa                                                               >  1:676846/1‑150 (MQ=255)
                                                                                        aTTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGAAAGGGTACGAAGAGGGTAACGAATAATCTAACCAATCAATACATCTGGGATCTTGATTGTGGTCTgaa                                                               >  1:775171/1‑148 (MQ=255)
                                                                                          tACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCacg                                                             <  1:1354049/150‑1 (MQ=255)
                                                                                           acacGAGATAAAGAACGCGAGCGCCAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATATCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCacga                                                            >  2:1118833/1‑150 (MQ=255)
                                                                                           acacGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCacga                                                            >  1:1422801/1‑150 (MQ=255)
                                                                                            cacGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCacgac                                                           >  1:1309377/1‑150 (MQ=255)
                                                                                            cacGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCacgac                                                           <  2:203457/150‑1 (MQ=255)
                                                                                              cGAGATAAAGAACTCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCAGAGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCAcgacga                                                         >  2:63717/1‑150 (MQ=255)
                                                                                                   tAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCtt                                                    >  1:639095/1‑150 (MQ=255)
                                                                                                     aaGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTg                                                                            <  2:955178/124‑1 (MQ=255)
                                                                                                     aaGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTg                                                                            >  1:955178/1‑124 (MQ=255)
                                                                                                         aCGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAg                                              <  2:556791/150‑1 (MQ=255)
                                                                                                         aCGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATAAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAg                                              <  2:1093576/150‑1 (MQ=255)
                                                                                                                   cAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGAATCTTCGCCAGCGAGAAAACg                                    >  1:543955/1‑150 (MQ=255)
                                                                                                                     gTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCa                                  <  1:1170971/150‑1 (MQ=255)
                                                                                                                     gTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCa                                  >  2:475025/1‑150 (MQ=255)
                                                                                                                      tAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCaa                                 >  1:395693/1‑150 (MQ=255)
                                                                                                                         aTTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGTGATCTTGATTTTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAgg                              <  2:838454/150‑1 (MQ=255)
                                                                                                                         aTTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACAACTCTTCGCCAGCGAGAAAACGCAAAgg                              >  1:417823/1‑150 (MQ=255)
                                                                                                                               tGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACctct                                                                                     <  1:174896/89‑1 (MQ=255)
                                                                                                                               tGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACctct                                                                                     >  2:174896/1‑89 (MQ=255)
                                                                                                                                 cGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAgg                      >  2:1502481/1‑150 (MQ=255)
                                                                                                                                 cGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAgg                      >  1:999454/1‑150 (MQ=255)
                                                                                                                                 cGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAgg                      <  2:52398/150‑1 (MQ=255)
                                                                                                                                 cGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCGTCGCCAGCGAGAAAACGCAAAGGAATGAAgg                      >  1:829931/1‑150 (MQ=255)
                                                                                                                                    gcaccaaCCTGCGCGGATGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACAAATCAATACCTATGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGaa                   >  2:984676/6‑150 (MQ=255)
                                                                                                                                    aaaaaaaCCTACGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTACAATGTGGAAACCACGACCACACATCGCCACCCAGACAAACCAAAGGAAGGAACGGaa                   >  2:564548/1‑150 (MQ=255)
                                                                                                                                     caccaaCCTGCGCGGATGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCAGGACGACTCTTCGCCAGCGAGAAAAAGCAAAGGAATGAAGGGaaa                  >  1:518386/5‑150 (MQ=255)
                                                                                                                                      caaacCCCGCGCACCCGCGCAGGTTGGCGCACGCCCCACGGCCCGACGCCCGCACCGAATAATCTCACCAATCAAAACCTCCCGGATCTTGATTGTGGCCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAAt                 <  2:287890/149‑1 (MQ=255)
                                                                                                                                      aaaaaCCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAAt                 <  1:845461/150‑1 (MQ=255)
                                                                                                                                      aaaaaCCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAAt                 <  1:1327079/150‑1 (MQ=255)
                                                                                                                                         aaCCTGCGCGGATGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCa              >  1:1151884/1‑150 (MQ=255)
                                                                                                                                           ccTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAAc            >  2:196407/1‑150 (MQ=255)
                                                                                                                                           ccTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAAc            >  1:499799/1‑150 (MQ=255)
                                                                                                                                            cTGCGCATCCGCGCAGGTTGGTGCTAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGTTCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGGGAGTAAACGCAAAGGAATGAAGGGAAATGCAACg           <  1:974179/150‑1 (MQ=255)
                                                                                                                                                 cATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGgtg      >  2:1240972/1‑150 (MQ=255)
                                                                                                                                                 cATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGgtg      >  1:1394941/1‑150 (MQ=255)
                                                                                                                                                  aTCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGgtgt     >  2:831899/1‑150 (MQ=255)
                                                                                                                                                    ccGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGgtgtgt   >  1:370818/1‑150 (MQ=255)
                                                                                                                                                     cgcgCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTa  >  1:856629/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATACAAAACACGGCACAAATTTCGAAAAAAGCTAAAAAGCATTCCAACTCCCTTTGCTCTGATTCAGTAAAAGCGAATGGAGGGAGATTACACGAGATAAAGAACGCGAGCGACAGTAAATTAGGTGCGAAAAAAAACCTGCGCATCCGCGCAGGTTGGTGCAAGAGACAGGGTACGAAGAGCGTACCGAATAATCTCACCAATCAATACCTCTGGGATCTTGATTGTGGTCTGCACGACGACTCTTCGCCAGCGAGAAAACGCAAAGGAATGAAGGGAAATGCAACGAGGTGTGTA  >  NC_000913/2975971‑2976269

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: