Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,397,505 G→A 6.4% E32K (GAA→AAA)  mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,397,5050GA 6.4% 379.8 / 10.5 111E32K (GAA→AAA) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (57/46);  new base A (4/3);  total (62/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

TACGGTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAG  >  NC_000913/4397366‑4397644
                                                                                                                                           |                                                                                                                                           
tACGGTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTgg                                                                                                                                             >  2:332512/1‑140 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCTATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTAGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         <  2:773278/140‑1 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACTACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         >  2:662747/1‑140 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         >  2:1148451/1‑140 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         >  2:1148486/1‑140 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         <  1:1657503/140‑1 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         >  2:201280/1‑140 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         <  1:143821/140‑1 (MQ=255)
    gTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGaaaa                                                                                                                                         <  2:378491/140‑1 (MQ=255)
       acgacgCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAg                                                                                                                                      >  2:109824/1‑140 (MQ=255)
          acgCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCATAGAACTAGTGGAAAACAGCCt                                                                                                                                   <  2:822854/140‑1 (MQ=255)
          acgCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCt                                                                                                                                   <  2:566749/140‑1 (MQ=255)
          acgCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCt                                                                                                                                   >  1:871537/1‑140 (MQ=255)
          acgCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCt                                                                                                                                   >  1:375076/1‑140 (MQ=255)
          acgCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCt                                                                                                                                   <  1:1176639/140‑1 (MQ=255)
             ccAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGa                                                                                                                                >  1:1675488/1‑140 (MQ=255)
             ccAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGa                                                                                                                                <  2:1462629/140‑1 (MQ=255)
             ccAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGa                                                                                                                                <  2:1462741/140‑1 (MQ=255)
             ccAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGa                                                                                                                                <  2:1462635/140‑1 (MQ=255)
               aGATCGCACGCTGCAAAACTAAGCACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATg                                                                                                                              <  2:1064236/140‑1 (MQ=255)
                gcattggcgatgagagaggaggattacaaaATGCCAATTCAGGTCTTACCGCCACAACTTGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGAAAAACAGCCTCGATGc                                                                                                                             <  1:1235413‑M2/110‑1 (MQ=255)
                     cACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTg                                                                                                                        >  1:571652/1‑140 (MQ=255)
                     cACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTg                                                                                                                        >  1:570581/1‑140 (MQ=255)
                      aCGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCTAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGc                                                                                                                       >  2:569821/1‑140 (MQ=255)
                      aCGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGACCCAGATTGCCGCAGGTGAGGTGGGCGAGCGACCTGCGTCGGTAGTCAAAGGACTAGTGGAAAACAGCCTCGATGCAGGTGc                                                                                                                       >  2:699125/1‑140 (MQ=255)
                      aCGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGAGTGCCGCAGGGGAGGTGGTCGAGCGACCTGCGGCGGTAGTCACAGACCGAGTGGAAAACAGCCTCGATGCAGGTGc                                                                                                                       >  2:699449/1‑140 (MQ=255)
                          tGCCAACCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGAcg                                                                                                                   <  1:320948/140‑1 (MQ=255)
                          tGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGAcg                                                                                                                   <  1:1163433/140‑1 (MQ=255)
                          tGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGAcg                                                                                                                   <  1:321093/140‑1 (MQ=255)
                          tGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGAcg                                                                                                                   <  1:895954/140‑1 (MQ=255)
                          tGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGAcg                                                                                                                   <  2:642240/140‑1 (MQ=255)
                                aCTAAGGACGTATGAAGCCAATTCAGGTCTTACCGCCACAACAGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:1680187/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:1196936/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:1196854/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:752093/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:874797/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:668668/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:511528/140‑1 (MQ=255)
                                aCTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATc                                                                                                             <  2:511874/140‑1 (MQ=255)
                                  tAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGa                                                                                                           <  2:1533605/140‑1 (MQ=255)
                                  tAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGa                                                                                                           >  1:300943/1‑140 (MQ=255)
                                    aGGACGATTGATGCGAATTCAGGTCTTACCGCCACAACTGGCGTACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAATGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGata                                                                                                         >  2:1637522/1‑140 (MQ=255)
                                    aGGACGATTGATGCCATTTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGGCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGAATCGata                                                                                                         >  2:1429521/1‑140 (MQ=255)
                                    aGGACGATTGATGCCATTTCAGGTCTTACCGACACAACTGGCGAACCAGATTGCCGCAGGTCAGGTGGTCGAGCGACCTGCGTCGGTAGTCCATGAACTAGTGGAAAACAGCCACGATGCAGGTGCGACGCGTATCGatc                                                                                                         >  2:1429532/1‑139 (MQ=255)
                                    aGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGata                                                                                                         <  1:109824/140‑1 (MQ=255)
                                    aGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGata                                                                                                         >  2:783998/1‑140 (MQ=255)
                                    aGGACGATTGATGACAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCTCAGGTGAGGTGGTCGAGCTACCTGCGTCGGTAGTCAAATAACTATTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGata                                                                                                         >  2:777899/1‑140 (MQ=255)
                                      gACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATAtt                                                                                                       >  1:1145271/1‑140 (MQ=255)
                                      gACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATAtt                                                                                                       <  1:228548/140‑1 (MQ=255)
                                      gACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATAtt                                                                                                       >  1:1145393/1‑140 (MQ=255)
                                      gACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATAtt                                                                                                       <  1:1282722/140‑1 (MQ=255)
                                        tacaaaATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGa                                                                                                     >  2:531111‑M2/7‑140 (MQ=255)
                                          aTTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGata                                                                                                   >  1:428226/1‑140 (MQ=255)
                                             aATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCg                                                                                                >  2:1330047‑M2/2‑140 (MQ=255)
                                              aTGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGa                                                                                               >  2:1371834/1‑140 (MQ=255)
                                              aTGCCAAATCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGa                                                                                               >  2:1594159/1‑140 (MQ=255)
                                                gCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAAc                                                                                             >  1:82537/1‑140 (MQ=255)
                                                gCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTAGGTAGTCAAAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAAc                                                                                             >  2:3701/1‑140 (MQ=255)
                                                     ttCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGt                                                                                        >  2:1328439/1‑140 (MQ=255)
                                                     ttCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGt                                                                                        >  1:556003/1‑140 (MQ=255)
                                                     ttCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGt                                                                                        >  1:1661729/1‑140 (MQ=255)
                                                     ttCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGt                                                                                        >  1:1661713/1‑140 (MQ=255)
                                                     ttCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGt                                                                                        >  1:1567327/1‑140 (MQ=255)
                                                     ttCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGGGGTCGAGCGACCTGCGTCGGTAGTCACAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATTTTGCTATCGAACGCGGt                                                                                        <  1:974624/140‑1 (MQ=255)
                                                     atCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGACGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGt                                                                                        >  1:1567299/2‑140 (MQ=255)
                                                               aCCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAAc                                                                              >  2:112612/1‑140 (MQ=255)
                                                                   ccACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTAt                                                                          <  1:930785/140‑1 (MQ=255)
                                                                    cacaACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATc                                                                         >  1:1401355/1‑140 (MQ=255)
                                                                    cacaACTGGCGAACCAGATTGCCGCAGGTGAGGTGGGCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATc                                                                         >  1:1330340/1‑140 (MQ=255)
                                                                        aCTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCa                                                                     <  1:1169842/140‑1 (MQ=255)
                                                                        aCTGGAGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTAAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCa                                                                     <  1:1297306/140‑1 (MQ=255)
                                                                           ggCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGTAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTc                                                                  <  1:156586/140‑1 (MQ=255)
                                                                           ggCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTc                                                                  <  1:153742/140‑1 (MQ=255)
                                                                                 ccAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATa                                                            >  1:208825/1‑140 (MQ=255)
                                                                                          cGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGAGGAAAAAAGCCTCGATGCAGGTGCGACGCGTATCGATATTGAAATCGAACGCTGAGGGGCGAAACTTATCCGCATACGTCATAACGGCTGCg                                                   >  2:185980/1‑140 (MQ=255)
                                                                                              ggtgaggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATAATGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTAt                                               >  2:1612477/1‑140 (MQ=255)
                                                                                                  aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                           <  1:596947/140‑1 (MQ=255)
                                                                                                  aggtgGTCGAGCGACCTGCGTCGGTAGTCAAAGAAATAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATTTTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaa                                           <  1:1683348/140‑1 (MQ=255)
                                                                                                    gtggtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCaaaaa                                         >  2:1428686/1‑140 (MQ=255)
                                                                                                       gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                      >  1:1423976/1‑140 (MQ=255)
                                                                                                       gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                      >  1:275116/1‑140 (MQ=255)
                                                                                                       gtCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGa                                      >  1:1162188/1‑140 (MQ=255)
                                                                                                             cGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgct                                <  2:850681/140‑1 (MQ=255)
                                                                                                               aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                              >  1:335373/1‑140 (MQ=255)
                                                                                                               aCCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctgg                              >  1:1153681/1‑140 (MQ=255)
                                                                                                                    cGTCGGTAGTCAAAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATAGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGATGCGGTATCAAAAAAGATGAgctggcgctg                         >  2:1579937/1‑140 (MQ=255)
                                                                                                                          tAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAgctggcgctggcgctg                   >  2:1545119/1‑140 (MQ=255)
                                                                                                                              cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc               <  2:206766/140‑1 (MQ=255)
                                                                                                                              cAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc               <  2:1478198/140‑1 (MQ=255)
                                                                                                                              cAAAGAACTAGTGAAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTc               <  2:311204/140‑1 (MQ=255)
                                                                                                                                aaGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGt             <  2:1250543/140‑1 (MQ=255)
                                                                                                                                  gAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCa           <  1:720394/140‑1 (MQ=255)
                                                                                                                                    aCTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATg         >  1:1319494/1‑140 (MQ=255)
                                                                                                                                       aGTGGAAATCAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca      >  2:1402293/1‑140 (MQ=255)
                                                                                                                                       aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca      >  1:1334265/1‑140 (MQ=255)
                                                                                                                                       aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca      >  1:1334287/1‑140 (MQ=255)
                                                                                                                                       aGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGcca      >  1:800754/1‑140 (MQ=255)
                                                                                                                                        gTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGccac     >  2:1284390/1‑140 (MQ=255)
                                                                                                                                           gAAAACAGCCTCGATGCAGGTGCGACGCGTTTCGATATTGATATCGAACGCGGTGGTGCGAAACTTATCCGCATTCGTGATAACGGATGCGGTATCAAAAAAGTTGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg  >  2:1222655/1‑140 (MQ=255)
                                                                                                                                           gAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGAGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAg  >  2:1222239/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TACGGTGACGACGCCAGATCGCACGCTGCCAAACTAAGGACGATTGATGCCAATTCAGGTCTTACCGCCACAACTGGCGAACCAGATTGCCGCAGGTGAGGTGGTCGAGCGACCTGCGTCGGTAGTCAAAGAACTAGTGGAAAACAGCCTCGATGCAGGTGCGACGCGTATCGATATTGATATCGAACGCGGTGGGGCGAAACTTATCCGCATTCGTGATAACGGCTGCGGTATCAAAAAAGATGAGCTGGCGCTGGCGCTGGCTCGTCATGCCACCAG  >  NC_000913/4397366‑4397644

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