Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,515,785 G→T 6.2% L43L (CTC→CTA) ‡ dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,7850GT 6.2% 348.0 / 11.2 129L43L (CTC→CTA) ‡damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (60/60);  new base T (5/3);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

Only 100 of 135 total aligned reads displayed.

TGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCT  >  NC_000913/3515636‑3515934
                                                                                                                                                     |                                                                                                                                                     
tGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTg                                                                                                                                                       <  1:929867/150‑1 (MQ=255)
  aTAGTAAACCGCGGCGCACTTTGGTTCGGGAACAACCAGCGCGCGTGCGGCCGGGACGTACTCCTCCGTACGCATCTTCACAATGTTATAGAGACTGATCAGGGCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                                                                     <  2:836264/150‑1 (MQ=255)
    aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATTGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                   >  2:707323/1‑150 (MQ=255)
     gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTAATGATATCGTCAAGGATATAACGAGAAAAGTCGGTGTTGAGaaa                                                                                                                                                  >  2:1040482/1‑150 (MQ=255)
       aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                >  2:918823/1‑150 (MQ=255)
          ccTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCg                                                                                                                                             <  2:984833/150‑1 (MQ=255)
             cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGAct                                                                                                                                          <  1:972265/150‑2 (MQ=255)
              ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGActt                                                                                                                                         <  2:602974/150‑3 (MQ=255)
              ggCGAAATTTGTTTCGGTAACAAACAGCTCGCGTGCGGCCTGTACGTACTCTTCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGActt                                                                                                                                         <  2:1073617/150‑3 (MQ=255)
               gcgcAATTTGTTTCGGTAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCg                                                                                                                                        <  1:1521007/150‑1 (MQ=255)
                cgcAATTTGTTTCGGGAACAAACAGCTCCCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTAAAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGATAAACACAGACACgg                                                                                                                                       >  1:1562908/1‑150 (MQ=255)
                  cAATTTGTTTCGGTAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                                                                                                                                     <  1:896263/150‑1 (MQ=255)
                   aaTTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGATATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCAc                                                                                                                                    >  2:966467/1‑150 (MQ=255)
                     tttGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCCGTACGCATCTTCACAAGGTAATAGAGCCTGCTCAGGTGGCGTTTGCTATCGGCAAGGATATAACGAGAAAAGTCGGTGTTCAGAAACACCGACCCgggccct                                                                                                                                  >  2:274231/1‑145 (MQ=255)
                            cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                                                                      <  2:1598788/123‑1 (MQ=255)
                            cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                                                                           >  2:1733968/1‑150 (MQ=255)
                            cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                                                                           >  2:873405/1‑150 (MQ=255)
                            cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGGCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                                                                           <  2:3707/150‑1 (MQ=255)
                              ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                         <  1:190197/150‑1 (MQ=255)
                              ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                         <  1:354145/150‑1 (MQ=255)
                               gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCt                                                                                                                        >  1:1434340/1‑150 (MQ=255)
                                aacaaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTc                                                                                                                       >  2:1217336/1‑150 (MQ=255)
                                    aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                                                                   <  1:1601064/150‑1 (MQ=255)
                                     acaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCa                                                                                                                  <  2:37441/150‑1 (MQ=255)
                                          tCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                                                                             >  1:931148/1‑150 (MQ=255)
                                             cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                          >  1:224737/1‑150 (MQ=255)
                                              gTGCGGCGTGTACGTAGTCATAAGTAGGCATATTCAAAATGTTATTGAGACTGTTCTGGTGGCTTTTTATATCGGCAATGATATAACGACAAACGGCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGc                                                                                                         <  2:1307958/150‑1 (MQ=255)
                                              gTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGc                                                                                                         <  1:398498/150‑1 (MQ=255)
                                              gTGCGGCCTGTACGTAATAATCAGTACGAATCTTCAAAATGTTATAGAGACTGATCAGGTCGCTATTTATATCGGAAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGc                                                                                                         <  1:1697056/150‑1 (MQ=255)
                                                    ccTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGATCAACCAGACATTCGCCCTTgg                                                                                                   >  2:1310857/1‑150 (MQ=255)
                                                     cTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                                                                  <  1:833680/150‑1 (MQ=255)
                                                     cTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                                                                  <  1:918823/150‑1 (MQ=255)
                                                      tGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                                                                 >  2:1203566/1‑150 (MQ=255)
                                                         acgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGCAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTAGGGCaaa                                                                                              >  1:1338731/1‑150 (MQ=255)
                                                          cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTTGCCCAAGGGCGAAt                                                                                             <  1:415958/150‑1 (MQ=255)
                                                            tacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGc                                                                                           <  1:1187074/150‑1 (MQ=255)
                                                            tacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGGCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGc                                                                                           <  1:718060/150‑1 (MQ=255)
                                                              cTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCg                                                                                         >  1:1147064/1‑150 (MQ=255)
                                                               tcatcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGt                                                                                        >  2:1242416/1‑150 (MQ=255)
                                                                catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                       >  2:1498703/1‑150 (MQ=255)
                                                                catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCAATTGATATCGGCAAGGATATAACGAGCAAAGTCGGAGTTGAGAAACACCGACCCGGCACCGACAAAAGGCTCAACCAGACATTCGCACTTGGGCAAATGCCGtt                                                                                       >  1:164957/1‑150 (MQ=255)
                                                                  tcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTa                                                                                     >  1:1400681/1‑150 (MQ=255)
                                                                   caGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATAAAACGAGAAAAGTCGGTGTTGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa                                                                                    <  1:707863/150‑1 (MQ=255)
                                                                     gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTAGGGCAAATGCCGTTTAata                                                                                  >  2:1512694/1‑150 (MQ=255)
                                                                        cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                               <  1:258881/150‑1 (MQ=255)
                                                                        cGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatca                                                                               >  2:466450/1‑150 (MQ=255)
                                                                          cATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcatc                                                                             >  2:1448292/1‑150 (MQ=255)
                                                                          cATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcatc                                                                             <  2:717047/150‑1 (MQ=255)
                                                                              ttCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                                                                         <  1:768601/150‑1 (MQ=255)
                                                                              ttCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                                                                         <  2:777389/150‑1 (MQ=255)
                                                                              ttCACAATGTTATAGAGACTGATCAGGTCGATATTGATATCGGCAAGGATTTAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                                                                         <  2:658611/150‑1 (MQ=255)
                                                                                 acaATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                                                                      <  1:1586314/150‑1 (MQ=255)
                                                                                 aaaaTTTTATAGAGACTTATCAGGTCGCTTTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgg                                                                      <  1:376075/148‑1 (MQ=255)
                                                                                  caATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                                                                     >  1:1144584/1‑150 (MQ=255)
                                                                                            gagaCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAATTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc                                                           <  1:966467/150‑1 (MQ=255)
                                                                                            gagaCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc                                                           <  1:1217336/150‑1 (MQ=255)
                                                                                              gaCTGATCAGGGCGCTATTGATCTCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcccc                                                         <  2:103081/150‑1 (MQ=255)
                                                                                                 tGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTg                                                      >  1:1471443/1‑150 (MQ=255)
                                                                                                  gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc                                                     >  1:540145/1‑150 (MQ=255)
                                                                                                   aTCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc                                                    <  1:873405/150‑1 (MQ=255)
                                                                                                   aTCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGGTGAGAAACACCGAACCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc                                                    <  1:1297140/150‑1 (MQ=255)
                                                                                                        gTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACtt                                               >  2:275332/1‑150 (MQ=255)
                                                                                                         gcGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAAGGGATACTTGCCCCCTGCCCACTTc                                              <  2:1062664/149‑1 (MQ=255)
                                                                                                            cTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaa                                           >  2:197026/1‑150 (MQ=255)
                                                                                                             tATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaa                                          >  2:196234/1‑150 (MQ=255)
                                                                                                              aTTGATATCGGCAAGGATATAACGAGAAAAGTCTGAGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCAGTTTAATATCATCAAGCAGGGTATACTTGCCCCCTGCCCACTTCaaaaa                                         <  1:1246164/150‑1 (MQ=255)
                                                                                                              aTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaaa                                         >  2:514090/1‑150 (MQ=255)
                                                                                                                tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                       <  2:248149/150‑1 (MQ=255)
                                                                                                                tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                       <  2:1722538/150‑1 (MQ=255)
                                                                                                                 tatatCGGCATGTATATAACTCGAAAAGTCGTTGTTGCTAACCACCGACCCGGCACCTACAAAATGCGCATCCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                      <  1:863887/149‑1 (MQ=255)
                                                                                                                   tatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgcgc                                    >  2:1674834/1‑150 (MQ=255)
                                                                                                                    atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAACCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCATGGGCAAAAGCCCTTTAATATCATCAAGCAGGGTATACTTGCCCCCTGCACACTTCAAAAAAgcgag                                   >  2:875584/1‑148 (MQ=255)
                                                                                                                       ggCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCCCGCAGGGGATACTTGCCCCCTGCACACTTCAAAAAAGCGAGAtt                                >  2:1529048/1‑150 (MQ=255)
                                                                                                                         cAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAtttt                              >  2:342697/1‑150 (MQ=255)
                                                                                                                         cAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAtttt                              <  2:1147064/150‑1 (MQ=255)
                                                                                                                           actacataACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttc                            <  2:621265/145‑1 (MQ=255)
                                                                                                                             gATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttctt                          <  2:1154941/150‑1 (MQ=255)
                                                                                                                              atatAACGAGAAAAGTCGGTGTTGACAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttcttc                         <  1:637954/150‑1 (MQ=255)
                                                                                                                               tataACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCCGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCa                        >  2:340701/1‑150 (MQ=255)
                                                                                                                                 taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATg                      >  1:1650577/1‑150 (MQ=255)
                                                                                                                                  aactataaatGTCGGTGTGGATAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGc                     <  1:1314129/144‑1 (MQ=255)
                                                                                                                                  aaCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCACCCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCGGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGc                     >  2:1012972/1‑150 (MQ=255)
                                                                                                                                  aaCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGc                     <  2:758509/150‑1 (MQ=255)
                                                                                                                                      agaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGac                 <  2:1662419/150‑1 (MQ=255)
                                                                                                                                         cgggTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaa              >  1:1359526/4‑150 (MQ=255)
                                                                                                                                          aaGTCGGTGTTGAGAAACACCGACCAGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTGCAAAAAAGCCCGCTTTTTCTTCATGCTGactaac             >  2:278098/1‑150 (MQ=255)
                                                                                                                                           ggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact            >  2:1275578/2‑150 (MQ=255)
                                                                                                                                           aGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact            >  2:748906/1‑150 (MQ=255)
                                                                                                                                            gTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaacta           >  2:1454864/1‑150 (MQ=255)
                                                                                                                                            gTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAAAAACATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaacga           >  2:1535552/1‑148 (MQ=255)
                                                                                                                                            gTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTACAAAAAAGCGCGATTATTCTTCATGCTGactaacta           >  1:863892/1‑150 (MQ=255)
                                                                                                                                                gttttGATAAAAACCGAACCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCTCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa       <  2:1562908/147‑1 (MQ=255)
                                                                                                                                                gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa       >  2:1243508/1‑150 (MQ=255)
                                                                                                                                                gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa       >  1:1707856/1‑150 (MQ=255)
                                                                                                                                                  gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTaca     >  1:442338/1‑150 (MQ=255)
                                                                                                                                                  gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATATTTCTTCATGCTGACTAACTTATTaca     >  2:1113884/1‑150 (MQ=255)
                                                                                                                                                   ttGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGTGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTacac    >  2:818619/1‑150 (MQ=255)
                                                                                                                                                    tGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACAcc   >  1:385612/1‑150 (MQ=255)
                                                                                                                                                     gagaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCt  >  2:815130/1‑150 (MQ=255)
                                                                                                                                                     gagaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCt  >  1:515047/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCT  >  NC_000913/3515636‑3515934

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: