Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007790 744 T→C D155D (GAT→GAC)  SAUSA300_RS14705 → replication protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077907440TC99.1% 280.7 / ‑3.4 106coding (465/717 nt)SAUSA300_RS14705replication protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (86/19);  minor base . (1/0);  total (87/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

Only 100 of 106 total aligned reads displayed.

AAAAGGCGATAAAGACATTGAATCGGCAATCAAAGAGACCTCAAAATATTCGGTTAAATCATCTGATTTTTTAACTGATGATGATG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTA  >  NC_007790/678‑841
                                                                  |                                                                                                  
caaaGGAGATAAAGATATTGCGTCTGCTATCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGaa                                                                  >  1:104/2‑100 (MQ=255)
      agATAAAGATATTGCGTCTGCTATCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAAttatc                                                            <  1:4542502/99‑4 (MQ=255)
                     tataacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaag                                             >  1:221725/6‑99 (MQ=255)
                      ataacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaagg                                            >  1:890634/5‑98 (MQ=255)
                       taacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggt                                           >  1:2750975/4‑97 (MQ=255)
                       gaccAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggt                                           >  1:3737961/4‑97 (MQ=255)
                       gaacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggt                                           >  1:5576103/4‑97 (MQ=255)
                       ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                           >  1:1713272/6‑100 (MQ=255)
                       ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                           >  1:3210335/6‑100 (MQ=255)
                       ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                           >  1:441617/6‑100 (MQ=255)
                       ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                           >  1:6018524/6‑100 (MQ=255)
                       ctgctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAgg                                           >  1:2943434/6‑100 (MQ=255)
                        aacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggtt                                          >  1:435161/3‑96 (MQ=255)
                        aacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggtt                                          >  1:5702444/3‑96 (MQ=255)
                        aacAATCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGaaaaggtt                                          >  1:2484140/3‑96 (MQ=255)
                         gctaTCAAGGTGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:5533489/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATGAAAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGt                                          >  1:1362130/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:2945307/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:629766/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:3267193/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:3393980/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:3945046/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:779373/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:4017731/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:426309/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:4370099/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:5213278/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:5367808/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:5603125/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:2027631/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:1595894/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGtt                                         >  1:1148623/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACAGATAGCTATG‑AAAAAAATCAAGAAATTATCAGTGTCCTAGAAAAAGGtt                                         >  1:5883051/4‑100 (MQ=255)
                         gctaTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGA‑CTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGttt                                        >  1:2773697/4‑100 (MQ=255)
                         ctgcccAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGttt                                        >  1:2356951/6‑100 (MQ=255)
                         acAATCAAAGAGACATCAAAATATTCCTTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTTAGTGATCTAG‑AAAAGGttt                                        <  1:2712678/99‑1 (MQ=255)
                          gggggagggagaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGttt                                        >  1:6015640/9‑100 (MQ=255)
                            ctCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGGGATCTAGAAAAAGGTTTGt                                      >  1:776831/2‑100 (MQ=255)
                            cgcAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTa                                     >  1:5761228/3‑100 (MQ=255)
                            cccAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTa                                     >  1:5887239/3‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATGAAAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTg                                       >  1:2803011/1‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGt                                      >  1:1091563/1‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGt                                      >  1:5415962/1‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGt                                      >  1:2336604/1‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGt                                      >  1:681172/1‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGt                                      >  1:3072070/1‑100 (MQ=255)
                            aTCAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGt                                      >  1:837156/1‑100 (MQ=255)
                            aTCAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTCTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTa                                     >  1:1087240/1‑100 (MQ=255)
                             ccAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAGAAAAAgatgagat                                     >  1:771387/2‑93 (MQ=255)
                             acagaaggaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTAt                                    >  1:4196737/8‑100 (MQ=255)
                              cAAGGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTAt                                    >  1:2825272/1‑100 (MQ=255)
                              cAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATc                                   <  1:5588955/100‑1 (MQ=255)
                              cAAAGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCCAGaaaagatcggaaga                                    >  1:5457938/1‑90 (MQ=255)
                              cAAAGAGACATCAAAATATTCCGTTAAATCATCTGACCTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGCTTGTATc                                   >  1:2826308/1‑100 (MQ=255)
                               ccagagACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCg                                  >  1:3232872/3‑100 (MQ=255)
                               atatgggggtCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCg                                  >  1:2982748/10‑100 (MQ=255)
                                aaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCg                                  >  1:1048700/1‑100 (MQ=255)
                                aaGAGACATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACA‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGa                                 >  1:2701390/1‑100 (MQ=255)
                                  gagaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGaaa                               >  1:2530701/1‑100 (MQ=255)
                                   tgaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGaaaa                              >  1:4534128/2‑100 (MQ=255)
                                   agaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGaaa                               >  1:489545/1‑100 (MQ=255)
                                   aaacatCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGaaaa                              >  1:2843240/6‑100 (MQ=255)
                                    gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGaaaa                              >  1:741493/1‑100 (MQ=255)
                                    gaCATCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGaaaa                              >  1:889749/1‑100 (MQ=255)
                                    cacatCAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAAc                             >  1:4276682/5‑100 (MQ=255)
                                         cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑ATAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATg                        >  1:4093208/1‑100 (MQ=255)
                                         cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTaat                        <  1:294743/100‑3 (MQ=255)
                                         cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATg                        >  1:3142804/1‑100 (MQ=255)
                                         cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATg                        >  1:1520127/1‑100 (MQ=255)
                                         cAAAATATTCCGTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTAt                         <  1:4273866/100‑1 (MQ=255)
                                           aaaTATTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGtt                      <  1:1312032/100‑1 (MQ=37)
                                           aaaTATTCCGTTAAAGCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGtt                      <  1:1792724/100‑1 (MQ=37)
                                            caggccccgTTAAATCATCTGACTTTTTAACTGATGACGACA‑AAAAAAATCAAGAAATACTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGCTa                     >  1:2393793/9‑100 (MQ=255)
                                              ttctccgTTAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTaa                    >  1:1355200/7‑100 (MQ=25)
                                               tgcccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAATTTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGt                  >  1:4123715/6‑100 (MQ=25)
                                               atTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGt                  >  1:1238017/1‑100 (MQ=25)
                                               atTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGt                  >  1:3731923/1‑100 (MQ=25)
                                               atTCCGTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGt                  >  1:3561498/1‑100 (MQ=25)
                                                 tccgTTAAATCATCTGACTTTTTAACTGATGACGACGAAAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGtt                 <  1:5900169/97‑1 (MQ=25)
                                                 tccgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTa                >  1:1145648/4‑100 (MQ=25)
                                                   cgTTAAATCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTAta              <  1:5316491/99‑2 (MQ=25)
                                                    ggggcctcatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtgg             >  1:4377072/7‑100 (MQ=25)
                                                      ttccgatcttGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:3754809/91‑1 (MQ=21)
                                                      tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:2486799/100‑1 (MQ=25)
                                                      tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:5911955/100‑1 (MQ=25)
                                                      tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:730791/100‑1 (MQ=25)
                                                      tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:107290/100‑1 (MQ=25)
                                                      tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:5398156/100‑1 (MQ=25)
                                                      tAAATCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTAACAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggt            <  1:3428258/100‑1 (MQ=25)
                                                       aggtcatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtgg          <  1:5591435/97‑1 (MQ=25)
                                                       aaaTCATCTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtgg          >  1:5751702/1‑100 (MQ=25)
                                                        cttcctcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGa         >  1:4084710/6‑100 (MQ=25)
                                                        aaTCATCTGACTTTTTAACTGATAGCGATG‑AAAAAAATCAAGAAATTATCAGTGATCTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTAtggtag          <  1:4184074/100‑3 (MQ=25)
                                                         atcatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGAt        >  1:4732997/1‑100 (MQ=25)
                                                          ccatcTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGAtt       >  1:2285589/2‑100 (MQ=25)
                                                              cTGACTTTTTAACTGATGACGACG‑AGAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACtt   >  1:5506430/1‑100 (MQ=25)
                                                              cTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACtt   >  1:3351071/1‑100 (MQ=25)
                                                              cTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACtt   >  1:1959118/1‑100 (MQ=25)
                                                              cTGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACtt   >  1:4902097/1‑100 (MQ=25)
                                                               tGACTTTTTAACTGATGACGACG‑AAAAAAATCAAGAAATTTTAAGTGATCTAG‑AAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTACTTa  >  1:2735786/1‑100 (MQ=25)
                                                                  |                                                                                                  
AAAAGGCGATAAAGACATTGAATCGGCAATCAAAGAGACCTCAAAATATTCGGTTAAATCATCTGATTTTTTAACTGATGATGATG‑AAAAAAATCAAGAAATTGTAAGTGATTTAGAAAAAGGTTTGTATCGAAAACGTATGTTAAGTTATGGTGGATTGCTTA  >  NC_007790/678‑841

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: ‑