Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,596,499 C→A 5.9% pseudogene (1489/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,4990CA 5.9% 317.9 / 5.9 105pseudogene (1489/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base C (50/48);  new base A (4/3);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

GTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTG  >  NC_000913/1596351‑1596648
                                                                                                                                                    |                                                                                                                                                     
gtATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCa                                                                                                                                                      >  1:1384212/1‑150 (MQ=255)
gtATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGAACTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAATTATTAACTATCCACTTACCGGAGGCa                                                                                                                                                      >  1:546619/1‑150 (MQ=255)
 tATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAt                                                                                                                                                     >  2:118500/1‑150 (MQ=255)
     aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGATATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                                 >  2:1051122/1‑150 (MQ=255)
     aCGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcc                                                                                                                                                 >  2:1773124/1‑150 (MQ=255)
      cgcgCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcct                                                                                                                                                >  1:881830/1‑150 (MQ=255)
        cgcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                              >  1:1127916/1‑150 (MQ=255)
        cgcTTTAACAATACCTCCGGTGGTAATATTCAGGTCACCCGTCCCTAATACACCGACCTGTAGCAATTGTGAGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtcctcc                                                                                                                                              <  2:297322/150‑1 (MQ=255)
         gcTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCg                                                                                                                                             <  2:1403732/150‑1 (MQ=255)
          gtttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCgg                                                                                                                                            <  2:164149/149‑2 (MQ=255)
              aaCTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAGCTATCCCCTTACCGGAGGCATTTCCTCCGATaaa                                                                                                                                        >  1:52749/1‑150 (MQ=255)
               aCTATACCTCCGGTGGTAATATTCAGATCACCCGTGCCTAATACACCGACCTGTAGCAATTCTGCGTTAGTAGATGATGTCTTTAAAATCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTaagg                                                                                                                                       <  1:104082/150‑3 (MQ=255)
                cTATACCTCCGGTGGTAATATTCTGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGTAGGCATTTCCTCCGATAAATc                                                                                                                                      <  1:920037/150‑1 (MQ=255)
                 tataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATcc                                                                                                                                     >  2:1861607/1‑150 (MQ=255)
                  ataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCa                                                                                                                                    >  2:1605777/1‑150 (MQ=255)
                       tCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGcc                                                                                                                               <  1:1554569/150‑1 (MQ=255)
                         cGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCg                                                                                                                             <  2:418344/150‑1 (MQ=255)
                         cGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCg                                                                                                                             <  2:461529/150‑1 (MQ=255)
                         cGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCg                                                                                                                             <  2:803733/150‑1 (MQ=255)
                           gtggtAATATTCAGTTCACCCGTGCCTAATCCCCCGACCGGTAGCAATTGTGCGGTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCAGATAAATCCATAGCCCGtt                                                                                                                           >  2:941868/1‑150 (MQ=255)
                            tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                          <  1:1080288/150‑1 (MQ=255)
                              gtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                                                                                        <  1:1103983/150‑1 (MQ=255)
                                aaTATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAtt                                                                                                                      >  2:236010/1‑150 (MQ=255)
                                aaTATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAtt                                                                                                                      <  1:104829/150‑1 (MQ=255)
                                    ttCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTTTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGAGTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                                                                                  <  1:607548/150‑1 (MQ=255)
                                    ttCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGCCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                                                                                  >  1:895785/1‑150 (MQ=255)
                                      cAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                                                                                >  2:108837/1‑150 (MQ=255)
                                       aGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATTATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                                                                               <  1:1486280/150‑1 (MQ=255)
                                        gTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                                                                                              >  2:1773003/1‑150 (MQ=255)
                                          tCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCa                                                                                                            >  1:825715/1‑150 (MQ=255)
                                          tCAACCGTGCCTAATACACCGAACTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCa                                                                                                            <  2:1089468/150‑1 (MQ=255)
                                               cGTGCCTTATACACCGACCTTTAGCAATTGTACGTTATTAGATTATGTCTTTAAATTCCATATACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCAtt                                                                                                       <  2:140545/150‑1 (MQ=255)
                                               cGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCAtt                                                                                                       <  1:1596646/150‑1 (MQ=255)
                                                 tGCCTAATACACCGACCTGTAGCAATTGGGCGTTAGTAGAAGATGTCTTTAAAGTCCAGAGACTGTCCGTAGAAATATTAACTAACCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTAc                                                                                                     <  2:1848769/150‑1 (MQ=255)
                                                  gCCTAATACACCGACCTGTAGCAATTCTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGATGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACt                                                                                                    <  2:1168598/150‑1 (MQ=255)
                                                   ccTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTg                                                                                                   <  1:1557723/150‑1 (MQ=255)
                                                    cTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGa                                                                                                  >  2:1269049/1‑150 (MQ=255)
                                                    cTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGa                                                                                                  <  1:439398/150‑1 (MQ=255)
                                                     tAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGAc                                                                                                 <  2:1374037/150‑1 (MQ=255)
                                                      aaTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACa                                                                                                <  1:1684163/150‑1 (MQ=255)
                                                      aaTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACa                                                                                                >  2:824230/1‑150 (MQ=255)
                                                       aTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCGTTAAATTCCAGAGACTGTCCGTTGAAATATTACCTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTGGAATTATGGGTGGCCCAATTACTGACAt                                                                                               >  2:1146228/1‑150 (MQ=255)
                                                          cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTACGGGTGGCACCATTACTGACATTTa                                                                                            >  1:413873/1‑150 (MQ=255)
                                                          cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTa                                                                                            >  1:561242/1‑150 (MQ=255)
                                                          cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATAAACTATCCCCTTACCGGAGGAATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTa                                                                                            <  2:1206198/150‑1 (MQ=255)
                                                            ccGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACg                                                                                          <  1:546603/150‑1 (MQ=255)
                                                               aCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCACCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATa                                                                                       >  2:1484564/1‑150 (MQ=255)
                                                                 cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATAcc                                                                                     >  1:592735/1‑150 (MQ=255)
                                                                 cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTAACCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGGTGAATTCAGGGTGGCACCATTACTGACATTTACGATAcc                                                                                     <  2:1127916/150‑1 (MQ=255)
                                                                  tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCAATTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCg                                                                                    >  1:777814/1‑150 (MQ=255)
                                                                     aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTa                                                                                 <  1:382555/150‑1 (MQ=255)
                                                                     aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTa                                                                                 >  2:1183258/1‑150 (MQ=255)
                                                                         aatgtgCTTTATTATATGATGTCTTTATATTCCAGAGACTGTCCTTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCAAAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCaa                                                                             <  1:1658468/148‑1 (MQ=255)
                                                                             tgCGTTAGTAGATGATGTCTTTAACTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAtt                                                                         >  1:1323701/1‑150 (MQ=255)
                                                                                gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGATGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                                      <  1:280932/150‑1 (MQ=255)
                                                                                 ttAGTTGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCg                                                                     <  1:1461735/150‑1 (MQ=255)
                                                                                 ttAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCg                                                                     <  1:636187/150‑1 (MQ=255)
                                                                                   agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                                   >  1:1830161/1‑150 (MQ=255)
                                                                                      agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                                <  1:118500/150‑1 (MQ=255)
                                                                                      agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                                >  1:1254979/1‑150 (MQ=255)
                                                                                      agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                                <  1:1783984/150‑1 (MQ=255)
                                                                                           atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATGGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAc                                                           >  1:1372296/1‑150 (MQ=255)
                                                                                               cTTTATATTCCAGAGACTGTCCGTTGAAATATTTACTATCCCCTTTCCGGAGGCATTTCCTCCGATAAATCAATAGCCCGTTTAATTAAGTGTGGCACCATTACTGACATTTACTATACCGCTACCAAAATTACCGTCATTGATACCcac                                                       <  1:1125824/150‑1 (MQ=255)
                                                                                               cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAATGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTTCCGACATTGATACCcac                                                       <  1:236010/150‑1 (MQ=255)
                                                                                               cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATacacac                                                       >  2:1801392/1‑150 (MQ=255)
                                                                                                tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcaca                                                      >  1:679962/1‑150 (MQ=255)
                                                                                                   aaaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGATGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTa                                                   <  2:546619/150‑1 (MQ=255)
                                                                                                   aaaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTa                                                   >  1:269003/1‑150 (MQ=255)
                                                                                                    aaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                                  >  1:404431/1‑150 (MQ=255)
                                                                                                       tCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGa                                               >  1:852573/1‑150 (MQ=255)
                                                                                                        ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                                              >  1:436884/1‑150 (MQ=255)
                                                                                                         cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGCCATTGATACCCAGAGTATCGAAc                                             >  1:473793/1‑150 (MQ=255)
                                                                                                         cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAc                                             >  1:1591950/1‑150 (MQ=255)
                                                                                                          agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                                            >  1:200437/1‑150 (MQ=255)
                                                                                                            agaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgat                                          <  1:35605/150‑1 (MQ=255)
                                                                                                                tGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGatgatg                                      <  1:1178741/150‑1 (MQ=255)
                                                                                                                 gTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAAAGATGATGt                                     >  2:90095/1‑150 (MQ=255)
                                                                                                                    cGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                                  <  1:1795234/150‑1 (MQ=255)
                                                                                                                       tGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatc                               <  1:1773124/150‑1 (MQ=255)
                                                                                                                        gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc                              >  2:1837082/1‑150 (MQ=255)
                                                                                                                        gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc                              <  2:1525889/150‑1 (MQ=255)
                                                                                                                           atatTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccatc                           >  1:474100/1‑150 (MQ=255)
                                                                                                                             atTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCCCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg                         >  1:1524141/1‑150 (MQ=255)
                                                                                                                             atTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCAGAGTATCGAACGATGATGTCGCTTCCATCAg                         >  1:1224397/1‑150 (MQ=255)
                                                                                                                                 aCTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTc                     <  1:1183258/150‑1 (MQ=255)
                                                                                                                                    atggatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgt                  >  2:406978/7‑150 (MQ=255)
                                                                                                                                     tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCTCCATTACTGAGATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCTTCAGGTTCAgtg                 >  1:739468/1‑150 (MQ=255)
                                                                                                                                     tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtg                 <  1:1598844/150‑1 (MQ=255)
                                                                                                                                     tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtg                 <  2:1254979/150‑1 (MQ=255)
                                                                                                                                           tatcGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGACGa           >  1:601539/4‑150 (MQ=255)
                                                                                                                                           tatcGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCGTCCAAAATTACCGACATTGATACCCACAGTAGCGAACGAGGATGTCGCATCCATCAGGGTACGTGTGCCGa           >  1:1336483/4‑150 (MQ=255)
                                                                                                                                           tACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGGTCAGTGTGCCGa           <  2:1384212/150‑1 (MQ=255)
                                                                                                                                              cGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAc        >  1:945157/1‑150 (MQ=255)
                                                                                                                                               ggAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc       >  2:803273/1‑150 (MQ=255)
                                                                                                                                                gAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCg      <  2:206064/150‑1 (MQ=255)
                                                                                                                                                  ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgct    <  2:399339/150‑1 (MQ=255)
                                                                                                                                                  ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgct    <  1:1126997/150‑1 (MQ=255)
                                                                                                                                                    cATTTCCTCCGATAAATCCATAGCCCGTTGAATTACGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttg  <  2:413873/150‑1 (MQ=255)
                                                                                                                                                    cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttg  >  1:124977/1‑150 (MQ=255)
                                                                                                                                                    cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttg  >  2:1755730/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTG  >  NC_000913/1596351‑1596648

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: