Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,563,518:1 (A)5→6 100% intergenic (+401/‑74) yffP → / → yffQ CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,563,5131.A97.9% 309.8 / ‑0.1 93intergenic (+396/‑79)yffP/yffQCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (50/41);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

TGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563417‑2563607
                                                                                                 |                                                                                               
tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                              >  1:2057040/1‑100 (MQ=255)
tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                              <  1:2371881/100‑1 (MQ=255)
       tgtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatgaa                                                                                      <  1:201671/100‑6 (MQ=255)
         tgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGATaaaaaataaaa                                                                                     >  1:93372/1‑93 (MQ=255)
          gATGCAGATCGCTTTTTCCCGTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                    >  2:1461690/1‑100 (MQ=255)
             gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTg                                                                                 >  2:2283395/1‑100 (MQ=255)
              cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGt                                                                                  >  1:192493/1‑98 (MQ=255)
               aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTggg                                                                               <  2:212469/100‑1 (MQ=255)
                   cGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGt                                                                           >  1:978319/1‑100 (MQ=255)
                    gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAg                                                                            >  2:1655502/1‑98 (MQ=255)
                    gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTg                                                                          <  1:1911728/100‑1 (MQ=255)
                    acTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTg                                                                          <  1:928745/99‑1 (MQ=255)
                      tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGc                                                                        <  1:952070/100‑1 (MQ=255)
                       ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCt                                                                       >  1:1841239/1‑100 (MQ=255)
                        tttCCCTTCAGAAATTTTTGTTTAATGGAATCAACGAAAGTCAGAAACCAACAACAAAATAAAGAAAAAAAATAAAAAATGAAACGGTg                                                                                 <  2:1221004/89‑1 (MQ=255)
                         ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                     >  1:1805236/1‑100 (MQ=255)
                          tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGATACTGTGGGAAGTGGCttt                                                                     <  2:209789/99‑1 (MQ=255)
                            ccTTCAGGCATTTTTTTTTTTTTTTTTTAACGAAAGACAGCTCCCAACAACAAAATAAAGATGGGGGGTAAAAAATGAATCGTTTGGAAGTCTCTTTACt                                                                  >  2:746557/1‑100 (MQ=255)
                            ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACt                                                                  >  1:1351874/1‑100 (MQ=255)
                             cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGCGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTg                                                                 >  1:413897/1‑100 (MQ=255)
                             cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTg                                                                 >  2:785033/1‑100 (MQ=255)
                               tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               >  1:2103496/1‑100 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               >  2:2345862/1‑99 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               <  1:2305192/99‑1 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  1:1641051/1‑100 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  2:4528/1‑100 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGAGAAGTGACTTTACTGGTg                                                              >  1:914503/1‑100 (MQ=255)
                                cAGGCATTTTTGTTTAATTCAATCAACAAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                                                               >  2:349272/1‑99 (MQ=255)
                                    cATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgg                                                          >  2:686343/1‑100 (MQ=255)
                                      tttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTa                                                        >  2:96145/1‑100 (MQ=255)
                                       ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAGATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       <  2:1678884/100‑1 (MQ=255)
                                       ttttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAAAGAAAAGGTGGGAAGTGGCTTTACAGGTGATGGTaa                                                       >  2:261232/1‑100 (MQ=255)
                                         ttGTTTAATTTAATCAACGAAAGGCAGTTAACAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGTTGGTAAGTGGCTTTACTGGTGATGGTAAt                                                      <  2:413085/99‑1 (MQ=255)
                                          tGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                                                    <  1:2323936/100‑1 (MQ=255)
                                          tGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                                                    <  2:492154/100‑1 (MQ=255)
                                           gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   >  1:1984480/1‑100 (MQ=255)
                                           gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAAATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   <  1:1649578/100‑1 (MQ=255)
                                            tttAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg                                                  >  1:1919900/1‑100 (MQ=255)
                                              tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGt                                                >  1:1061506/1‑100 (MQ=255)
                                               aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               <  1:1282502/100‑1 (MQ=255)
                                               aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               >  2:1678360/1‑100 (MQ=255)
                                               aaTTCAATCAACGAAAGGCAGTTCCCAAAAACAAAATAAAGATGAGAGGTAAAAAATGAAACGGTGGGAAGTGGCTCTACTGGTGATGGTAATg                                                     >  2:2347557/1‑94 (MQ=255)
                                                 ttCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg                                             <  2:716465/100‑1 (MQ=255)
                                                   caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 >  2:1604355/1‑94 (MQ=255)
                                                   caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                           >  1:782805/1‑100 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTGCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  1:1144643/100‑1 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  2:1357595/100‑1 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  2:2358230/100‑1 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  1:1596965/1‑100 (MQ=255)
                                                    aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          >  2:1619975/1‑100 (MQ=255)
                                                      tcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa                                        >  1:207202/1‑100 (MQ=255)
                                                       gatctaaaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  2:1861069/95‑1 (MQ=255)
                                                       caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       <  2:1863405/100‑1 (MQ=255)
                                                         aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     >  1:2254042/1‑100 (MQ=255)
                                                         aCGAAAGACAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     <  2:861255/100‑1 (MQ=255)
                                                          cGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGTGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGa                                    >  2:2387687/1‑100 (MQ=255)
                                                            aaaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtagatc                                  >  1:1861090/1‑95 (MQ=255)
                                                             aaGGCAGTTCCCAACAACAAAAGAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACt                                 <  2:119269/100‑1 (MQ=255)
                                                              aagCAGTTCCCAACAACAAAATAAAGATGAGAGGTAAAAAAGGAAACGGTGTGTAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGACCtg                                <  2:1446743/98‑1 (MQ=255)
                                                              aGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtg                                <  1:217843/100‑1 (MQ=255)
                                                                  aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGaa                            <  1:1059597/100‑1 (MQ=255)
                                                                   gTTCCCAACAATAAAATAAAGATTAGAGGTAAAAAATGAAACTGTGGGAAGTGGCTGTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  2:2307501/100‑1 (MQ=255)
                                                                   gTTCCCAACAACAAAATAAAGATGGGGGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  1:457047/100‑1 (MQ=255)
                                                                   gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  2:93601/100‑1 (MQ=255)
                                                                    ttttcaacaacaaAATAATGATTGGAGGTAAAATATGAAACGGTGGGAAGTTGCTTTATTGGTTATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  2:1162802/96‑1 (MQ=255)
                                                                    ttCCCAACAACAAAATAAAGATGGGAGTTAAAAAATGAATCGCTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  2:1686118/100‑1 (MQ=255)
                                                                    ttCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  1:2086054/100‑1 (MQ=255)
                                                                     tCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACTGTAGGAAGTGGCTTTACTGTAGATGGTAATGTTGGTTTGCTGTATTGAACTGTGGag                           >  2:1303750/1‑98 (MQ=255)
                                                                      cccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                        >  1:558008/1‑100 (MQ=255)
                                                                      cccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTAGCTTTACTGGTGATGGTAATGATGGTTTGCTGTATGGAACTGTGAagga                         >  2:1263762/1‑99 (MQ=255)
                                                                       ccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTTCTGTATTGAACTGTGAaggagg                       <  2:1018742/100‑1 (MQ=255)
                                                                       ccAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                       >  2:751450/1‑100 (MQ=255)
                                                                        caacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                      >  2:3009/1‑100 (MQ=255)
                                                                        caacaacaaAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGAc                     >  1:724540/1‑100 (MQ=255)
                                                                            aacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGcc                  >  2:351330/1‑100 (MQ=255)
                                                                             acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  1:244115/1‑100 (MQ=255)
                                                                             acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:752345/1‑100 (MQ=255)
                                                                             acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 <  2:1351860/100‑1 (MQ=255)
                                                                              caaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                >  2:1272678/1‑100 (MQ=255)
                                                                               aaaaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCAt                <  1:830750/99‑1 (MQ=255)
                                                                               aaaaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATg               >  1:1146344/1‑100 (MQ=255)
                                                                               aaaaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATg               <  1:179343/100‑1 (MQ=255)
                                                                                aaaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAGGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:1422963/1‑97 (MQ=255)
                                                                                 aaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             <  1:1347660/100‑1 (MQ=255)
                                                                                 aaTAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGaa             >  2:880682/1‑100 (MQ=255)
                                                                                     aaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTAATGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         >  2:79772/1‑100 (MQ=255)
                                                                                      aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        >  2:2087094/1‑100 (MQ=255)
                                                                                      aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        >  2:464187/1‑100 (MQ=255)
                                                                                       gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       <  2:2164470/100‑1 (MQ=255)
                                                                                       gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       >  1:450463/1‑100 (MQ=255)
                                                                                       gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACGTGTGATGGTAATGTTTGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        <  2:2392155/99‑1 (MQ=255)
                                                                                        aTGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa      <  2:1379298/100‑1 (MQ=255)
                                                                                        aTGCGAGGTAAAAAATGAAATGGGGGCAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTGTTa      <  2:2069564/100‑1 (MQ=255)
                                                                                          gggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    <  1:682781/100‑1 (MQ=255)
                                                                                           ggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATc   >  1:1174941/1‑100 (MQ=255)
                                                                                           ggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATc   <  2:978308/100‑1 (MQ=255)
                                                                                            gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  2:176248/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563417‑2563607

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