Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 596,072 A→C 6.9% Q158P (CAG→CCG)  cusC → copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913596,0720AC 6.9% 1523.4 / 17.4 535Q158P (CAG→CCG) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base A (219/262);  new base C (13/37);  total (235/300)
Fisher's exact test for biased strand distribution p-value = 1.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 542 total aligned reads displayed.

GCTGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGTTGACAGGTA  >  NC_000913/595925‑596221
                                                                                                                                                   |                                                                                                                                                     
 cTGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAgcg                                                                                                                                                    >  1:2903834/1‑150 (MQ=255)
  tGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAgcgc                                                                                                                                                   <  2:1484079/150‑1 (MQ=255)
  tGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAgcgc                                                                                                                                                   >  2:7313165/1‑150 (MQ=255)
   ggAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAgcgcg                                                                                                                                                  >  1:5889247/1‑150 (MQ=255)
     aGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGcgcgcg                                                                                                                                                >  1:6101878/1‑150 (MQ=255)
      gcgggcaTCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGACGCTCCGCGCGCgg                                                                                                                                               >  2:1754318/7‑150 (MQ=255)
      gCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTTACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCgg                                                                                                                                               <  2:3864365/150‑1 (MQ=255)
      gCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCgg                                                                                                                                               >  2:2185939/1‑150 (MQ=255)
       cGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGt                                                                                                                                              >  2:2244506/1‑150 (MQ=255)
        ggTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTg                                                                                                                                             <  2:4828946/150‑1 (MQ=255)
             tCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCatat                                                                                                                                        >  2:3750275/1‑150 (MQ=255)
              cTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATAtt                                                                                                                                       >  2:6058724/1‑150 (MQ=255)
               ttAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAATATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTc                                                                                                                                      >  2:4177034/1‑150 (MQ=255)
                  aaGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTctgc                                                                                                                                   <  1:2580457/150‑1 (MQ=255)
                  aaGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTctgc                                                                                                                                   >  1:7776517/1‑150 (MQ=255)
                  aaGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTctgc                                                                                                                                   <  1:3678474/150‑1 (MQ=255)
                   aGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCATCGCGCGGTGCATATTctgat                                                                                                                                  >  2:4397612/1‑148 (MQ=255)
                    gggAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTctgctg                                                                                                                                 >  2:7472832/1‑150 (MQ=255)
                      gAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGt                                                                                                                               >  1:3146244/1‑150 (MQ=255)
                       aaTACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGACACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGtt                                                                                                                              >  1:7021070/1‑150 (MQ=255)
                           ccGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTa                                                                                                                          <  1:1912037/150‑1 (MQ=255)
                            cGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTaa                                                                                                                         >  1:1103126/1‑150 (MQ=255)
                              ccACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATg                                                                                                                       >  1:5725676/1‑150 (MQ=255)
                              ccACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATg                                                                                                                       <  1:6647984/150‑1 (MQ=255)
                                   acacGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTcgcg                                                                                                                  <  2:3263688/150‑1 (MQ=255)
                                       gggAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCaaa                                                                                                              <  2:5894534/150‑1 (MQ=255)
                                       gggAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCaaa                                                                                                              <  2:7968325/150‑1 (MQ=255)
                                       gggAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCaaa                                                                                                              <  1:1004157/150‑1 (MQ=255)
                                             tttCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTAtt                                                                                                        <  1:3271670/150‑1 (MQ=255)
                                             tttCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTAtt                                                                                                        <  2:4733737/150‑1 (MQ=255)
                                             tttCGACCGGCCTTAACGCCAGCATTGACCTCGATATTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTAtt                                                                                                        <  2:6601628/150‑1 (MQ=255)
                                                  aCCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAAt                                                                                                   <  1:3486071/150‑1 (MQ=255)
                                                  aCCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAAt                                                                                                   <  1:1203199/150‑1 (MQ=255)
                                                           aaCGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCATCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTg                                                                                          <  2:5305088/150‑1 (MQ=255)
                                                             cGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGc                                                                                        <  2:3461759/150‑1 (MQ=255)
                                                             cGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGc                                                                                        <  1:6241290/150‑1 (MQ=255)
                                                               ccAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACAAGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGt                                                                                      >  1:7673296/1‑150 (MQ=255)
                                                                cAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTa                                                                                     <  1:7472832/150‑1 (MQ=255)
                                                                    tttGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCCGCAACTGGCGTATgcg                                                                                 >  1:441792/1‑150 (MQ=255)
                                                                    tttGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATgcg                                                                                 <  2:2715960/150‑1 (MQ=255)
                                                                        aCCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAAt                                                                             >  2:6259043/1‑150 (MQ=255)
                                                                         ccTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAAtt                                                                            <  2:6735915/150‑1 (MQ=255)
                                                                          cTCGATTTTTTCGGTCGCTTTAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGACGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTg                                                                           >  1:7576920/1‑150 (MQ=255)
                                                                          cTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTg                                                                           >  2:2374841/1‑150 (MQ=255)
                                                                              aTTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCaaa                                                                       <  1:7836161/150‑1 (MQ=255)
                                                                              aTTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCaaa                                                                       <  1:6479588/150‑1 (MQ=255)
                                                                               ttttttCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAAt                                                                      >  1:3035378/1‑150 (MQ=255)
                                                                                tttttCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATa                                                                     <  2:2567235/150‑1 (MQ=255)
                                                                                 ttttCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAg                                                                    <  2:290929/150‑1 (MQ=255)
                                                                                 ttttCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAg                                                                    <  2:2028824/150‑1 (MQ=255)
                                                                                 ttttCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAg                                                                    <  1:2358352/150‑1 (MQ=255)
                                                                                  tttCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGc                                                                   <  1:4374336/150‑1 (MQ=255)
                                                                                      ggTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCgaa                                                               <  2:4245262/150‑1 (MQ=255)
                                                                                        tCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTCTGCGCACTTGCAAATAGCCgaaga                                                             >  2:6385651/1‑150 (MQ=255)
                                                                                          gCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATAGCAAATAGCCGAAGaaa                                                           >  2:6106959/1‑150 (MQ=255)
                                                                                           cTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCTACTGGCGTATGCGCAATTGCAAATAGCCGATGAAAc                                                          >  2:8141615/1‑150 (MQ=255)
                                                                                              aaaGAACATGAGCGAAGCCGAGCGACAATATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCt                                                       <  2:51030/150‑1 (MQ=255)
                                                                                              aaaGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCt                                                       <  2:1014080/150‑1 (MQ=255)
                                                                                               aaGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTg                                                      <  2:7235529/150‑1 (MQ=255)
                                                                                               aaGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTg                                                      <  2:90825/150‑1 (MQ=255)
                                                                                                 gAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCg                                                    >  2:1545475/1‑150 (MQ=255)
                                                                                                    cATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTaa                                                 <  2:5637970/150‑1 (MQ=255)
                                                                                                     aTGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAAt                                                <  2:1642247/150‑1 (MQ=255)
                                                                                                        aGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTAt                                             >  2:3657054/1‑150 (MQ=255)
                                                                                                        aGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTATGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTAt                                             <  1:5416822/150‑1 (MQ=255)
                                                                                                           gAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATcag                                          >  2:2444171/1‑150 (MQ=255)
                                                                                                            aaGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATcagc                                         <  2:6850846/150‑1 (MQ=255)
                                                                                                            aaGCCGAGCGACAAAATTATTTAGCCAATGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATcagc                                         <  1:1575421/150‑1 (MQ=255)
                                                                                                             aGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATcagca                                        <  1:6835528/150‑1 (MQ=255)
                                                                                                             aGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATcagca                                        <  2:5701045/150‑1 (MQ=255)
                                                                                                             aGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATcagca                                        <  1:693452/150‑1 (MQ=255)
                                                                                                               ccGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGt                                      <  1:2261697/150‑1 (MQ=255)
                                                                                                                cGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTc                                     <  2:4281807/150‑1 (MQ=255)
                                                                                                                cGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTc                                     <  2:7094336/150‑1 (MQ=255)
                                                                                                                    cGACAAAATTATTTAGCCTCTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTCATTATCCGCAGTCatat                                 <  2:5383885/150‑1 (MQ=255)
                                                                                                                    cGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCatat                                 <  1:7369579/150‑1 (MQ=255)
                                                                                                                       cAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCg                              <  1:3927258/150‑1 (MQ=255)
                                                                                                                       cAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCg                              <  1:1400349/150‑1 (MQ=255)
                                                                                                                        aaaaTTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGt                             <  1:8119155/150‑1 (MQ=255)
                                                                                                                        aaaaTTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGt                             <  1:3897881/150‑1 (MQ=255)
                                                                                                                        aaaaTTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGt                             <  2:7596073/150‑1 (MQ=255)
                                                                                                                         aaaTTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGtt                            >  1:5260698/1‑150 (MQ=255)
                                                                                                                            ttattTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGt                         <  1:7195250/150‑1 (MQ=255)
                                                                                                                             tattTAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTc                        <  1:4772359/150‑1 (MQ=255)
                                                                                                                             tattTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTc                        <  1:1471246/150‑1 (MQ=255)
                                                                                                                               tttAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGa                      >  1:5144861/1‑150 (MQ=255)
                                                                                                                               tttAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGa                      <  1:412272/150‑1 (MQ=255)
                                                                                                                               tttAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGa                      <  2:7737684/150‑1 (MQ=255)
                                                                                                                                ttAGCCACTGAGGAAGCTCCGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGaa                     <  1:3069306/150‑1 (MQ=255)
                                                                                                                                   gCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGaaaaa                  <  2:777939/150‑1 (MQ=255)
                                                                                                                                    ccACTGAGGAAGCTCGGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAaac                 >  1:7172654/1‑150 (MQ=255)
                                                                                                                                    ccACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAaac                 <  2:4464316/150‑1 (MQ=255)
                                                                                                                                     cACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAaaca                >  2:4792497/1‑150 (MQ=255)
                                                                                                                                       cTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATTTGCGTTTGTCGAAAaacaac              <  2:6941988/150‑1 (MQ=255)
                                                                                                                                        tGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACt             >  1:465707/1‑150 (MQ=255)
                                                                                                                                          aGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGt           >  1:7972095/1‑150 (MQ=255)
                                                                                                                                          aGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGt           >  2:3089787/1‑150 (MQ=255)
                                                                                                                                             aaGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGTTGa        <  2:122730/150‑1 (MQ=255)
                                                                                                                                              aGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGTTGAc       <  2:2323010/150‑1 (MQ=255)
                                                                                                                                              aGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGTTGAc       >  2:4611672/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTGGAGCGGTAATCTTAAAGGGAATACCGCCACGACACGGGAGTTTTCGACCGGCCTTAACGCCAGCTTTGACCTCGATTTTTTCGGTCGCTTAAAGAACATGAGCGAAGCCGAGCGACAAAATTATTTAGCCACTGAGGAAGCTCAGCGCGCGGTGCATATTCTGCTGGTTTCTAATGTCGCGCAAAGCTATTTCAATCAGCAACTGGCGTATGCGCAATTGCAAATAGCCGAAGAAACGCTGCGTAATTATCAGCAGTCATATGCGTTTGTCGAAAAACAACTGTTGACAGGTA  >  NC_000913/595925‑596221

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