Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 586,465 T→A V304V (GTT→GTA)  SAUSA300_RS02820 → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793586,4650TA99.4% 1058.1 / ‑3.8 359V304V (GTT→GTA) SAUSA300_RS02820DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (354/2);  total (357/2)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 395 total aligned reads displayed.

TATAAAACAAACAAAAAATTACATTTAAAACATCGTTTATTTAATCAAAAATTAGCTGAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTTGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTAGC  >  NC_007793/586370‑586561
                                                                                               |                                                                                                
tataAAACAAACAAAAAATTACATTTAAAACATCGTTTATTTAATCAAAAATTAGCTGAGTCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAgaaa                                                                                              <  1:2157385/100‑2 (MQ=255)
                                                   ttAGCTGAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTg                                           >  1:833829/1‑100 (MQ=255)
                                                    tAGCTGAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGa                                          >  1:593860/1‑100 (MQ=255)
                                                     aGCTGAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGaa                                         >  1:2316427/1‑100 (MQ=255)
                                                        tGAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCa                                      >  1:1381708/1‑100 (MQ=255)
                                                        ctaGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCa                                      >  1:4485722/3‑100 (MQ=255)
                                                         taGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:3040776/2‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:2889895/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:770687/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:796970/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:5770083/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:1340430/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:1174283/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:5540662/1‑100 (MQ=255)
                                                         gAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:354019/1‑100 (MQ=255)
                                                         atgCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaa                                     >  1:3407174/3‑100 (MQ=255)
                                                          aGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaaa                                    >  1:5576073/1‑100 (MQ=255)
                                                          aGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCCTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCaaa                                    >  1:37829/1‑100 (MQ=255)
                                                           gCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAAt                                   >  1:1700225/1‑100 (MQ=255)
                                                                 tGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGTTGTACTTGAATCAAATGCAAAc                             >  1:1744101/1‑100 (MQ=255)
                                                                 tGTAAATACTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAAc                             >  1:5118742/1‑100 (MQ=255)
                                                                        aCTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAg                      >  1:2820695/1‑100 (MQ=255)
                                                                         ttGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgt                     >  1:3871143/2‑100 (MQ=255)
                                                                         cTGAAACTGGTGAAGTTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgt                     >  1:1490503/1‑100 (MQ=255)
                                                                         cTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgt                     >  1:1663679/1‑100 (MQ=255)
                                                                         cTGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgt                     >  1:4180869/1‑100 (MQ=255)
                                                                          tGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtg                    >  1:5285603/1‑100 (MQ=255)
                                                                          tGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtg                    >  1:2685767/1‑100 (MQ=255)
                                                                          tGAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGAACTTGAATCAAATGCAAACAGCGAAgtg                    >  1:3096545/1‑100 (MQ=255)
                                                                           taaaCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:5592240/2‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:3339166/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:3131318/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:1123986/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:4846204/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:926827/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:1074648/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:4110413/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:5238469/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:2933492/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:5019959/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:4407884/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:2594184/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:4590451/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:2037518/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:2467895/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:2404198/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgtgt                   >  1:4951738/1‑100 (MQ=255)
                                                                           gAAACTGGTGAAATTGTAGTAGAAGAAGGTACAATGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGTGAAgtgt                   >  1:3714/1‑100 (MQ=255)
                                                                            aaaCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGtt                  >  1:2223550/1‑100 (MQ=255)
                                                                             taCTGGTGAAATTGTAGTAGAAGAAGGAACAGTGCTTGATCGTCGGAAAATCGACGAAATCATGGATATACTTGAATCAAATGCAAACAGCGAAGAGttt                 >  1:1233403/2‑100 (MQ=255)
                                                                             gaCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGttt                 >  1:2727095/2‑100 (MQ=255)
                                                                             aaCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGttt                 >  1:887374/1‑100 (MQ=255)
                                                                              tatGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTg                >  1:5476569/3‑100 (MQ=255)
                                                                              ccTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAgatcgga                >  1:1579048/2‑94 (MQ=255)
                                                                              ccTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTg                >  1:4635035/2‑100 (MQ=255)
                                                                              aCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTg                >  1:4723624/1‑100 (MQ=255)
                                                                              aCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTg                >  1:5140313/1‑100 (MQ=255)
                                                                              aCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTg                >  1:3311388/1‑100 (MQ=255)
                                                                              aCTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTg                >  1:2915671/1‑100 (MQ=255)
                                                                               cTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGa               >  1:3559707/1‑100 (MQ=255)
                                                                               cTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGa               >  1:653407/1‑100 (MQ=255)
                                                                               cTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGa               >  1:1966935/1‑100 (MQ=255)
                                                                               cTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGa               >  1:1815230/1‑100 (MQ=255)
                                                                               cTGGTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGa               >  1:1559127/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCGAATGCAAACAGCGAAGTGTTTGAAt             >  1:1009542/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:5068080/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:5004022/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:57101/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:104778/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:4427716/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:4350250/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:4036887/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:1855635/1‑100 (MQ=255)
                                                                                 ggTGAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAAt             >  1:2492349/1‑100 (MQ=255)
                                                                                    gAAATTGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGc          >  1:1509758/1‑100 (MQ=255)
                                                                                       gctGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATg       >  1:2995660/3‑100 (MQ=255)
                                                                                        ttGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATgg      >  1:3062341/1‑100 (MQ=255)
                                                                                        ttGTAGTAGAAGAAGGTACAGTGCCTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATgg      >  1:3646425/1‑100 (MQ=255)
                                                                                        atGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATgg      >  1:4868236/2‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:728532/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:3223680/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:5432576/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:5310710/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:2364719/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:3905768/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:2091596/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:4550749/1‑100 (MQ=255)
                                                                                         tGTAGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGt     >  1:2603038/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAATGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:2555303/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTCCTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:5465183/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:4653491/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:1101347/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:5608149/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:5196489/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:3971252/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:4005107/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:4322190/1‑100 (MQ=255)
                                                                                          gtagtaGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTa    >  1:4817874/1‑100 (MQ=255)
                                                                                            aGTAGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTAGc  >  1:1268079/1‑100 (MQ=255)
                                                                                            aGTAGAAGAAGGTACAGTGCTTGATCGCCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGCAGTGTTTGAATTGCATGGTAGc  >  1:3697393/1‑100 (MQ=255)
                                                                                               |                                                                                                
TATAAAACAAACAAAAAATTACATTTAAAACATCGTTTATTTAATCAAAAATTAGCTGAGCCAATTGTAAATACTGAAACTGGTGAAATTGTAGTTGAAGAAGGTACAGTGCTTGATCGTCGTAAAATCGACGAAATCATGGATGTACTTGAATCAAATGCAAACAGCGAAGTGTTTGAATTGCATGGTAGC  >  NC_007793/586370‑586561

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