Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 835,211 G→A 5.9% R714R (CGG→CGA dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913835,2110GA 5.9% 266.8 / 5.5 85R714R (CGG→CGAdinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (46/34);  new base A (2/3);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

AGTGGTTATCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAATGGTGAAGTCC  >  NC_000913/835067‑835359
                                                                                                                                                |                                                                                                                                                    
aGTGGTTATCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCggcggc                                                                                                                                                  <  2:168177/149‑1 (MQ=255)
 gTGGTTATCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGgcggcg                                                                                                                                                 >  2:4021/1‑149 (MQ=255)
     ttATCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTaa                                                                                                                                             >  4:27652/1‑149 (MQ=255)
     ttATCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGCCGACGGCGtggg                                                                                                                                             <  3:504258/149‑4 (MQ=255)
       aTCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGCCGACGGCGtgggga                                                                                                                                           <  3:144589/149‑6 (MQ=255)
       aTCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGCCGACGGCGtgggga                                                                                                                                           >  4:144589/1‑144 (MQ=255)
        tCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatga                                                                                                                                          <  4:482536/149‑1 (MQ=255)
        tCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGAATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatga                                                                                                                                          >  1:103688/1‑149 (MQ=255)
         cTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatgat                                                                                                                                         <  4:17782/149‑1 (MQ=255)
          tACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatgatg                                                                                                                                        >  2:574575/1‑149 (MQ=255)
          tACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatgatg                                                                                                                                        <  2:77076/149‑1 (MQ=255)
             gACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGa                                                                                                                                     <  4:244512/149‑1 (MQ=255)
             gACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGa                                                                                                                                     <  4:22622/149‑1 (MQ=255)
                        gCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAgg                                                                                                                          >  4:180103/1‑149 (MQ=255)
                            aCCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTc                                                                                                                      <  1:425356/149‑1 (MQ=255)
                             ccAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGCCGACGGCGtgg                                                                                                                                              <  3:155139/124‑3 (MQ=255)
                             ccAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGCCGACGGCGtgg                                                                                                                                              >  4:155139/1‑122 (MQ=255)
                              cAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGt                                                                                                                    >  2:299263/1‑149 (MQ=255)
                               aaaaaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatga                                                                                                                                          <  3:404651/126‑1 (MQ=255)
                               aaaaaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAatga                                                                                                                                          >  4:404651/1‑126 (MQ=255)
                               aaaaaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTa                                                                                                                   >  2:569874/1‑149 (MQ=255)
                               aaaaaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTa                                                                                                                   >  2:459818/1‑149 (MQ=255)
                               aaaaaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTa                                                                                                                   >  1:147038/1‑149 (MQ=255)
                                 aaaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTaaa                                                                                                                 <  1:464936/149‑1 (MQ=255)
                                  aaTTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAAt                                                                                                                <  1:574575/149‑1 (MQ=255)
                                           aaGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCg                                                                                                       >  3:104682/1‑149 (MQ=255)
                                               gACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCaaa                                                                                                   >  1:9758/1‑149 (MQ=255)
                                                  tactGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATc                                                                                                >  1:357443/1‑149 (MQ=255)
                                                    ctGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCAt                                                                                              >  2:448115/1‑149 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGa                                                                                                             >  3:557365/1‑131 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGa                                                                                                             <  4:557365/131‑1 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCa                                                                                               >  3:469572/1‑145 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCa                                                                                               <  4:469572/145‑1 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTaa                                                                                           >  2:414656/1‑149 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTaa                                                                                           >  3:115350/1‑149 (MQ=255)
                                                       gATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGACCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTaa                                                                                           <  3:314018/149‑1 (MQ=255)
                                                        aTGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTaaa                                                                                          <  3:420225/149‑1 (MQ=255)
                                                             ttACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGAt                                                                                     >  4:238145/1‑149 (MQ=255)
                                                              tACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATc                                                                                    <  1:585917/149‑1 (MQ=255)
                                                                 cGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCggg                                                                                 >  2:509568/1‑149 (MQ=255)
                                                                   gTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAaga                                                                                        <  2:584333/140‑1 (MQ=255)
                                                                   gTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAaga                                                                                        >  1:584333/1‑140 (MQ=255)
                                                                   gTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGgcg                                                                               >  2:577511/1‑149 (MQ=255)
                                                                     aTTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGcgcg                                                                             <  1:289950/149‑1 (MQ=255)
                                                                            aTAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGaaaaa                                                                      <  1:187816/149‑1 (MQ=255)
                                                                                agCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCAc                                                                  >  2:433529/1‑149 (MQ=255)
                                                                                            tCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTa                                                      <  2:147038/149‑1 (MQ=255)
                                                                                            tCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTa                                                      <  3:308332/149‑1 (MQ=255)
                                                                                               cTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGAc                                                   >  1:137131/1‑149 (MQ=255)
                                                                                               cTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGAc                                                   <  2:9758/149‑1 (MQ=255)
                                                                                                tGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGAcc                                                  >  1:416680/1‑149 (MQ=255)
                                                                                                   aGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCggg                                               <  1:4021/149‑1 (MQ=255)
                                                                                                    ggTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGa                                              <  1:299263/149‑1 (MQ=255)
                                                                                                      tataGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATa                                            >  2:473590/1‑149 (MQ=255)
                                                                                                       ataGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATAc                                           >  1:381105/1‑149 (MQ=255)
                                                                                                       ataGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATAc                                           >  1:4717/1‑149 (MQ=255)
                                                                                                        taGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGAGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACt                                          >  1:268011/1‑149 (MQ=255)
                                                                                                          gTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTgc                                        >  4:250306/1‑149 (MQ=255)
                                                                                                          gTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTgc                                        >  4:527215/1‑149 (MQ=255)
                                                                                                           ttAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTgcg                                       >  2:349011/1‑149 (MQ=255)
                                                                                                           ttAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTgcg                                       >  2:14471/1‑149 (MQ=255)
                                                                                                            tAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTgcgc                                      >  4:506923/1‑149 (MQ=255)
                                                                                                             aaaaaGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTgcgcg                                     >  4:364689/1‑149 (MQ=255)
                                                                                                                aagaaagaAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCgg                                  >  3:329471/1‑149 (MQ=255)
                                                                                                                   aaagaaaAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCt                               <  4:104682/149‑1 (MQ=255)
                                                                                                                     agaaaAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCt                               >  3:65741/1‑147 (MQ=255)
                                                                                                                     agaaaAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCt                               <  4:430697/147‑1 (MQ=255)
                                                                                                                     agaaaAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCt                               <  4:65741/147‑1 (MQ=255)
                                                                                                                     agaaaAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCt                               >  3:430697/1‑147 (MQ=255)
                                                                                                                     agaaaAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGt                             <  2:137131/149‑1 (MQ=255)
                                                                                                                       aaaaaaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATgg                                                                                                              <  2:231873/66‑1 (MQ=255)
                                                                                                                       aaaaaaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATgg                                                                                                              >  1:231873/1‑66 (MQ=255)
                                                                                                                       aaaaaaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCg                                                                                   >  3:62228/1‑93 (MQ=255)
                                                                                                                       aaaaaaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCg                                                                                   <  4:62228/93‑1 (MQ=255)
                                                                                                                       aaaaaaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGtat                           >  4:30031/1‑149 (MQ=255)
                                                                                                                           aaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGtatgctc                       >  1:285916/1‑149 (MQ=255)
                                                                                                                           aaCGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCAGTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGtatgctc                       >  4:308476/1‑149 (MQ=255)
                                                                                                                              gAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGtatgctcata                    >  1:468936/1‑149 (MQ=255)
                                                                                                                               aaaTCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGtatgctcatat                   <  3:27652/149‑1 (MQ=255)
                                                                                                                                    cccACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGtatgctcatatgctca              >  3:578850/1‑149 (MQ=255)
                                                                                                                                     ccACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCaa             >  4:562000/1‑149 (MQ=255)
                                                                                                                                      cACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAAt            <  1:569874/149‑1 (MQ=255)
                                                                                                                                      cACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAAt            <  2:103688/149‑1 (MQ=255)
                                                                                                                                      cACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAAt            <  2:448316/149‑1 (MQ=255)
                                                                                                                                       aCGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAATg           >  3:503407/1‑149 (MQ=255)
                                                                                                                                       aCGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAATg           <  3:30031/149‑1 (MQ=255)
                                                                                                                                         gCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAATGGt         <  1:509568/149‑1 (MQ=255)
                                                                                                                                                gcggcgTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAATGGTGAAGTcc  <  2:357443/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGTGGTTATCTACGACAAACGCTTGCTGACCAAAAATTATGGCAAGCGACTACTGGATGCATTACCGGTATTTCCGATAGAGCAACCGGAAGTCCCTGAAGGTATAGTTAAAAAGAAAGAAAAAACGAAATCCCCACGCCGTCGGCGGCGTTAATGATGTGAGTCAGGTAAGGAGTCGTAAATGGACTATCGCAAAATCATTAAAGAGATCGGGCGCGGGAAAAACCACGCGCGCGATTTAGACCGGGATACTGCGCGCGGTCTGTATGCTCATATGCTCAATGGTGAAGTCC  >  NC_000913/835067‑835359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: