Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007790 844 C→G Q189E (CAA→GAA)  SAUSA300_RS14705 → replication protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077908440CG89.8% 1598.2 / 153.1 655Q189E (CAA→GAA) SAUSA300_RS14705replication protein
Reads supporting (aligned to +/- strand):  ref base C (63/4);  new base G (551/37);  total (614/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-02

Only 100 of 655 total aligned reads displayed.

TTTTAACTGATGATGATGAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAG  >  NC_007790/746‑931
                                                                                                  |                                                                                       
agttAACTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGCGGATTACTTAAAGaa                                                                                       >  1:3338585/3‑100 (MQ=21)
 cttAACTGATGACGACGAAAAAAATCAAGAAATTTTACGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGAGGATTACTTAAAGaaa                                                                                      >  1:875542/2‑100 (MQ=21)
    taCTGATGACGACAAAAAAAATCAAGAAACTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGaaaagc                                                                                   >  1:557317/2‑98 (MQ=21)
    atcTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTATGTTATGGTGGATTACTTAAAGaaaagc                                                                                   >  1:1228187/3‑98 (MQ=25)
     aCTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCa                                                                                  >  1:4776408/1‑100 (MQ=25)
     aCTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCa                                                                                  >  1:3021319/1‑100 (MQ=25)
      cTGATGACGACGAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCAt                                                                                 >  1:5190448/1‑100 (MQ=25)
        tatagcgatgaAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATa                                                                                >  1:922305/7‑100 (MQ=25)
        gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:3790360/1‑100 (MQ=21)
        gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:803382/1‑100 (MQ=21)
        gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:3534454/1‑100 (MQ=21)
        gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:2721919/1‑100 (MQ=21)
        gatgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTAGTATCGAAAACGAATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:3928237/1‑100 (MQ=255)
        gatgacgacgaAAAAAATCAAGAAATTTTAAGAGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaa                                                                               >  1:4410877/1‑100 (MQ=21)
         ctaccgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaa                                                                              >  1:6032139/9‑100 (MQ=21)
         atgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaa                                                                              >  1:1294213/1‑100 (MQ=21)
          tgacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATaaaa                                                                             >  1:2141456/1‑100 (MQ=21)
          agacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGGATGTTAAGATATGGGGGATTACTTAAAGAAAAGCATaaaa                                                                             >  1:666560/8‑100 (MQ=255)
           gacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                            >  1:4569625/7‑100 (MQ=21)
           gacgacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                            >  1:2777537/7‑100 (MQ=21)
           gacgacgaAAAAAATCAAGAAATTTTAAGAGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAt                                                                            >  1:5043532/7‑100 (MQ=18)
              gacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAtttt                                                                         >  1:5858239/4‑100 (MQ=21)
              cacgaAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAAtttt                                                                         >  1:1457214/4‑100 (MQ=21)
                    aaaaaTCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACtt                                                                   <  1:847562/100‑1 (MQ=21)
                         tCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAGAAGCATAAAATTTTAAACTCAGACg                                                              <  1:4927043/100‑1 (MQ=21)
                          cAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGa                                                             >  1:4022919/1‑100 (MQ=21)
                                aTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGa                                                       <  1:1617645/100‑1 (MQ=21)
                                aTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGa                                                       <  1:3978035/100‑1 (MQ=21)
                                aTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGa                                                       <  1:2356436/100‑1 (MQ=21)
                                          gatctagaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTaa                                              >  1:3217636/8‑100 (MQ=21)
                                            ctAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAtt                                            >  1:2620836/2‑100 (MQ=21)
                                             ctagaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGGGGATTACTTAAAGAACAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGGAAttt                                           >  1:5350119/5‑100 (MQ=21)
                                              aGAAAAAGTTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:1140344/1‑99 (MQ=21)
                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:3611780/1‑99 (MQ=21)
                                              aGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTAGATTGCTTAAACAAAAGCATAAAATTTTAAATTTAGATAATGTTGAAGATGGTAAttta                                          >  1:1519091/1‑99 (MQ=21)
                                               agaagaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:318026/6‑100 (MQ=25)
                                               agaaaaGGTTTGTGTCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5714197/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGTTTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1192773/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGGAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1412427/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1801710/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2410589/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2028612/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4497811/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4586818/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:465813/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:472271/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:483333/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1932627/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4994327/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:54562/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5902834/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1754782/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5291275/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5765082/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4351406/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5430619/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:6071149/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:5505805/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3451401/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1387180/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2526399/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2334802/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4221578/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3757897/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2221924/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:100123/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGGGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:3731660/3‑100 (MQ=25)
                                               agaaaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGCGGATTACTTAAAGAAAAGCATAAACTTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:1328355/3‑100 (MQ=21)
                                               agaaaaGGTTTGTATCAAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:4527225/3‑100 (MQ=25)
                                               aaaaagGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGa                                         >  1:2957154/2‑100 (MQ=25)
                                                 caaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAACTTTAAACTTAGACGATGCCGAAGATGGTAATTTGAtt                                       >  1:735224/2‑100 (MQ=25)
                                                  aaaGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTa                                      >  1:307452/1‑100 (MQ=25)
                                                   aaGGTTTGTATCGAAAACGTATGTTAAGCTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTaa                                     >  1:4971090/1‑100 (MQ=25)
                                                      acttGTCTCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAGAGCATAAAATTTTAAACTTAGTCGATGCCGAAGATGGTAATTTGATTAAtac                                  >  1:2829192/3‑98 (MQ=25)
                                                       cttGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACa                                 >  1:3301031/2‑100 (MQ=25)
                                                         tGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAg                               >  1:4283474/1‑100 (MQ=37)
                                                                 taaCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAg                       >  1:198036/2‑100 (MQ=37)
                                                                  aaCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGa                      >  1:2184842/1‑100 (MQ=38)
                                                                  aaCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGa                      >  1:3386570/1‑100 (MQ=38)
                                                                   ccGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGAt                     >  1:4184402/2‑100 (MQ=255)
                                                                   aCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGAt                     <  1:188021/100‑1 (MQ=38)
                                                                   aCGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGAt                     >  1:2272097/1‑100 (MQ=38)
                                                                    tgTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATa                    >  1:5107982/2‑100 (MQ=255)
                                                                    cGTATGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATa                    >  1:297088/1‑100 (MQ=38)
                                                                       aTGTTAAGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAAGACAAGagatcggaagagca                 >  1:2061554/1‑90 (MQ=255)
                                                                            caGTTATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGGTTAATACAAGTGATGAAGATAAAacgaca            >  1:1674667/2‑100 (MQ=255)
                                                                               ttATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGAc         <  1:4086432/100‑1 (MQ=255)
                                                                                tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGACGATGGTAATTTGATTAATACAAGTGATGAAGATAAAAGGACAGaga        >  1:5938291/1‑98 (MQ=255)
                                                                                tATGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAATGCGACAGatc        >  1:4142086/1‑98 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAGAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAATCGACAGatcg       >  1:3880655/1‑97 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaactgacagatc       >  1:5146050/1‑91 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAacaacaaagat       >  1:1035788/1‑95 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGatcg       >  1:1066781/1‑97 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACga       >  1:24982/1‑100 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACAACAGatcg       >  1:5818753/1‑97 (MQ=255)
                                                                                 aTGGTGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACAACAGatcg       >  1:5235762/1‑97 (MQ=255)
                                                                                 aTGGGGGATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAaacggacagatc       >  1:5628513/1‑91 (MQ=255)
                                                                                   cgtggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAGCGACAGatcgga     >  1:1643209/2‑95 (MQ=255)
                                                                                      ggATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACgaagaag  >  1:4889856/1‑100 (MQ=255)
                                                                                      agATTACTTAAAGAAAAGCATAAAATTTTAAACTTAGACGATGCCGAAGATGGTAATTTGATTAATACAAGTGATGAAGATAAAACGACAGACgaagaag  >  1:3138682/2‑100 (MQ=255)
                                                                                                  |                                                                                       
TTTTAACTGATGATGATGAAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTAAACAAAAACATAAAATTTTAAACTTAGACGATGTCGAAGATGGTAATTTGATTAATGCAAGTGATGAAGATAAAACAACAGACGAAGAAG  >  NC_007790/746‑931

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