Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,867,269 T→A 16.3% I95F (ATC→TTC)  rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,2690TA16.3% 255.7 / 44.4 86I95F (ATC→TTC) rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (10/4);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAA  >  NC_000913/2867173‑2867367
                                                                                                |                                                                                                  
gCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatca                                                                                                  <  1:660046/99‑1 (MQ=255)
 ccccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatcat                                                                                                 >  2:1590847/1‑99 (MQ=255)
   ccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATcc                                                                                               >  1:2393142/1‑99 (MQ=255)
     aGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcgg                                                                                             <  2:989269/99‑1 (MQ=255)
     aGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcgg                                                                                             <  2:2367381/99‑1 (MQ=255)
      ggTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggc                                                                                            >  2:393294/1‑99 (MQ=255)
      ggTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCcggc                                                                                            >  1:1202282/1‑99 (MQ=255)
      ggTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCcggc                                                                                            >  1:913922/1‑99 (MQ=255)
        ttGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcgg                                                                                          >  1:1304606/1‑99 (MQ=255)
        ttGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggcgg                                                                                          >  1:403358/1‑99 (MQ=255)
         tGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACTAGACGCAAGTTACTCTCGATCATCCggcggc                                                                                         >  1:991759/1‑99 (MQ=255)
            ccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCgag                                                                                      >  1:103189/1‑99 (MQ=255)
            ccTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCgag                                                                                      <  2:1679727/99‑1 (MQ=255)
             ctctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCgaga                                                                                     >  1:2030925/1‑99 (MQ=255)
              tctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGagag                                                                                    >  1:1285067/1‑99 (MQ=255)
                ttCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGc                                                                                  >  2:682833/1‑99 (MQ=255)
                  cGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGa                                                                                <  1:2251956/99‑1 (MQ=255)
                    aTAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACa                                                                              >  1:1192691/1‑99 (MQ=255)
                    aTAAAGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACa                                                                              <  2:304569/99‑1 (MQ=255)
                       aGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtct                                                                           >  1:2152802/1‑99 (MQ=255)
                       aGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACAtct                                                                           >  1:737153/1‑99 (MQ=255)
                        ggTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtctc                                                                          <  2:1787741/99‑1 (MQ=255)
                        ggTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtctc                                                                          <  1:655837/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                        <  1:858368/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                        <  1:638866/99‑1 (MQ=255)
                          tCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCa                                                                        >  2:351757/1‑99 (MQ=255)
                            cAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACg                                                                      >  1:343079/1‑99 (MQ=255)
                             aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGc                                                                     <  1:457634/99‑1 (MQ=255)
                             aGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGc                                                                     >  1:410859/1‑99 (MQ=255)
                              gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCa                                                                    <  1:1726228/99‑1 (MQ=255)
                              gCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCa                                                                    >  1:1676305/1‑99 (MQ=255)
                               cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                   <  1:17956/99‑1 (MQ=255)
                               cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                   <  2:1507001/99‑1 (MQ=255)
                               cAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                                                                   <  2:1379156/99‑1 (MQ=255)
                                 aCGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTg                                                                 >  2:1212193/1‑99 (MQ=255)
                                      aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa                                                            <  1:2401865/99‑1 (MQ=255)
                                      aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa                                                            <  1:146345/99‑1 (MQ=255)
                                        aCCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcg                                                          <  2:414539/99‑1 (MQ=255)
                                         ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc                                                         <  2:369413/99‑1 (MQ=255)
                                          cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgcg                                                        >  1:2228785/1‑99 (MQ=255)
                                           aCGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACgcgc                                                       <  2:1996173/99‑1 (MQ=255)
                                            cGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCa                                                      >  2:6671/1‑99 (MQ=255)
                                              aTTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaa                                                    >  2:2261097/1‑99 (MQ=255)
                                               ttGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaaa                                                   >  1:2530046/1‑99 (MQ=255)
                                                tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAAt                                                  >  1:723089/1‑99 (MQ=255)
                                                 gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATa                                                 <  2:1979262/99‑1 (MQ=255)
                                                 gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATa                                                 >  2:260333/1‑99 (MQ=255)
                                                    aTAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAAc                                              <  1:1929939/99‑1 (MQ=255)
                                                     tAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAAct                                             >  1:919098/1‑99 (MQ=255)
                                                      aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActt                                            >  2:2237200/1‑99 (MQ=255)
                                                       aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtc                                           >  1:2158009/1‑99 (MQ=255)
                                                       aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtc                                           >  1:2306363/1‑99 (MQ=255)
                                                        cggcggGCAATTTTTACCACCGGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtct                                          <  2:1469796/99‑1 (MQ=255)
                                                        cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtct                                          <  1:2395682/99‑1 (MQ=255)
                                                        cggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtct                                          <  2:2084920/99‑1 (MQ=255)
                                                         ggcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtctt                                         <  2:46229/99‑1 (MQ=255)
                                                          gcggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttc                                        <  1:758443/99‑1 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                      <  1:1396024/99‑1 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                      >  2:1993141/1‑99 (MQ=255)
                                                            gggCAATTTTTACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtcttctt                                      >  1:729257/1‑99 (MQ=255)
                                                                  tttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCg                                >  2:2078204/1‑99 (MQ=255)
                                                                   ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGt                               <  1:554219/99‑1 (MQ=255)
                                                                    tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGtt                              >  2:2030009/1‑99 (MQ=255)
                                                                      tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTaa                            >  1:1782168/1‑99 (MQ=255)
                                                                      tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTaa                            >  2:408387/1‑99 (MQ=255)
                                                                      tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTaa                            >  1:2168029/1‑99 (MQ=255)
                                                                      tACCACCAGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTaa                            >  1:1010921/1‑99 (MQ=255)
                                                                       accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAAc                           >  2:217274/1‑99 (MQ=255)
                                                                         caccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAg                         <  2:1202282/99‑1 (MQ=255)
                                                                          accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgt                        >  2:1686392/1‑99 (MQ=255)
                                                                          accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgt                        <  1:149820/99‑1 (MQ=255)
                                                                          accaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgt                        <  2:1999568/99‑1 (MQ=255)
                                                                           ccaGACGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtg                       <  2:1612257/99‑1 (MQ=255)
                                                                             aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAgtggt                     <  1:473595/99‑1 (MQ=255)
                                                                               aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGa                   <  2:1522512/99‑1 (MQ=255)
                                                                               aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGa                   <  1:802875/99‑1 (MQ=255)
                                                                                cGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                  >  1:2518864/1‑99 (MQ=255)
                                                                                cGCAAGTTACTCTCGAACATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGaa                  >  1:1647412/1‑99 (MQ=255)
                                                                                   aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATaa               <  1:53995/99‑1 (MQ=255)
                                                                                   aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATaa               >  2:839808/1‑99 (MQ=255)
                                                                                    aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAc              >  2:1416414/1‑99 (MQ=255)
                                                                                    aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAAc              <  1:2341008/99‑1 (MQ=255)
                                                                                      ttACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCa            >  1:654755/1‑99 (MQ=255)
                                                                                       tACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCaa           >  2:871645/1‑99 (MQ=255)
                                                                                        aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAt          <  2:680240/99‑1 (MQ=255)
                                                                                         ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtc         >  1:418152/1‑99 (MQ=255)
                                                                                          tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAtct        >  2:419185/1‑99 (MQ=255)
                                                                                            tcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCa      <  2:2335458/99‑1 (MQ=255)
                                                                                               atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa   >  1:1690828/1‑99 (MQ=255)
                                                                                               atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa   >  2:124219/1‑99 (MQ=255)
                                                                                               atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCa   >  2:1874280/1‑99 (MQ=255)
                                                                                                tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCaa  >  1:1693209/1‑99 (MQ=255)
                                                                                                |                                                                                                  
GCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAA  >  NC_000913/2867173‑2867367

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