Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 59,302 C→A 5.2% intergenic (+23/+385) yabP → / ← rluA pseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,3020CA 5.2% 390.1 / 5.9 123intergenic (+23/+385)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (62/54);  new base A (4/3);  total (66/57)
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 124 total aligned reads displayed.

TTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGCCCAT  >  NC_000913/59155‑59451
                                                                                                                                                   |                                                                                                                                                     
ttattattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGACGGGGCTGGAACGc                                                                                                                                                     >  2:283533/1‑150 (MQ=255)
ttattattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGc                                                                                                                                                     <  2:1454988/150‑1 (MQ=255)
  attattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGccc                                                                                                                                                   <  2:437778/150‑1 (MQ=255)
   ttattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGcccc                                                                                                                                                  <  1:1257339/150‑1 (MQ=255)
    tattaATAAAATCCACCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCt                                                                                                                                                 <  2:295136/150‑1 (MQ=255)
      ttaATAAAATCCCCCCCAAAAAAATTTATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTc                                                                                                                                               <  2:1395687/150‑1 (MQ=255)
       taataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCa                                                                                                                                              >  2:1540910/1‑150 (MQ=255)
       taataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCa                                                                                                                                              >  1:619823/1‑150 (MQ=255)
       taataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCa                                                                                                                                              <  1:489315/150‑1 (MQ=255)
       taataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAAAGCCCCTTCa                                                                                                                                              >  1:1088881/1‑150 (MQ=255)
         ataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                            >  1:1048995/1‑150 (MQ=255)
         ataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                            >  1:576486/1‑150 (MQ=255)
          taaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAAcc                                                                                                                                           <  2:1649839/150‑1 (MQ=255)
          taaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAAcc                                                                                                                                           <  1:1029665/150‑1 (MQ=255)
              aTCCCCCACAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAg                                                                                                                                       <  1:801804/150‑1 (MQ=255)
               tCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGc                                                                                                                                      >  2:1489523/1‑150 (MQ=255)
                cccccccAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCa                                                                                                                                     <  2:655453/150‑1 (MQ=255)
                    cccAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAg                                                                                                                                 >  1:679875/1‑150 (MQ=255)
                     ccAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGc                                                                                                                                >  2:520554/1‑150 (MQ=255)
                      cAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCg                                                                                                                               >  1:940699/1‑150 (MQ=255)
                       aaaaaaaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAAcc                                                                                                                                           >  2:282233/1‑137 (MQ=39)
                       aaaaaaaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAAcc                                                                                                                                           <  1:282233/137‑1 (MQ=39)
                            aaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTTAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGgat                                                                                                                         <  2:1534977/150‑1 (MQ=255)
                            aaTTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGgat                                                                                                                         <  1:47543/150‑1 (MQ=255)
                                gATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGAtt                                                                                                                     >  2:56461/1‑150 (MQ=255)
                                 aTACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGAttt                                                                                                                    >  2:1719303/1‑150 (MQ=255)
                                    cTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTcac                                                                                                                 <  2:921262/150‑1 (MQ=255)
                                     aaaaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGTAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTcaca                                                                                                                <  1:967658/149‑1 (MQ=255)
                                      aaaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACaa                                                                                                               >  2:1053659/1‑150 (MQ=255)
                                       aaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAAt                                                                                                              >  2:1587953/1‑150 (MQ=255)
                                            tttGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGATCAAATTTCTGACAACGGAATTATTGAATGAGCCAGAATAAGCTAATGATGAAGGGGCTTGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCTCAATTAGaa                                                                                                         <  1:256035/149‑1 (MQ=255)
                                              tGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTTAACTAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAATGTTGAAGGTGCTGGAACGACCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAg                                                                                                       <  1:1217245/150‑1 (MQ=255)
                                              tGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAg                                                                                                       <  1:204781/150‑1 (MQ=255)
                                               gCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGa                                                                                                      >  2:953877/1‑150 (MQ=255)
                                                  aCTGAAGAATTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCt                                                                                                   <  2:1061926/150‑1 (MQ=255)
                                                  aCTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGAcat                                                                                                   >  1:556273/1‑148 (MQ=255)
                                                    tGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGc                                                                                                 >  2:595160/1‑150 (MQ=255)
                                                          agTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatgatg                                                                                           <  2:1336752/150‑1 (MQ=255)
                                                          agTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatgatg                                                                                           >  2:1166522/1‑150 (MQ=255)
                                                               ctactaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTa                                                                                      <  2:147087/150‑1 (MQ=255)
                                                                 actaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTaga                                                                                    <  2:1653895/150‑1 (MQ=255)
                                                                  ctaAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTagag                                                                                   >  1:966258/1‑150 (MQ=255)
                                                                   taAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAgaga                                                                                  <  1:832849/150‑1 (MQ=255)
                                                                     aaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGCCCTGCATGATGAGCTAGagaag                                                                                >  1:1126590/1‑150 (MQ=255)
                                                                     aaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGagaat                                                                                >  1:1495433/1‑149 (MQ=255)
                                                                     aaaaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATCAGCTAGagaag                                                                                >  1:655161/1‑150 (MQ=255)
                                                                       aaaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAagag                                                                              <  2:1370676/150‑1 (MQ=255)
                                                                        aaGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAgg                                                                             <  2:650690/150‑1 (MQ=255)
                                                                         aGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGc                                                                            >  1:39270/1‑150 (MQ=255)
                                                                            aCTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAg                                                                         >  2:1409018/1‑150 (MQ=255)
                                                                               gctgAACAAATTACTGACAATGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGa                                                                      >  1:1462710/1‑150 (MQ=255)
                                                                                 tgAACAAATTCCTGACATCGCAATTATTGAATGAGCAAGATTAAGCTAAGGTTGAAGTGGATGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACg                                                                    <  1:80846/150‑1 (MQ=255)
                                                                                 tgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAATACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACg                                                                    >  1:1376037/1‑150 (MQ=255)
                                                                                  gAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGc                                                                   >  2:1191344/1‑150 (MQ=255)
                                                                                  gAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGc                                                                   >  2:1299217/1‑150 (MQ=255)
                                                                                    aCAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCATCAACCTTAGCAGTAGCGTTGGATGATTTCACAATTAGAAAGACGTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCaa                                                                 >  2:847740/1‑150 (MQ=255)
                                                                                     cAAATTACTGACAACGCAATTATTGAATGAGCCTGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTCCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAg                                                                >  2:1116604/1‑150 (MQ=255)
                                                                                     aaaaTTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAg                                                                <  1:592508/149‑1 (MQ=255)
                                                                                       aaTTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTTACGCAAGGc                                                              >  2:848400/1‑150 (MQ=255)
                                                                                       aaTTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGc                                                              <  2:403979/150‑1 (MQ=255)
                                                                                       aaTTACTGAAAACGCAATTTTTGAATGAGCGAGATTATGCTAAGGTTGAAGGGGCTGGAACGCCCGTTCTACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGc                                                              <  1:739892/150‑1 (MQ=255)
                                                                                            cTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGcgtcgt                                                         >  2:1451869/1‑150 (MQ=255)
                                                                                              gACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGc                                                       <  1:1489523/150‑1 (MQ=255)
                                                                                              gACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGc                                                       >  1:1696605/1‑150 (MQ=255)
                                                                                               aCAACGCAATTTTTGTATGAGCCAGTATAAGCTAAGGTAGAAGGGGCTGGAACGCCCGTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCTAGGCGTCGTGCa                                                      <  1:1630297/150‑1 (MQ=255)
                                                                                                    gCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGaca                                                 >  1:769828/1‑150 (MQ=255)
                                                                                                     cAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGacac                                                >  1:1499438/1‑150 (MQ=255)
                                                                                                      aaTTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACg                                               >  2:1535959/1‑150 (MQ=255)
                                                                                                      aaTTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACg                                               >  1:1624155/1‑150 (MQ=255)
                                                                                                       attattGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGACCGCCCCTTCAACCTTAGCAGGAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGCCGCAAGGCGTCGGGCAGGACCCgg                                              >  2:329504/1‑150 (MQ=255)
                                                                                                          tttGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATc                                           <  2:1427259/149‑1 (MQ=255)
                                                                                                          attGCAAGAGCCAGAATAAGCGAAGGTTGTAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTATCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATc                                           <  2:1368820/150‑1 (MQ=255)
                                                                                                          attGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATc                                           <  1:683130/150‑1 (MQ=255)
                                                                                                           ttGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCa                                          >  2:1230275/1‑150 (MQ=255)
                                                                                                           ttGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCa                                          >  1:1482511/1‑150 (MQ=255)
                                                                                                            tGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAc                                         <  1:1540910/150‑1 (MQ=255)
                                                                                                             gAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc                                        <  1:1452651/150‑1 (MQ=255)
                                                                                                             gAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc                                        >  2:731706/1‑150 (MQ=255)
                                                                                                             gAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGAGGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc                                        >  2:1362170/1‑150 (MQ=255)
                                                                                                                  aGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGAt                                   >  1:937340/1‑150 (MQ=255)
                                                                                                                   cccACAATACGCTACGGTTCACGGGGCTGCACCGCCCCTTCAACATTAGCAGTAGCGTGGGATTATTCCACAATAACAACCACCTCCAGGTGGCGCTAGACAAGCGGCTAGGGACGCAAGGCGTCGGGCAGGACACGGATCACCGAGATg                                  <  2:949409/149‑1 (MQ=255)
                                                                                                                    ccAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATgg                                 >  2:1644612/1‑150 (MQ=255)
                                                                                                                      aGAATAAGCTAAGGTTGAAGGGGCAGGAACGCCGCTTCAACCTGTGCAGTAGCGTGGGATGATTTCACGATTAGAAAGCCCTGCTTGATGAGCTAGAGAAGAGGCTAGTGACGCACGGCGTCGTGCCGGACACGGATCACCGAGATGGGc                               >  1:288747/1‑150 (MQ=255)
                                                                                                                        aaTAAGCTACGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCTGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCAt                             <  1:1084098/150‑1 (MQ=255)
                                                                                                                           aaGCTAAGGTTGATGGGGCTGGAACGCCCCTTCAACATTAGTAGTTGCGTTGGATGATTTCACAATTAGAAAGACCTGCTTGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGc                          <  2:556273/150‑1 (MQ=255)
                                                                                                                             gCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCa                        >  2:248297/1‑150 (MQ=39)
                                                                                                                              cTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAGCCTTAGCAGTAGCGTGGGAGGATTTCACAATTAGAAAGACCTGCATGATGAGGTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCaa                       >  1:1241604/1‑150 (MQ=38)
                                                                                                                               tAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAAc                      <  1:186885/150‑1 (MQ=255)
                                                                                                                                aaGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAAcc                     <  2:679875/150‑1 (MQ=255)
                                                                                                                                aaGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAAcc                     <  2:940699/150‑1 (MQ=255)
                                                                                                                                    ttcaaGGGGCTGGAACGCCCCGTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGAc                 <  1:1551727/147‑1 (MQ=255)
                                                                                                                                        aGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCt             >  1:859383/1‑150 (MQ=255)
                                                                                                                                        aGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCt             <  1:612079/150‑1 (MQ=255)
                                                                                                                                        aGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCt             >  1:1698118/1‑150 (MQ=255)
                                                                                                                                             cTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTa        >  2:572313/1‑150 (MQ=255)
                                                                                                                                              tGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAg       >  2:432190/1‑150 (MQ=255)
                                                                                                                                               ggAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGc      <  2:576486/150‑1 (MQ=255)
                                                                                                                                                gAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGcc     <  1:1729502/150‑1 (MQ=255)
                                                                                                                                                 aaCGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGccc    <  2:1501471/150‑1 (MQ=255)
                                                                                                                                                   cGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGCCCAt  <  2:619823/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGCCCAT  >  NC_000913/59155‑59451

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: