Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,515,787 G→T 5.0% L43I (CTC→ATC) ‡ dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,7870GT 5.0% 339.1 / 5.4 99L43I (CTC→ATC) ‡damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (48/46);  new base T (2/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 101 total aligned reads displayed.

AGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTC  >  NC_000913/3515640‑3515936
                                                                                                                                                   |                                                                                                                                                     
aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                     >  2:631382/1‑150 (MQ=255)
 gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaaa                                                                                                                                                    >  2:948557/1‑150 (MQ=255)
         cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAcc                                                                                                                                            <  2:957935/150‑1 (MQ=255)
             gcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGc                                                                                                                                        >  1:1992914/1‑150 (MQ=255)
               aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGActtttctc                                                                                                                                      <  1:1323473/150‑8 (MQ=255)
               aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGGCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGActtttctc                                                                                                                                      <  1:196938/150‑8 (MQ=255)
                        cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                                                                             >  1:634281/1‑150 (MQ=255)
                         gggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                                                                            <  1:1646990/150‑1 (MQ=255)
                         gggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                                                                            >  1:218104/1‑150 (MQ=255)
                         gggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                                                                            >  1:18055/1‑150 (MQ=255)
                          ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                           >  2:1390344/1‑150 (MQ=255)
                          ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                           <  2:116748/150‑1 (MQ=255)
                          cgAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                           >  2:1997889/2‑150 (MQ=255)
                             acaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCATGTCGCTATTGTTATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCa                                                                                                                        <  2:1480701/150‑1 (MQ=255)
                              caaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCaa                                                                                                                       >  2:1521925/1‑150 (MQ=255)
                                     cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACa                                                                                                                <  2:693551/150‑1 (MQ=255)
                                     cTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACa                                                                                                                >  2:923085/1‑150 (MQ=255)
                                      tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGCCTCAACCAGACAt                                                                                                               >  2:1291258/1‑150 (MQ=255)
                                        gcgTGCGGCCTTTACGTACTCATCAGTACTCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTc                                                                                                             <  1:1089162/150‑1 (MQ=255)
                                         cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                            >  2:1132415/1‑150 (MQ=255)
                                         cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                            >  2:720712/1‑150 (MQ=255)
                                          gTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGc                                                                                                           >  2:693030/1‑150 (MQ=255)
                                           tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGcc                                                                                                          >  2:745933/1‑150 (MQ=255)
                                           tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGcc                                                                                                          >  1:521570/1‑150 (MQ=255)
                                           tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGcc                                                                                                          >  2:136781/1‑150 (MQ=255)
                                               gCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTg                                                                                                      <  2:1144844/150‑1 (MQ=255)
                                               gCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTg                                                                                                      <  2:684225/150‑1 (MQ=255)
                                                ccTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTgg                                                                                                     <  2:1325321/150‑1 (MQ=255)
                                                 cTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                                                                    <  2:577924/150‑1 (MQ=255)
                                                  tGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                                                                   <  2:1886425/150‑1 (MQ=255)
                                                  gatgcgtacTGATGAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                                                                   >  1:1886284/5‑150 (MQ=255)
                                                   gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                                                                                                  <  2:1912253/150‑1 (MQ=255)
                                                   gtacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCa                                                                                                  <  2:387099/150‑1 (MQ=255)
                                                    tgcgtacTGATGAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCaa                                                                                                 >  2:1725788/3‑150 (MQ=255)
                                                       gtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATg                                                                                              >  2:388546/1‑150 (MQ=255)
                                                         acTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGcc                                                                                            <  1:369282/150‑1 (MQ=255)
                                                           tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGt                                                                                          >  2:51744/1‑150 (MQ=255)
                                                            catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                         <  2:1458636/150‑1 (MQ=255)
                                                            catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                         >  1:1644702/1‑150 (MQ=255)
                                                            catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                         >  2:1193308/1‑150 (MQ=255)
                                                             atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGttt                                                                                        <  2:335857/150‑1 (MQ=255)
                                                             atcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGttt                                                                                        <  2:1002164/150‑1 (MQ=255)
                                                               caGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa                                                                                      >  1:574892/1‑150 (MQ=255)
                                                                 gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata                                                                                    >  2:1984338/1‑150 (MQ=255)
                                                                   aCGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatc                                                                                  >  2:315208/1‑150 (MQ=255)
                                                                    cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                                 <  2:218104/150‑1 (MQ=255)
                                                                    cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                                 <  2:1644702/150‑1 (MQ=255)
                                                                    cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                                 <  2:1549392/150‑1 (MQ=255)
                                                                    cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                                 >  1:42091/1‑150 (MQ=255)
                                                                    cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                                 <  2:1422305/150‑1 (MQ=255)
                                                                          ttCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                                                                           >  2:890155/1‑150 (MQ=255)
                                                                              caATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                                                                       >  2:667611/1‑150 (MQ=255)
                                                                               aaTGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATATCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgggg                                                                      <  2:1592545/150‑1 (MQ=255)
                                                                                aTGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGa                                                                     <  2:1992914/150‑1 (MQ=255)
                                                                                aTGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGa                                                                     <  2:634281/150‑1 (MQ=255)
                                                                                   ttATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATAc                                                                  <  1:1521925/150‑1 (MQ=255)
                                                                                   ttATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATAc                                                                  >  2:1059090/1‑150 (MQ=255)
                                                                                    tataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGTCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACt                                                                 >  2:1978965/1‑150 (MQ=255)
                                                                                    tataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                                                                              >  2:1722480/1‑137 (MQ=255)
                                                                                    tataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                                                                              <  1:1722480/137‑1 (MQ=255)
                                                                                      taGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTg                                                               <  1:136781/150‑1 (MQ=255)
                                                                                       agagACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGc                                                              >  2:236561/1‑150 (MQ=255)
                                                                                            cTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCt                                                         >  1:47732/1‑150 (MQ=255)
                                                                                            cTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCt                                                         >  1:860325/1‑150 (MQ=255)
                                                                                             tGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg                                                        >  1:269198/1‑150 (MQ=255)
                                                                                              gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc                                                       >  1:1608892/1‑150 (MQ=255)
                                                                                               aTCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc                                                      >  1:76707/1‑150 (MQ=255)
                                                                                                   ggTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACt                                                  <  1:948557/150‑1 (MQ=255)
                                                                                                     tCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTc                                                <  2:574892/150‑1 (MQ=255)
                                                                                                       gCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaa                                              >  2:1496236/1‑150 (MQ=255)
                                                                                                             tatatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                        <  2:1236812/149‑1 (MQ=255)
                                                                                                             gATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                        >  1:1914920/1‑150 (MQ=255)
                                                                                                              atatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgcg                                       <  1:1390344/150‑1 (MQ=255)
                                                                                                                atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                     <  2:872856/150‑1 (MQ=255)
                                                                                                                        ggATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttct                             <  2:1223760/150‑1 (MQ=255)
                                                                                                                             taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                                                                        >  2:777497/1‑102 (MQ=255)
                                                                                                                             taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                                                                        <  1:777497/102‑1 (MQ=255)
                                                                                                                             taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                     >  1:599996/1‑137 (MQ=255)
                                                                                                                             taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                     <  2:599996/137‑1 (MQ=255)
                                                                                                                             taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATg                        >  2:1329544/1‑150 (MQ=255)
                                                                                                                              aaCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCTGGTGATACTTGCCCCCTTCCCACTTCAAAAAAGCGCGATTTTTCTTCTTgg                       <  2:798950/150‑2 (MQ=255)
                                                                                                                              aaCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGc                       >  1:1102305/1‑150 (MQ=255)
                                                                                                                               aCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCt                      <  1:1496236/150‑1 (MQ=255)
                                                                                                                               aCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCt                      <  1:1725788/150‑1 (MQ=255)
                                                                                                                                cGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTg                     >  2:1219480/1‑150 (MQ=255)
                                                                                                                                 gagaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTGATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGa                    <  1:1971776/150‑1 (MQ=255)
                                                                                                                                 gagaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGa                    >  2:876/1‑150 (MQ=255)
                                                                                                                                   gaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGact                  >  1:1041882/1‑150 (MQ=255)
                                                                                                                                   gaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGact                  <  2:1266249/150‑1 (MQ=255)
                                                                                                                                      gggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaac               >  1:1663379/3‑150 (MQ=255)
                                                                                                                                       aGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact              <  1:720712/150‑1 (MQ=255)
                                                                                                                                        gTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGact                  >  1:41807/1‑145 (MQ=255)
                                                                                                                                        gTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGact                  <  2:41807/145‑1 (MQ=255)
                                                                                                                                         tCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGActaactaa            <  2:269198/150‑1 (MQ=255)
                                                                                                                                         tCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGCCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGActaactaa            <  1:693030/150‑1 (MQ=255)
                                                                                                                                             tgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTac        >  1:1256353/1‑150 (MQ=255)
                                                                                                                                              gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTaca       <  1:745933/150‑1 (MQ=255)
                                                                                                                                               ttGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTacac      <  2:18055/150‑1 (MQ=255)
                                                                                                                                                  agaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCtt   >  1:1961835/1‑150 (MQ=255)
                                                                                                                                                   gaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTtc  >  1:914091/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTTC  >  NC_000913/3515640‑3515936

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: