Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,576,847 G→A 5.2% pseudogene (2713/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,8470GA 5.2% 306.5 / 4.0 95pseudogene (2713/2721 nt)yjiVputative uncharacterized protein YjiV
Reads supporting (aligned to +/- strand):  ref base G (43/47);  new base A (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 = 4.26e-01

CGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGG  >  NC_000913/4576700‑4576985
                                                                                                                                                   |                                                                                                                                          
cgcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcc                                                                                                                                          >  1:662629/1‑150 (MQ=255)
 gcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcct                                                                                                                                         <  1:412022/150‑1 (MQ=255)
 gcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcct                                                                                                                                         <  2:1353779/150‑1 (MQ=255)
 gcATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcct                                                                                                                                         >  1:1631012/1‑150 (MQ=255)
  cATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctc                                                                                                                                        >  2:892327/1‑150 (MQ=255)
  cATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctc                                                                                                                                        <  2:1662951/150‑1 (MQ=255)
  cATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGcctc                                                                                                                                        >  2:1824044/1‑150 (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:427160/150‑1 (MQ=255)
       aCTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgat                                                                                                                                   <  1:419891/150‑1 (MQ=255)
       aCTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgat                                                                                                                                   >  1:1579558/1‑150 (MQ=255)
        cTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGCGATCAcctccaccag                                                                                                                                  <  2:136841/150‑6 (MQ=255)
         ttttCATCGCAATGGCATTACTTTGATTTCCGCTATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgatga                                                                                                                                 <  1:704589/150‑1 (MQ=255)
         ttttCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgatga                                                                                                                                 <  2:916250/150‑1 (MQ=255)
         ttttCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCtgatga                                                                                                                                 <  2:899917/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)
                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                                                                                                              <  1:191676/150‑1 (MQ=255)
                            aCTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATgtgg                                                                                                              <  2:662629/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                                                                                           >  1:551193/1‑150 (MQ=255)
                                               cGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTc                                                                                           <  2:285663/150‑1 (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:550039/150‑1 (MQ=255)
                                                              aGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGc                                                                            >  1:236207/1‑150 (MQ=255)
                                                              aGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGc                                                                            <  1:1337078/150‑1 (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                                                                                                                           <  2:651902/93‑1 (MQ=255)
                                                                        gCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGcc                                                                                                                           >  1:651902/1‑93 (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:366849/1‑150 (MQ=255)
                                                                                    aaggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa                                                      >  1:1245717/1‑150 (MQ=255)
                                                                                     aggaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaaa                                                     >  1:667605/1‑150 (MQ=255)
                                                                                       gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATg                                                   >  1:399370/1‑150 (MQ=255)
                                                                                       gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATg                                                   >  1:1697878/1‑150 (MQ=255)
                                                                                         aggTAATGTCGCCAACAGCGAAGAGAGTTGGGGAAGGGATGTGGTGGTTGAGGTGATCGCCTCCTGATTTTGAGCCGCTCCGGATGTGGTGTCGGGATCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                 <  1:1297921/150‑1 (MQ=255)
                                                                                           ggAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa                                               >  2:612564/1‑150 (MQ=255)
                                                                                           ggAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa                                               <  1:1626966/150‑1 (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)
                                                                                                                                                   |                                                                                                                                          
CGCATCAACTTTTCATCGCAATGGCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGG  >  NC_000913/4576700‑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: