Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,596,503 T→G 6.3% pseudogene (1485/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,596,5030TG 6.3% 271.7 / 6.8 91pseudogene (1485/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base G (3/3);  total (45/46)
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.49e-01

Only 100 of 101 total aligned reads displayed.

TCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTT  >  NC_000913/1596354‑1596651
                                                                                                                                                     |                                                                                                                                                    
tCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAttt                                                                                                                                                      <  2:1435452/150‑1 (MQ=255)
 cACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtc                                                                                                                                                     <  1:1398072/150‑1 (MQ=255)
  aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                                    <  1:615220/150‑1 (MQ=255)
  aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                                    <  2:1370055/150‑1 (MQ=255)
   cgcgCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcct                                                                                                                                                   >  1:1176499/1‑150 (MQ=255)
     cgcTTTAACTATACCTCCGGTGTTAATATTCAGTTCACCCGTGACTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGATATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                                 <  1:221745/150‑1 (MQ=255)
     cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                                 <  2:1354936/150‑1 (MQ=255)
     cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                                 <  2:229194/150‑1 (MQ=255)
      gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                                <  2:555897/150‑1 (MQ=255)
      gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                                <  1:1271347/150‑1 (MQ=255)
      gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                                <  1:487911/150‑1 (MQ=255)
      gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACACGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                                <  1:1178538/150‑1 (MQ=255)
        tttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACCCCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTATCGGAGGAAATGCCTCCggt                                                                                                                                              <  2:114356/150‑3 (MQ=255)
        tttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGAt                                                                                                                                              >  1:1014878/1‑150 (MQ=255)
        tttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGAt                                                                                                                                              <  2:409838/150‑1 (MQ=255)
        tttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCggt                                                                                                                                              <  1:614521/150‑3 (MQ=255)
         ttACCTCTACCTCCGGTTGTAATATTCAGTTCACCCGTGCCTAATTCTCCGACCTGTAGCATTTGTGCGGTAGTAGATGATTTCTGTAAATTCCAGAGACTGTGCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTa                                                                                                                                             <  2:679764/150‑1 (MQ=255)
              tataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATcc                                                                                                                                        <  2:639994/150‑1 (MQ=255)
              aataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATcc                                                                                                                                        <  2:1310364/149‑1 (MQ=255)
                taCCTCCGGTGGTAATATTCTGTTCACGCGTGCCTAATACACCGACCTGTAGCATTTGTGCGTTAGTAGATGATGTCTTTAAATTGCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCAt                                                                                                                                      <  2:385011/150‑1 (MQ=255)
                 aCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCAGATAAATCCATa                                                                                                                                     >  1:1435994/1‑150 (MQ=255)
                   cTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTACATGATGTCTTTAAATTCCAGCGACTGTCCGTTGAAATATTAACAATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGc                                                                                                                                   <  2:1306537/150‑1 (MQ=255)
                     ccGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGccc                                                                                                                                 >  2:133610/1‑150 (MQ=255)
                        gtggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                                                                                              <  2:1371954/150‑1 (MQ=255)
                        gtggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                                                                                              >  2:1294584/1‑150 (MQ=255)
                         tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                             <  1:1445047/150‑1 (MQ=255)
                         tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                             <  1:1314785/150‑1 (MQ=255)
                            tAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAt                                                                                                                          <  1:426547/150‑1 (MQ=255)
                            tAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAt                                                                                                                          >  2:1301338/1‑150 (MQ=255)
                            tAATAGTCAGTTCACCCGGGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAt                                                                                                                          <  2:50966/150‑1 (MQ=255)
                                 ttCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                                                                                     <  1:1632539/150‑1 (MQ=255)
                                    aGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                                                                                  >  2:1713208/1‑150 (MQ=255)
                                    aGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATAAAGGGTg                                                                                                                  <  2:1164281/150‑1 (MQ=255)
                                     gTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                                                                                                 >  1:531253/1‑150 (MQ=255)
                                        cACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCAc                                                                                                              <  1:549206/150‑1 (MQ=255)
                                         aCCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGAAcc                                                                                                             >  1:822801/1‑150 (MQ=255)
                                          cccGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCa                                                                                                            >  1:531735/1‑150 (MQ=255)
                                            cGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCAtt                                                                                                          >  1:1509798/1‑150 (MQ=255)
                                             gTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTa                                                                                                         <  1:142014/150‑1 (MQ=255)
                                                  tAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGAc                                                                                                    >  2:515373/1‑150 (MQ=255)
                                                       cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTa                                                                                               >  1:1182041/1‑150 (MQ=255)
                                                       cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATACATAGACCGTTGAATTAAGGGTGGAACCATTACTGACATTTa                                                                                               >  2:50423/1‑150 (MQ=255)
                                                             ccTGTAGCAATTGTGTGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATAc                                                                                         >  2:289553/1‑150 (MQ=255)
                                                             ccTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATAc                                                                                         <  2:422918/150‑1 (MQ=255)
                                                              cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATAcc                                                                                        <  2:1058980/150‑1 (MQ=255)
                                                                   gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTAc                                                                                   <  2:372894/150‑1 (MQ=255)
                                                                   gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTAc                                                                                   >  2:1676439/1‑150 (MQ=255)
                                                                        tgtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGATACCaaaa                                                                              >  2:1068039/1‑150 (MQ=255)
                                                                         gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAt                                                                             <  1:1713208/150‑1 (MQ=255)
                                                                           gCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTa                                                                           >  2:1720223/1‑150 (MQ=255)
                                                                            cGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCAATAGCACGTTGAATTAAGGGTTGCACCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                                          >  1:1434258/1‑150 (MQ=255)
                                                                             gTTTGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCAGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACAATTACGATACCGCTACCAAAATTAcc                                                                         >  1:1430603/1‑150 (MQ=255)
                                                                             gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                                         >  2:376176/1‑150 (MQ=255)
                                                                                agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                                      <  2:1529552/150‑1 (MQ=255)
                                                                                 ttagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                                     <  1:496880/149‑1 (MQ=255)
                                                                                  tagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAt                                                                    >  2:1175185/1‑150 (MQ=255)
                                                                                   agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCATTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                                   >  2:204308/1‑150 (MQ=255)
                                                                                        atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                                              <  2:985200/150‑1 (MQ=255)
                                                                                         tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                                             <  2:531735/150‑1 (MQ=255)
                                                                                           tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGc                                                                                      <  2:551712/123‑1 (MQ=255)
                                                                                           tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGc                                                                                      >  1:551712/1‑123 (MQ=255)
                                                                                            cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                          <  1:1417136/150‑1 (MQ=255)
                                                                                            cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                          >  2:1408475/1‑150 (MQ=255)
                                                                                               tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                                       >  1:1669198/1‑150 (MQ=255)
                                                                                                  aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACAAAAATTACCGACATTGATACCCACAGTATc                                                    >  2:1552086/1‑150 (MQ=255)
                                                                                                  aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACAGTAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAATGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAACGACATTGATACCCACAGAATc                                                    >  1:1058352/1‑150 (MQ=255)
                                                                                                   ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                                                   <  1:1101033/150‑1 (MQ=255)
                                                                                                      cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                                                >  2:212562/1‑150 (MQ=255)
                                                                                                      cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                                                >  2:1405873/1‑150 (MQ=255)
                                                                                                       agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                                               >  2:1406191/1‑150 (MQ=255)
                                                                                                       agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                                               <  1:1720223/150‑1 (MQ=255)
                                                                                                        gagaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                                              <  2:1323740/150‑1 (MQ=255)
                                                                                                         agaCTGTCCGTTGAAATATTAACTATCCCCTAACCGGAGGCATTTCCTCCGATAAATCCTTAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                                             <  2:666657/150‑1 (MQ=255)
                                                                                                          gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGCCATTTACGCTACCGCTACCAAAATTACCGACATTGATAACCACAGTATCGAACgaag                                            >  2:253003/1‑148 (MQ=255)
                                                                                                          gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                                            >  2:25080/1‑150 (MQ=255)
                                                                                                          gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                                            <  2:1509798/150‑1 (MQ=255)
                                                                                                                      aaaTATTAACTATCCCATTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca                                <  1:50423/150‑1 (MQ=255)
                                                                                                                         tatTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa                             <  1:133610/150‑1 (MQ=255)
                                                                                                                         tatTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa                             >  2:1034322/1‑150 (MQ=255)
                                                                                                                          atTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg                            >  1:422179/1‑150 (MQ=255)
                                                                                                                           ttAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg                           <  1:988073/150‑1 (MQ=255)
                                                                                                                              aCTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc                        <  1:1408475/150‑1 (MQ=255)
                                                                                                                               cTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCa                       >  1:1438692/1‑150 (MQ=255)
                                                                                                                                  tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtg                    <  1:25080/150‑1 (MQ=255)
                                                                                                                                  tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtg                    <  1:692423/150‑1 (MQ=255)
                                                                                                                                  tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtg                    <  1:1033346/150‑1 (MQ=255)
                                                                                                                                    cccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCCCCATTACTGACATTTACGATACAGCTACCAAATTTACCGACATTGATACCCAAAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGtgtg                  >  1:299222/1‑150 (MQ=255)
                                                                                                                                      attACCGTAGGAATTTCCTCAGATAAATCCATAGAACGTTGAATTAAGTGTGGCACAATTACTTACATTTTCGATACAGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGATCAGTGTGcc                <  1:80290/149‑1 (MQ=255)
                                                                                                                                       ttACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCg               <  1:953199/150‑1 (MQ=255)
                                                                                                                                        tatcGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAATTTACCGACAATGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGa              >  1:290466/4‑150 (MQ=255)
                                                                                                                                        tACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGa              >  1:593621/1‑150 (MQ=255)
                                                                                                                                         atcGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGAt             >  1:1010610/3‑150 (MQ=255)
                                                                                                                                          ccGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATa            >  2:1608247/1‑150 (MQ=255)
                                                                                                                                          ccGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATa            <  2:1014878/150‑1 (MQ=255)
                                                                                                                                            ggaggaAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc          >  1:1547045/1‑150 (MQ=255)
                                                                                                                                              aGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTTCCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgc        >  1:454101/1‑150 (MQ=255)
                                                                                                                                               ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgct       >  2:1023243/1‑150 (MQ=255)
                                                                                                                                                 cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGTTTCAGTGTGCCTATACCgcttg     >  1:151686/1‑150 (MQ=255)
                                                                                                                                                   tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgct   <  1:515373/150‑1 (MQ=255)
                                                                                                                                                    ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgctt  <  1:1301338/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTT  >  NC_000913/1596354‑1596651

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