Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,418 A→G 100% R42R (CGT→CGC ykfF ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,4180AG97.9% 313.1 / ‑2.7 94R42R (CGT→CGCykfFCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (54/38);  total (55/39)
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.46e-01

CTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGC  >  NC_000913/265321‑265514
                                                                                                 |                                                                                                
cTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATcgcg                                                                                                <  1:2175526/100‑1 (MQ=255)
 tGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGa                                                                                               <  1:527333/100‑1 (MQ=255)
  gTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGaa                                                                                              >  1:2023843/1‑100 (MQ=255)
    aCGGATGCCTGAGGTGCGGATATAGCGGTTACGACCTTCGCCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAacn                                                                                            <  2:353005/100‑2 (MQ=25)
     cGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAacca                                                                                           <  1:227295/100‑1 (MQ=255)
      ggATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAaccac                                                                                          >  2:111750/1‑100 (MQ=255)
      ggATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCAATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAaccac                                                                                          >  1:1225676/1‑100 (MQ=255)
       gATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAaccacc                                                                                         <  2:7380/100‑1 (MQ=255)
          gCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGa                                                                                      >  2:1968513/1‑100 (MQ=255)
           ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                     >  1:2058321/1‑100 (MQ=255)
            cTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACg                                                                                    <  2:819310/100‑1 (MQ=255)
                ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                                <  1:539892/100‑1 (MQ=255)
                ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                                <  1:2164106/100‑1 (MQ=255)
                 gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGt                                                                               >  1:745557/1‑100 (MQ=255)
                   gCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGa                                                                             >  2:1449910/1‑100 (MQ=255)
                   gCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGa                                                                             >  1:774548/1‑100 (MQ=255)
                     ggATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACt                                                                           >  2:1068478/1‑100 (MQ=255)
                         ataGCGGTTAAGACCTTCACCGGCATCGGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGACCCCCCAGACGGAAATGACTGcc                                                                        >  2:2125347/1‑99 (MQ=25)
                         ataGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGccc                                                                       >  1:1756757/1‑100 (MQ=255)
                          taGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCt                                                                      >  1:1373243/1‑100 (MQ=255)
                          taGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCt                                                                      <  1:927497/100‑1 (MQ=255)
                          taGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCt                                                                      >  1:1768471/1‑100 (MQ=255)
                          taGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCt                                                                      >  1:2430896/1‑100 (MQ=255)
                           aGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTGTCGCGAACCACCAGACGGAAGTGACTGcc                                                                        >  2:514304/1‑97 (MQ=255)
                           aGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTg                                                                     <  2:1660104/100‑1 (MQ=255)
                            gCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCt                                                                      <  1:927406/98‑1 (MQ=255)
                            gCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTgg                                                                    <  2:1862633/100‑1 (MQ=255)
                             cGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgt                                                                   <  1:2407590/100‑1 (MQ=255)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCCCCAGACGAAAGTGACTGCCCTGgtcg                                                                 >  2:1657791/1‑100 (MQ=255)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGCCGGAAGTGACTGCCCTGgtcg                                                                 >  2:19996/1‑100 (MQ=255)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                 <  1:2285230/100‑1 (MQ=255)
                               gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                 <  2:378704/100‑1 (MQ=255)
                                 tAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATACACCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTACCCTGgccg                                                                 >  2:271326/1‑98 (MQ=25)
                                  aaGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCt                                                              >  1:2349160/1‑100 (MQ=255)
                                   aGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCtt                                                             >  1:1373203/1‑100 (MQ=255)
                                   aGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCtt                                                             >  2:1903787/1‑100 (MQ=255)
                                   aGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCtt                                                             >  2:1874127/1‑100 (MQ=255)
                                    gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTc                                                            <  1:2250825/100‑1 (MQ=255)
                                    gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTc                                                            >  1:2184319/1‑100 (MQ=255)
                                      ccTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCg                                                           >  2:1023000/1‑99 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                         <  1:1165239/100‑1 (MQ=255)
                                       cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                         <  1:97534/100‑1 (MQ=255)
                                        ttCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgt                                                        <  1:692/100‑1 (MQ=255)
                                        ttCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgt                                                        <  2:2153456/100‑1 (MQ=255)
                                         tCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtg                                                       <  2:666140/100‑1 (MQ=255)
                                          cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                      <  2:481909/100‑1 (MQ=255)
                                           aCCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGtgtg                                                     <  1:243211/100‑1 (MQ=255)
                                            ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGa                                                    >  1:410162/1‑100 (MQ=255)
                                            ccGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGa                                                    <  1:344858/100‑1 (MQ=255)
                                                cATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCCCCAGACGGAAGTGACTGCCCTGGTCGTCTTCCAGCGTGacct                                                 >  2:777909/1‑96 (MQ=25)
                                                  tCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGc                                              >  1:826644/1‑100 (MQ=255)
                                                   ccGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCt                                             <  2:2042864/100‑1 (MQ=255)
                                                    cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCt                                             <  2:2042832/99‑1 (MQ=255)
                                                    cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                            >  1:1865047/1‑100 (MQ=255)
                                                     ggCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGCCTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGt                                           >  2:554309/1‑100 (MQ=255)
                                                      gCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                          <  1:1067195/100‑1 (MQ=255)
                                                      gCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                          >  2:1405245/1‑100 (MQ=255)
                                                      gCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                          >  1:1272619/1‑100 (MQ=255)
                                                      gCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                          >  1:1615280/1‑100 (MQ=255)
                                                       cTCAAAGTTCCATGCCCGCCAGACCCTCCGGCCTTCAGTATCGCGAACCACCAGACCGACGTGACTGCCCTGGTCGTCCCCGAGTGTGATATTGCTGTAc                                         >  2:2154749/1‑100 (MQ=37)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGCACCACCAGACGGAAGTGACTGCCCTGGTCGCCTTCGAGTCTGACATCGCt                                             >  2:361205/1‑96 (MQ=25)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTAc                                         <  2:745550/100‑1 (MQ=255)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTAc                                         <  1:1574362/100‑1 (MQ=255)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTAc                                         >  1:1745819/1‑100 (MQ=255)
                                                       cTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGCGACTGCCCGGGTCGCCTTCGAGTGTGATACTGCTGTAc                                         >  2:945032/1‑100 (MQ=37)
                                                         cAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGt                                       >  2:215886/1‑100 (MQ=255)
                                                          aaaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTg                                      >  2:1229266/1‑100 (MQ=255)
                                                          aaaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTg                                      >  1:1552941/1‑100 (MQ=255)
                                                           cttccgatctcccGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                     <  2:1145614/90‑1 (MQ=38)
                                                           aaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                     <  1:645959/100‑1 (MQ=255)
                                                             gtgcgctGCCCGCCAGCCCGTCCGGCCTTCAGTGTCGCGAATCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                     <  2:2018772/92‑1 (MQ=255)
                                                              ttCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                  <  1:1712939/100‑1 (MQ=255)
                                                               tCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGt                                 <  2:1372704/100‑1 (MQ=255)
                                                                    gCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGCGTGTGATATTGCTGTACGTGGTAGTGACcg                            >  2:353345/1‑100 (MQ=255)
                                                                    gCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcg                            >  1:1594546/1‑100 (MQ=255)
                                                                     cccGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGGGTGATATTGCTGTACGTggagatcggaag                           >  1:1145622/1‑90 (MQ=38)
                                                                      ccGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgct                          >  2:2053068/1‑100 (MQ=255)
                                                                       cGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgctt                         >  2:810099/1‑100 (MQ=255)
                                                                       cGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgctt                         >  1:2052661/1‑100 (MQ=255)
                                                                          cAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcg                      >  1:2159512/1‑100 (MQ=255)
                                                                          cAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcg                      <  1:1657803/100‑1 (MQ=255)
                                                                           aGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgc                     <  1:2152012/100‑1 (MQ=255)
                                                                            gACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                    <  2:2419964/100‑1 (MQ=255)
                                                                            gACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                    <  1:1222456/100‑1 (MQ=255)
                                                                            gACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGACTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                    >  2:885354/1‑100 (MQ=255)
                                                                              ccATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgcttg                  >  2:235840/1‑100 (MQ=255)
                                                                               cATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGt                 >  1:1045032/1‑100 (MQ=255)
                                                                                aTCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtc                >  1:2120976/1‑100 (MQ=255)
                                                                                aTCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATGGCTGTACGTGGTAGTGACCGCTTGCGCTTGtc                >  2:890574/1‑100 (MQ=255)
                                                                                 tCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGCTATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtct               >  2:2389034/1‑100 (MQ=255)
                                                                                         tCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGGAGTGACCGCTTGCGCTTGTCTCCGGGTGa       >  1:1113867/1‑100 (MQ=255)
                                                                                          cAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGaa      <  2:1031177/100‑1 (MQ=255)
                                                                                          cAGTATCGCGAAACACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGaa      >  1:1690552/1‑100 (MQ=255)
                                                                                            gTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAg    <  2:2023829/100‑1 (MQ=255)
                                                                                              aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc  >  1:2194469/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGC  >  NC_000913/265321‑265514

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