Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,559,351 A→C 6.4% I59S (ATC→AGC) ‡ iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,559,3510AC 6.4% 461.9 / 15.1 165I59S (ATC→AGC) ‡iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (83/71);  new base C (5/6);  total (88/77)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 166 total aligned reads displayed.

TAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGT  >  NC_000913/4559202‑4559499
                                                                                                                                                     |                                                                                                                                                    
tAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGa                                                                                                                                                      <  2:412064/150‑1 (MQ=255)
 agagTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAATTGGACGTGTTGATCAATAAAACCTGGGCAGAGGAt                                                                                                                                                     >  2:410845/1‑150 (MQ=255)
   agTCGGTGCCCAGCAGACCAACCACTGACGTGCCGCCCGCTTCCGTCAGGCGACTTAGCGCCCCTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCCCCGCCCATCAAATGGCCGTGTTGATCAATAAAACCTGGGCAGCGGATcg                                                                                                                                                   >  1:191283/1‑149 (MQ=255)
   agTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCt                                                                                                                                                   <  1:1735983/150‑1 (MQ=255)
    gTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTg                                                                                                                                                  <  2:843996/150‑1 (MQ=255)
    gTCGGTGCCCAGCAGACCAACCACTGACGTGACGCACGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAAAGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTg                                                                                                                                                  >  1:825720/1‑150 (MQ=255)
      cGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGcc                                                                                                                                                <  2:1564611/150‑1 (MQ=255)
       ggTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGTCGAATTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCATCGCCCCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGCTCTGccc                                                                                                                                               >  1:535283/1‑150 (MQ=255)
         ggCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATCCTCTGCCCag                                                                                                                                             <  2:466279/149‑2 (MQ=255)
          gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAACCCTGGGCAGAGGATCTGCCCACt                                                                                                                                            >  1:616400/1‑150 (MQ=255)
           cccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGccacctg                                                                                                                                           >  2:1594177/1‑145 (MQ=255)
            ccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgg                                                                                                                                          >  1:301022/1‑149 (MQ=255)
            ccAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTga                                                                                                                                          >  1:22754/1‑150 (MQ=255)
            ccAGCAGACCAACCACTGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTga                                                                                                                                          <  1:382829/150‑1 (MQ=255)
             cagcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgag                                                                                                                                         >  2:1496152/1‑150 (MQ=255)
             cagcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCACAGGATCTGCCCACTgag                                                                                                                                         >  2:1510006/1‑150 (MQ=255)
              agcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgaga                                                                                                                                        <  1:1526067/150‑1 (MQ=255)
                cagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTAAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGGTGATCAATAACACCTGGGCAGAGGATCTGCCCACTGAGATc                                                                                                                                      <  1:738305/150‑1 (MQ=255)
                 agACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGATGATCTGCCACCTGAGATCg                                                                                                                                     >  2:136902/1‑150 (MQ=255)
                       accaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAAcc                                                                                                                               >  2:270056/1‑150 (MQ=255)
                        ccaCTGCCCTGCCGCCCGCCTCCGGCAGGCGACTTCGCGCCACTTCCGGCGGGCGCGGCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACCACCg                                                                                                                              <  2:1767501/150‑1 (MQ=255)
                         caCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGcc                                                                                                                                                >  2:1801422/1‑131 (MQ=255)
                          aCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTg                                                                                                                            >  2:382140/1‑150 (MQ=255)
                          aCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTg                                                                                                                            >  1:1311539/1‑150 (MQ=255)
                           cTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGc                                                                                                                           <  1:1603290/150‑1 (MQ=255)
                             gacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAg                                                                                                                         >  1:999304/1‑150 (MQ=255)
                             gacgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAg                                                                                                                         <  2:1582578/150‑1 (MQ=255)
                               cgtgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGtt                                                                                                                       >  2:748862/1‑150 (MQ=255)
                                    cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                                  >  1:871290/1‑150 (MQ=255)
                                    cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                                  >  1:136006/1‑150 (MQ=255)
                                    cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                                  >  1:446555/1‑150 (MQ=255)
                                    cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGGGTGGATCAAGAAAACCTGGGCAGAGGATCTGCCCACTGAGAGCGCCAACCGTGCAGTTCGGTa                                                                                                                  >  1:1045624/1‑150 (MQ=255)
                                     gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTAc                                                                                                                 <  2:734066/150‑1 (MQ=255)
                                       ccGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACaa                                                                                                               <  2:1059773/150‑1 (MQ=255)
                                        cGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAAt                                                                                                              <  2:34097/150‑1 (MQ=255)
                                          cTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGt                                                                                                            <  2:586812/150‑1 (MQ=255)
                                           ttCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTc                                                                                                           <  1:642787/150‑1 (MQ=255)
                                                 gaGGCGACTTAGCGCGACTTCCGGCGTGCGCGTCGTGGGACCGGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATATAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAgg                                                                                                     <  1:1696591/149‑1 (MQ=255)
                                                 cAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAgg                                                                                                     >  1:870628/1‑150 (MQ=255)
                                                      gACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGatat                                                                                                >  2:637432/1‑150 (MQ=255)
                                                        cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCTATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                              >  2:923024/1‑150 (MQ=255)
                                                        cTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTg                                                                                              >  1:264273/1‑150 (MQ=255)
                                                          tAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCt                                                                                            >  2:1438567/1‑150 (MQ=255)
                                                           aGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGACAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTg                                                                                           <  1:299577/150‑1 (MQ=255)
                                                             cgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGc                                                                                         >  1:572406/1‑150 (MQ=255)
                                                             cgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGc                                                                                         <  1:1400739/150‑1 (MQ=255)
                                                                   ttCCGGCGTTAGCGGCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGCTCTGACCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGc                                                                                   <  2:154233/150‑1 (MQ=255)
                                                                      cGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCTATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATTTTGCTGGCAACGGCGAt                                                                                <  1:1263731/150‑1 (MQ=255)
                                                                       ggCGTGCACGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATa                                                                               <  2:577463/150‑1 (MQ=255)
                                                                             cgcgTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTg                                                                         >  2:1295644/1‑150 (MQ=255)
                                                                                gtcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCg                                                                      <  2:109543/150‑1 (MQ=255)
                                                                                gtcgtGGGACCCGCTTCGCCGCCACCGCCAATCAAATGGCCGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCg                                                                      <  1:1006238/150‑1 (MQ=255)
                                                                                   gtGGGACCTGCTTCGCCGCCACCGCCAATCAAATCGACGTGTTGATCAATAAACCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACACTGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTa                                                                   <  2:111129/150‑1 (MQ=255)
                                                                                      ggACCTGCTTCGCCGCAACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTTGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGGACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcg                                                                <  2:1642860/150‑1 (MQ=255)
                                                                                       gACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcga                                                               >  2:875395/1‑150 (MQ=255)
                                                                                        aCCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgac                                                              >  1:1027908/1‑150 (MQ=255)
                                                                                         ccTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgacg                                                             >  2:1221472/1‑150 (MQ=255)
                                                                                            ccTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAgg                                                          <  2:128610/149‑1 (MQ=255)
                                                                                                cgccgcCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACAt                                                      <  2:823688/150‑1 (MQ=255)
                                                                                                 gccgccACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATc                                                     >  2:748729/1‑150 (MQ=255)
                                                                                                  ccgccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCg                                                    >  1:770116/1‑150 (MQ=255)
                                                                                                    gccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCa                                                  <  2:1715472/150‑1 (MQ=255)
                                                                                                        ccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAtt                                              >  1:165450/1‑150 (MQ=255)
                                                                                                        ccgcaaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAtt                                              >  1:609758/1‑150 (MQ=255)
                                                                                                          gccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTcc                                            >  1:1479464/1‑150 (MQ=255)
                                                                                                          gccaATCAAATGGACGTGTTGATCAAAAAAACCTGGGCAGAGGATCAGCCCACTGATATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTcc                                            <  1:427681/150‑1 (MQ=255)
                                                                                                           ccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTccc                                           <  1:1252579/150‑1 (MQ=255)
                                                                                                            caATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTcacc                                          >  1:484981/1‑150 (MQ=255)
                                                                                                                cataTGGACGTGTTGATCAATAAAAGCTCGGCAGAGGATCTGCCCACTGAGATCGACACCCGTGCTGTTCGGTACAATGTCAGAATGGATATTGCTGGCTACGGCGATATTTTTGCAGTTAGCGACGAGGACATCGCAAATTCCCCGATc                                      <  1:988184/147‑1 (MQ=255)
                                                                                                                cAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATc                                      >  1:44125/1‑150 (MQ=255)
                                                                                                                cAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAAGTCCCCGATc                                      >  2:715191/1‑150 (MQ=255)
                                                                                                                cAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCAGTTAGCGACGAGGACATCGCAAATTCCCCGATc                                      <  2:1248347/150‑1 (MQ=255)
                                                                                                                 aaaTGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCt                                     >  1:1726857/1‑150 (MQ=255)
                                                                                                                    tGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTGCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTTGCGACGAGGACATCGCAAATTCCCCGATCTTcc                                  >  2:916613/1‑150 (MQ=255)
                                                                                                                      gACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCgg                                >  2:1215332/1‑150 (MQ=255)
                                                                                                                       aCGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATAGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACAGCGCAAATTCCCCGATCTTCCGgc                               >  1:1146684/1‑150 (MQ=255)
                                                                                                                        cGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcg                              >  2:209599/1‑150 (MQ=255)
                                                                                                                          ggtTGATCAATAAACCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCa                            <  2:351643/149‑1 (MQ=255)
                                                                                                                              gATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACa                        >  1:1631679/1‑150 (MQ=255)
                                                                                                                                 cAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGCCGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt                     >  2:29714/1‑150 (MQ=255)
                                                                                                                                 cAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt                     >  2:700880/1‑150 (MQ=255)
                                                                                                                                 cAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGCGATATTGCTGGCAACGGCGATAATTTTGCCCTTAGCGACGAGGCCATCGCAAATTGCCCGATCTTCCGGCGCATACAAAt                     >  1:1381361/1‑150 (MQ=255)
                                                                                                                                  aaTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTTCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATg                    <  1:915639/150‑1 (MQ=255)
                                                                                                                                   aTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGCAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACAGCGCAAATTGCCCGATCTTCCGGCGCATACAAATgc                   >  2:907327/1‑150 (MQ=255)
                                                                                                                                        aCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTcc              <  2:795954/150‑1 (MQ=255)
                                                                                                                                         ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTccc             <  2:301022/150‑1 (MQ=255)
                                                                                                                                         ccTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTccc             >  1:787977/1‑150 (MQ=255)
                                                                                                                                          gtGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCt            >  1:149013/2‑150 (MQ=255)
                                                                                                                                          gtGGGCAGATCCTCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAAGCGCTCCCt            >  1:652970/2‑150 (MQ=255)
                                                                                                                                           tGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCtt           >  2:1838518/1‑149 (MQ=255)
                                                                                                                                           tGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTg           <  2:210215/150‑1 (MQ=255)
                                                                                                                                             ggCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgca         <  2:454006/150‑1 (MQ=255)
                                                                                                                                              gCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCACCGGCGATAATTTTGCCGTTTGCGCCGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcag        >  2:1633984/1‑150 (MQ=255)
                                                                                                                                              gCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcag        <  1:1052549/150‑1 (MQ=255)
                                                                                                                                                agagGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTgcagca      <  1:666973/150‑1 (MQ=255)
                                                                                                                                                   ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATGCAAATGCGCTCCCTGCAGCAggg   >  1:1796070/1‑150 (MQ=255)
                                                                                                                                                   ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAggg   <  1:1228537/150‑1 (MQ=255)
                                                                                                                                                   ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAggg   <  2:1851082/150‑1 (MQ=255)
                                                                                                                                                   ggATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAggg   <  1:1594177/150‑1 (MQ=255)
                                                                                                                                                    gATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGCCGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGGGCTCCCTGCAGCAGGGt  >  2:143105/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGCTCCCTGCAGCAGGGT  >  NC_000913/4559202‑4559499

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