Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 145,268 G→A A118T (GCG→ACG)  SAUSA300_RS00650 → siderophore biosynthesis protein SbnI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793145,2680GA99.7% 930.9 / ‑5.6 296A118T (GCG→ACG) SAUSA300_RS00650siderophore biosynthesis protein SbnI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (290/5);  total (291/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

Only 100 of 298 total aligned reads displayed.

TTAGTACATGGCAACATAAAGTTCCATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAGCGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGAAAAGGTTGT  >  NC_007793/145173‑145362
                                                                                               |                                                                                              
         atacgcATAAAGTTCCATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTAc                                                                                   >  1:4541363/6‑100 (MQ=255)
          gCAACATAAAGTTCCATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACa                                                                                  >  1:5663458/1‑100 (MQ=255)
                       ccATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTg                                                                     >  1:3040207/1‑100 (MQ=255)
                       ccATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTg                                                                     >  1:3353696/1‑100 (MQ=255)
                       ccATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTg                                                                     >  1:3482507/1‑100 (MQ=255)
                       ccATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTg                                                                     >  1:5967851/1‑100 (MQ=255)
                       ccATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTg                                                                     >  1:5590966/1‑100 (MQ=255)
                        taTTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGa                                                                    >  1:5867324/2‑100 (MQ=255)
                         acttGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGAt                                                                   >  1:4734434/3‑100 (MQ=255)
                           ttGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATAc                                                                 >  1:5323969/1‑100 (MQ=255)
                           ttGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATAc                                                                 >  1:1004859/1‑100 (MQ=255)
                           ttGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATAc                                                                 >  1:8490/1‑100 (MQ=255)
                           ttGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATAc                                                                 >  1:516351/1‑100 (MQ=255)
                           ctGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATAc                                                                 >  1:6031657/2‑100 (MQ=255)
                            tGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACa                                                                >  1:2135277/1‑100 (MQ=255)
                            tGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACa                                                                >  1:1582765/1‑100 (MQ=255)
                             ggTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAg                                                               >  1:5439768/1‑100 (MQ=255)
                             ggTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAg                                                               >  1:1971242/1‑100 (MQ=255)
                                 tgtgGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACa                                                           >  1:4157710/1‑100 (MQ=255)
                                  gtggtgGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:4651477/1‑100 (MQ=255)
                                  gtggtgGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:4603886/1‑100 (MQ=255)
                                  gtggtgGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:4071492/1‑100 (MQ=255)
                                  gtggtgGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:3901501/1‑100 (MQ=255)
                                  gtggtgGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:3247315/1‑100 (MQ=255)
                                  gtggtgGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:1493810/1‑100 (MQ=255)
                                  atggtggGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACaa                                                          >  1:988519/2‑100 (MQ=255)
                                    ggtggGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAata                                                        >  1:5691189/1‑100 (MQ=255)
                                     ctggGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAatat                                                       >  1:3760656/2‑100 (MQ=255)
                                       gggAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATAttt                                                     >  1:4334207/1‑100 (MQ=255)
                                        ggAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTg                                                    >  1:2839948/1‑100 (MQ=255)
                                                 acaacaAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACaaaa                                           >  1:89112/1‑100 (MQ=255)
                                                 acaacaAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACaaaa                                           >  1:3341552/1‑100 (MQ=255)
                                                    acaaGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGAt                                        >  1:2869543/1‑100 (MQ=255)
                                                             acccgtGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGaa                               >  1:5997655/6‑100 (MQ=255)
                                                               ccTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGc                             >  1:4022126/2‑100 (MQ=255)
                                                               ccTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGAGGCATGGTGATACAGAACAATCTTTGCATACAAAGGATATAgatcggaa                             >  1:539742/2‑93 (MQ=255)
                                                                cTTGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGca                            >  1:1359194/1‑100 (MQ=255)
                                                                agtGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGca                            >  1:1436419/3‑100 (MQ=255)
                                                                 atGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGcac                           >  1:3824853/2‑100 (MQ=255)
                                                                 aggCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGcac                           >  1:1715569/3‑100 (MQ=255)
                                                                  tGCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGcaca                          >  1:2949974/1‑100 (MQ=255)
                                                                   gCCATGGACTACTTAGGCAAGACAAGAAACGCCATTTATTACAATGGTTCATGGGGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:1198307/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:2177076/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:4865279/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:5106776/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:4141591/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:1179134/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:146444/1‑100 (MQ=255)
                                                                   gCCATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAt                         >  1:2308441/1‑100 (MQ=255)
                                                                    ccATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAtt                        >  1:2066295/1‑100 (MQ=255)
                                                                     cgtGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGaaagatcgga                       >  1:5223759/3‑92 (MQ=255)
                                                                     cATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAttt                       >  1:4182765/1‑100 (MQ=255)
                                                                     cATGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAttt                       >  1:2601202/1‑100 (MQ=255)
                                                                      gtGGACTACTGAGACAAGACAAGAAGCGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAtttt                      >  1:264406/2‑100 (MQ=255)
                                                                      ctGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAtttt                      >  1:760537/2‑100 (MQ=255)
                                                                      aTGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACAtttt                      >  1:1621409/1‑100 (MQ=255)
                                                                       tGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTc                     >  1:5260018/1‑100 (MQ=255)
                                                                       tGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTc                     >  1:4755658/1‑100 (MQ=255)
                                                                       tGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTc                     >  1:711689/1‑100 (MQ=255)
                                                                       tGGACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTc                     >  1:1141696/1‑100 (MQ=255)
                                                                        tgACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:5887988/2‑100 (MQ=255)
                                                                        gtactactGAGACAAGACAAGAAACGCCATTTATTACCATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:642488/3‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAGACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:4701006/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAgatcggaagag                    >  1:93652/1‑90 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:2436700/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:2359801/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:330709/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:3664623/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:1518977/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:5993622/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:4302631/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:3091071/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:5037611/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:4446006/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:4478519/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:1911359/1‑100 (MQ=255)
                                                                        ggACTACTGAGACAAGACAAGAAACGCCATTTACTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCa                    >  1:5746674/1‑100 (MQ=255)
                                                                         gACTACTGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCaa                   >  1:4379500/1‑100 (MQ=255)
                                                                             actGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTAt               >  1:5361831/1‑100 (MQ=255)
                                                                             actGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTAt               >  1:5493377/1‑100 (MQ=255)
                                                                             actGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTAt               >  1:957037/1‑100 (MQ=255)
                                                                             actGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGGGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTAt               >  1:3611611/1‑100 (MQ=255)
                                                                              ctGAGACAAGACAAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATg              <  1:4508523/100‑1 (MQ=255)
                                                                                gagaCAAGACAAGAAACGCCATTTATTACAATGTGTTATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATggg            >  1:1277218/1‑100 (MQ=255)
                                                                                accacaagacaagAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATggg            >  1:746010/4‑100 (MQ=255)
                                                                                 ctgccagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGa           >  1:5567394/6‑100 (MQ=255)
                                                                                 agacaagacaAGAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATCGGa           >  1:5984290/1‑100 (MQ=255)
                                                                                  gtgcagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGaa          >  1:5045777/5‑100 (MQ=255)
                                                                                  cacgagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGaa          >  1:4043220/5‑100 (MQ=255)
                                                                                   tgcagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGGATgagatc         >  1:5296824/4‑98 (MQ=255)
                                                                                   ccaagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGaaa         >  1:2046424/2‑100 (MQ=255)
                                                                                   acaagacaagAAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATgggcca         >  1:312298/1‑97 (MQ=255)
                                                                                    gcagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGaaga        >  1:1032802/3‑98 (MQ=255)
                                                                                    ccagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGaaaa        >  1:177976/3‑100 (MQ=255)
                                                                                    caagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGaaaa        >  1:3859437/1‑100 (MQ=255)
                                                                                     tagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGAAAAg       >  1:3542526/2‑100 (MQ=255)
                                                                                     cagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATATAGGCGCAGCACATTTTCAAGTACGGGaaaga       >  1:5612916/2‑98 (MQ=255)
                                                                                     aagacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGagatc       >  1:1931325/1‑96 (MQ=255)
                                                                                      agacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTGCGGGaagatc      >  1:456850/1‑96 (MQ=255)
                                                                                      agacaagaAACGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATggggagatc      >  1:270618/1‑94 (MQ=255)
                                                                                               |                                                                                              
TTAGTACATGGCAACATAAAGTTCCATTTGGTGTGTGGTGGGAAACGTTACAACAAGAACATCGCTTGCCATGGACTACTGAGACAAGACAAGAAGCGCCATTTATTACAATGTGTCATGGTGATACAGAACAATATTTGTATACAAAAGATTTAGGCGAAGCACATTTTCAAGTATGGGAAAAGGTTGT  >  NC_007793/145173‑145362

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑