Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,576,842 T→C 5.2% pseudogene (2708/2721 nt) yjiV → putative uncharacterized protein YjiV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,576,8420TC 5.2% 292.7 / 4.8 95pseudogene (2708/2721 nt)yjiVputative uncharacterized protein YjiV
Reads supporting (aligned to +/- strand):  ref base T (43/47);  new base C (2/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGG  >  NC_000913/4576693‑4576985
                                                                                                                                                     |                                                                                                                                               
cTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGt                                                                                                                                                 >  1:1861013/1‑150 (MQ=255)
       cgcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcc                                                                                                                                          >  1:662629/1‑150 (MQ=255)
        gcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcct                                                                                                                                         <  1:412022/150‑1 (MQ=255)
        gcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcct                                                                                                                                         >  1:1631012/1‑150 (MQ=255)
        gcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcct                                                                                                                                         <  2:1353779/150‑1 (MQ=255)
         cATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctc                                                                                                                                        >  2:892327/1‑150 (MQ=255)
         cATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctc                                                                                                                                        >  2:1824044/1‑150 (MQ=255)
         cATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctc                                                                                                                                        <  2:1662951/150‑1 (MQ=255)
          aTCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctcc                                                                                                                                       >  2:1121008/1‑150 (MQ=255)
           tCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctcct                                                                                                                                      <  2:1107671/150‑1 (MQ=255)
            cAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtg                                                                                                                                     <  2:1861013/150‑1 (MQ=255)
             aaCTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGCGATCAcctccacc                                                                                                                                    >  1:136841/1‑147 (MQ=255)
              aCTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgat                                                                                                                                   <  1:419891/150‑1 (MQ=255)
              aCTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgat                                                                                                                                   >  1:1579558/1‑150 (MQ=255)
              aCTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgat                                                                                                                                   <  1:427160/150‑1 (MQ=255)
               cTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGCGATCAcctccaccag                                                                                                                                  <  2:136841/150‑6 (MQ=255)
                ttttCATCGCAATGGCATTACTTTGATTTCCGCTATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgatga                                                                                                                                 <  1:704589/150‑1 (MQ=255)
                ttttCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgatga                                                                                                                                 <  2:899917/150‑1 (MQ=255)
                ttttCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgatga                                                                                                                                 <  2:916250/150‑1 (MQ=255)
                ttttCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGCGATCAcctccaccagc                                                                                                                                 <  1:1614522/150‑7 (MQ=255)
                  ttCATCGCAATGGCATTACTTTGATTTCCGCTATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGTTGTGCTGGTGGAGGTGATCGCCTCCtgatgatg                                                                                                                               <  1:1562115/150‑1 (MQ=255)
                     aTCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTGAAAATGGCTCAGGCCTGCATGGAAGGGAATGTCGCCAACAGCGATGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGc                                                                                                                            <  1:1793063/150‑1 (MQ=255)
                     aTCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGt                                                                                                                                                 >  1:346332/1‑129 (MQ=255)
                     aTCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGt                                                                                                                                                 <  2:346332/129‑1 (MQ=255)
                       cGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGTCg                                                                                                                          >  2:1041111/1‑150 (MQ=255)
                        gCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGc                                                                                                                         >  2:1907694/1‑150 (MQ=255)
                          aaTGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTc                                                                                                                       <  1:1660097/150‑1 (MQ=255)
                            tGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTccc                                                                                                                     <  2:329866/150‑1 (MQ=255)
                               cATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCGTGATGATGAGCCGCTCCCGAt                                                                                                                  >  2:995475/1‑150 (MQ=255)
                               cATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCGTGATGATGAGCCGCTCCCGAt                                                                                                                  >  2:1002180/1‑150 (MQ=255)
                               cATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGAt                                                                                                                  >  1:1246254/1‑150 (MQ=255)
                                  tACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGAtgtg                                                                                                               <  2:956180/150‑1 (MQ=255)
                                   aCTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATgtgg                                                                                                              <  2:662629/150‑1 (MQ=255)
                                   aCTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATgtgg                                                                                                              <  1:191676/150‑1 (MQ=255)
                                   aCTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTCATCGCCTCCTGATGATGAGCCGCTCCCGATgtgg                                                                                                              >  2:435013/1‑150 (MQ=255)
                                      atGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACTGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGgtgt                                                                                                           >  2:1297921/2‑150 (MQ=255)
                                       tGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTc                                                                                                          <  1:1313857/150‑1 (MQ=255)
                                        gATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCg                                                                                                         <  1:631749/150‑1 (MQ=255)
                                          tttCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCggg                                                                                                       <  1:127384/150‑1 (MQ=255)
                                             ccGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGc                                                                                                    <  1:1907694/150‑1 (MQ=255)
                                               gcgATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCgg                                                                                                  <  2:564567/150‑1 (MQ=255)
                                                cgATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGt                                                                                                 >  1:1978188/1‑150 (MQ=255)
                                                    gACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTAttt                                                                                             >  2:1535364/1‑150 (MQ=255)
                                                     aCGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTAtttt                                                                                            >  1:1622713/1‑150 (MQ=255)
                                                      cGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTc                                                                                           <  2:285663/150‑1 (MQ=255)
                                                      cGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTc                                                                                           >  1:551193/1‑150 (MQ=255)
                                                        tttctttcCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCta                                                                                         <  1:558299/150‑1 (MQ=255)
                                                           ctttcCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATaa                                                                                      >  1:1364954/1‑150 (MQ=255)
                                                              tcCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAAc                                                                                   <  2:1851108/150‑1 (MQ=255)
                                                               ccAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACt                                                                                  >  2:701548/1‑150 (MQ=255)
                                                                 aGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTa                                                                                >  2:416339/1‑150 (MQ=255)
                                                                     aGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGc                                                                            <  1:1337078/150‑1 (MQ=255)
                                                                     aGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGc                                                                            <  1:550039/150‑1 (MQ=255)
                                                                     aGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGc                                                                            >  1:236207/1‑150 (MQ=255)
                                                                     aGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGc                                                                            >  1:1170381/1‑150 (MQ=255)
                                                                      gTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCt                                                                           >  2:1626966/1‑150 (MQ=255)
                                                                          aaaTGGCTCAGGTCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTAtt                                                                       >  1:984902/1‑150 (MQ=255)
                                                                              ggCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgag                                                                   <  2:236207/150‑1 (MQ=255)
                                                                              ggCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgag                                                                   <  1:1899425/150‑1 (MQ=255)
                                                                               gCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGcc                                                                                                                           >  1:651902/1‑93 (MQ=255)
                                                                               gCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGcc                                                                                                                           <  2:651902/93‑1 (MQ=255)
                                                                               gCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgaga                                                                  >  1:1307622/1‑150 (MQ=255)
                                                                                cTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGat                                                                 <  2:1246254/150‑1 (MQ=255)
                                                                                     gCCTGCAAGGAAGGGAATGTCTCCAACAGCGAAGAGAGTTGGCCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGATCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCCCTTATTTGAGATATTCa                                                            >  2:117572/1‑150 (MQ=255)
                                                                                     gCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCa                                                            >  1:1147321/1‑150 (MQ=255)
                                                                                           aaggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa                                                      >  1:1245717/1‑150 (MQ=255)
                                                                                           aaggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa                                                      >  1:366849/1‑150 (MQ=255)
                                                                                            aggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaaa                                                     >  1:667605/1‑150 (MQ=255)
                                                                                              gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATg                                                   >  1:1697878/1‑150 (MQ=255)
                                                                                              gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATg                                                   >  1:399370/1‑150 (MQ=255)
                                                                                                aggTAATGTCGCCAACAGCGAAGAGAGTTGGGGAAGGGATGTGGTGGTTGAGGTGATCGCCTCCTGATTTTGAGCCGCTCCGGATGTGGTGTCGGGATCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                 <  1:1297921/150‑1 (MQ=255)
                                                                                                  ggAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa                                               <  1:1626966/150‑1 (MQ=255)
                                                                                                  ggAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa                                               >  2:612564/1‑150 (MQ=255)
                                                                                                    aaTGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCTTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt                                             <  1:995475/150‑1 (MQ=255)
                                                                                                    aaTGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt                                             <  2:366849/150‑1 (MQ=255)
                                                                                                    aaTGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGt                                             <  1:1002180/150‑1 (MQ=255)
                                                                                                          gCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGAGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCt                                       >  1:1527465/1‑150 (MQ=255)
                                                                                                            cAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGAGGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTg                                     <  1:1824044/150‑1 (MQ=255)
                                                                                                              aCAGCGAAGAGAGTTGGGCAACGGATGTGCTGGGGGTGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCGTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGAt                                   <  1:1032484/150‑1 (MQ=255)
                                                                                                               cAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATg                                  >  1:467162/1‑150 (MQ=255)
                                                                                                                 gCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGta                                <  1:701548/150‑1 (MQ=255)
                                                                                                                  cGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtat                               <  2:1170381/150‑1 (MQ=255)
                                                                                                                     agagagTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGt                                 >  1:875401/1‑145 (MQ=255)
                                                                                                                     agagagTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGt                                 <  2:875401/145‑1 (MQ=255)
                                                                                                                             gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCAtt                    >  2:223552/1‑150 (MQ=255)
                                                                                                                              ggCAACGGAAGTGCTGGTGCAGGTGATCGCCTCCTGATGATGATCCGCTCCCGATGTCGAGTCCGGAGCGGTATTATCTATAAAACTAACCGCTTATTTGAGATAATCATCGAAAATGTCCAGTAATTCTTGATGTATACACGCCCAtta                   >  2:351151/1‑149 (MQ=255)
                                                                                                                               gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGTCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc                  <  1:1041111/150‑1 (MQ=255)
                                                                                                                               gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc                  <  1:416339/150‑1 (MQ=255)
                                                                                                                               gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCACCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc                  >  2:1482782/1‑150 (MQ=255)
                                                                                                                                        gtgATGCTGGAGGTGATCGCCTAATCATGATGACCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTaaa         >  2:1891385/1‑150 (MQ=255)
                                                                                                                                           cTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTg      <  1:1535364/150‑1 (MQ=255)
                                                                                                                                            atcaggAGGCGATCACCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc                  <  2:442795/133‑1 (MQ=255)
                                                                                                                                            atcaggAGGCGATCACCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc                  >  1:442795/5‑137 (MQ=255)
                                                                                                                                             ggtggAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGAc    <  2:1147321/150‑1 (MQ=255)
                                                                                                                                             ggtggAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGAc    >  2:380956/1‑150 (MQ=255)
                                                                                                                                               tggAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata                              <  2:1852478/122‑1 (MQ=255)
                                                                                                                                               tggAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGtata                              >  1:1852478/1‑122 (MQ=255)
                                                                                                                                               tggAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACgg  <  2:1978188/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CTACTGGCGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGG  >  NC_000913/4576693‑4576985

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