Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 59,301 A→C 6.0% intergenic (+22/+386) 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,3010AC 6.0% 324.7 / 7.8 122intergenic (+22/+386)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (62/52);  new base C (4/4);  total (66/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 124 total aligned reads displayed.

TAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGCCC  >  NC_000913/59152‑59449
                                                                                                                                                     |                                                                                                                                                    
tAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTTCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGaa                                                                                                                                                      <  1:425896/150‑1 (MQ=255)
   ttattattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGACGGGGCTGGAACGc                                                                                                                                                   >  2:283533/1‑150 (MQ=255)
   ttattattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGc                                                                                                                                                   <  2:1454988/150‑1 (MQ=255)
     attattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGccc                                                                                                                                                 <  2:437778/150‑1 (MQ=255)
      ttattaATAAAATCCCCCCCAAAAAAATTGATACTAAAGGCGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCGTAACAAATTACTCACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTGGAAGGGGCTGGAACGcccc                                                                                                                                                <  2:1637333/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                                                                                                                                            >  1:619823/1‑150 (MQ=255)
          taataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCa                                                                                                                                            >  2:1540910/1‑150 (MQ=255)
          taataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAAAGCCCCTTCa                                                                                                                                            >  1:1088881/1‑150 (MQ=255)
            ataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                          >  1:1048995/1‑150 (MQ=255)
            ataaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAAc                                                                                                                                          >  1:576486/1‑150 (MQ=255)
             taaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAAcc                                                                                                                                         <  1:1029665/150‑1 (MQ=255)
             taaAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCGTTCCAGCCCCTTCAAcc                                                                                                                                         <  2:1649839/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)
                         cAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCg                                                                                                                             >  2:1452651/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)
                                       cTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTcac                                                                                                               <  2:921262/150‑1 (MQ=255)
                                        aaaaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGTAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTcaca                                                                                                              <  1:967658/149‑1 (MQ=255)
                                         aaaGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACaa                                                                                                             >  2:1053659/1‑150 (MQ=255)
                                            ggAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTa                                                                                                          >  2:1392929/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)
                                                 tGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCACTTCAACCTTAGCAGTAGCGTCGGATGATTTCACAATTAGAAAg                                                                                                     >  2:1084098/1‑150 (MQ=255)
                                                  gCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGa                                                                                                    >  2:953877/1‑150 (MQ=255)
                                                     aCTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGAcat                                                                                                 >  1:556273/1‑148 (MQ=255)
                                                       tGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGc                                                                                               >  2:595160/1‑150 (MQ=255)
                                                          agaagTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatg                                                                                            >  2:1729502/1‑150 (MQ=255)
                                                             agTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatgatg                                                                                         >  2:1166522/1‑150 (MQ=255)
                                                             agTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCatgatg                                                                                         <  2:1336752/150‑1 (MQ=255)
                                                                tGCTACTAAAAAAGTACAGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTTGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGc                                                                                      <  1:825415/150‑1 (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)
                                                                           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)
                                                                                    tgAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAATACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACg                                                                  >  1:1376037/1‑150 (MQ=255)
                                                                                     gAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGc                                                                 >  2:1191344/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)
                                                                                                 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                                             >  1:1624155/1‑150 (MQ=255)
                                                                                                         aaTTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACg                                             >  2:1535959/1‑150 (MQ=255)
                                                                                                             tttGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATc                                         <  2:1427259/149‑1 (MQ=255)
                                                                                                             attGCAAGAGCCAGAATAAGCGAAGGTTGTAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTATCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATc                                         <  2:1368820/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                                      >  2:731706/1‑150 (MQ=255)
                                                                                                                gAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc                                      <  1:1452651/150‑1 (MQ=255)
                                                                                                                gAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGAGGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCAcc                                      >  2:1362170/1‑150 (MQ=255)
                                                                                                                 aaTGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGGGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCg                                     <  1:949884/150‑1 (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)
                                                                                                                          gAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCa                            >  2:1319134/1‑150 (MQ=255)
                                                                                                                           aaTAAGCTACGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCTGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCAt                           <  1:1084098/150‑1 (MQ=255)
                                                                                                                            aTAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATc                          >  2:552781/1‑150 (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                   <  1:541275/150‑1 (MQ=255)
                                                                                                                                   aaGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAAcc                   <  2:940699/150‑1 (MQ=255)
                                                                                                                                     ggTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAg                 >  1:616298/1‑150 (MQ=255)
                                                                                                                                     ggTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAg                 >  2:392699/1‑150 (MQ=255)
                                                                                                                                       ttcaaGGGGCTGGAACGCCCCGTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGAc               <  1:1551727/147‑1 (MQ=255)
                                                                                                                                          aaGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGAGCACCGAGATGGGCATCGCCAACCAGACTGc            <  1:710818/150‑1 (MQ=255)
                                                                                                                                           aGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCt           >  1:1698118/1‑150 (MQ=255)
                                                                                                                                           aGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCt           >  1:859383/1‑150 (MQ=255)
                                                                                                                                            ggggCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGTCACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTa          >  2:582368/1‑150 (MQ=255)
                                                                                                                                                cTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTa      >  2:572313/1‑150 (MQ=255)
                                                                                                                                                 tGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAg     >  2:432190/1‑150 (MQ=255)
                                                                                                                                                  ggAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGc    <  2:1532367/150‑1 (MQ=255)
                                                                                                                                                   gAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGcc   <  1:1729502/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTGCTACTAAAAAAGTACTGCTGAACAAATTACTGACAACGCAATTATTGAATGAGCCAGAATAAGCTAAGGTTGAAGGGGCTGGAACGCCCCTTCAACCTTAGCAGTAGCGTGGGATGATTTCACAATTAGAAAGACCTGCATGATGAGCTAGAGAAGAGGCTAGTGACGCAAGGCGTCGTGCAGGACACGGATCACCGAGATGGGCATCGCCAACCAGACTGCTAATTAGCCC  >  NC_000913/59152‑59449

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: