Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,563,518 (A)5→6 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.5% 287.8 / 0.3 81intergenic (+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 (43/36);  total (44/37)
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.98e-01

AGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATC  >  NC_000913/2563414‑2563606
                                                                                                    |                                                                                              
aGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGt                                                                                                 <  1:2075682/100‑1 (MQ=255)
 ggTGAATTTTGTGATGCAGATCGCTTTTTTTTTTTTTTCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a                                                                                               >  2:144617/1‑100 (MQ=255)
 ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a                                                                                               >  1:511604/1‑100 (MQ=255)
 ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a                                                                                               <  2:147452/100‑1 (MQ=255)
 ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑a                                                                                               >  1:1279203/1‑100 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCTATTCAGGCATTTTTGTTTATTTCAATAAACGAACGGCCGTTCACAACAACAAAATGAAGATGGGAGGa‑aa                                                                                              >  2:1044220/1‑99 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaa                                                                                             >  1:873045/1‑100 (MQ=255)
   tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaa                                                                                             <  2:2259107/100‑1 (MQ=255)
    gAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaa‑a                                                                                           >  1:754167/1‑100 (MQ=255)
       ttttGTGATGCAGATCGCTTTTTCCCTTCAGTCCTGTTTGTTTAATTCAATCAACGAAAGAAAGTTCCCAACAACAAAATAAAGATGGGAGGT‑gaa‑aa                                                                                          <  2:2001259/98‑1 (MQ=255)
        tttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑aaa‑aaatg                                                                                       <  1:717244/100‑4 (MQ=255)
               tGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGt                                                                                 >  2:1146690/1‑100 (MQ=255)
                 cAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCATTCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGt                                                                                 >  2:714965/1‑98 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGa                                                                             >  2:1130309/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGa                                                                             >  2:1130263/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGa                                                                             >  1:979862/1‑100 (MQ=255)
                   gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGa                                                                             <  2:511602/100‑1 (MQ=255)
                    aTCGATTTTTCCGTTGAGGCATTTTTGTTTAATTCAATCGACGTGAGGAAGTTTCTATCAACACACTCAAGATGGGAGGT‑AAAGATTGCAAAGGTGGGaa                                                                            >  2:1634976/1‑100 (MQ=255)
                     tCGCTTTTTCTCTTTATTCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAAAAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAg                                                                           >  2:11137/1‑100 (MQ=255)
                     tCGCTTTTTCCCTTCAGGCATTTTTTTTTTTTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAg                                                                           >  2:1174923/1‑100 (MQ=255)
                      cGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGt                                                                          >  1:1469469/1‑100 (MQ=255)
                       gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTg                                                                         <  1:1980515/100‑1 (MQ=255)
                       gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTg                                                                         >  2:115184/1‑100 (MQ=255)
                       gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTg                                                                         >  2:1987137/1‑100 (MQ=255)
                           tttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCtt                                                                     <  2:1680173/100‑1 (MQ=255)
                              cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTAc                                                                  <  2:1145550/100‑1 (MQ=255)
                               ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACt                                                                 <  2:624745/100‑1 (MQ=255)
                               ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACt                                                                 <  1:1444814/100‑1 (MQ=255)
                                 ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTgg                                                               <  1:1222666/100‑1 (MQ=255)
                                  tCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGt                                                              <  1:1341782/100‑1 (MQ=255)
                                   cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                             <  1:1146692/100‑1 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                           >  1:838461/1‑100 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                           >  1:1245269/1‑100 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                           <  1:208118/100‑1 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                           <  2:2259060/100‑1 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                           >  2:333569/1‑100 (MQ=255)
                                     ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AAGGAAACGGTGGGAAGTGGCTTTACTGGTGa                                                            >  2:758020/1‑99 (MQ=255)
                                      gCATTTTTGTTTCATTCAATTAACGAAAGGAAGTTCCCAACAACAAACTACTGATGGGAGGTAAAA‑AATGAAAAGGTGGGAAGTGGCTTTACTGGTGATg                                                          >  1:926214/1‑100 (MQ=255)
                                      gCATTTTTGTTTAATTCAATCAACGAAAGTCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTg                                                                >  1:926206/1‑94 (MQ=255)
                                      gCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATg                                                          <  2:754166/100‑1 (MQ=255)
                                            ttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                                                    <  1:1130265/100‑1 (MQ=255)
                                            ttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                                                    <  1:1130311/100‑1 (MQ=255)
                                            ttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                                                    <  1:758021/100‑1 (MQ=255)
                                              gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                  >  1:704597/1‑100 (MQ=255)
                                              gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                  >  2:1688268/1‑100 (MQ=255)
                                              gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                  <  2:30990/100‑1 (MQ=255)
                                                  aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                              <  1:714627/100‑1 (MQ=255)
                                                  aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                              <  2:22165/100‑1 (MQ=255)
                                                  aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                              >  2:698754/1‑100 (MQ=255)
                                                   aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtat                                             >  1:1211497/1‑98 (MQ=255)
                                                    ttCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg                                            <  1:508112/100‑1 (MQ=255)
                                                    ttCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg                                            >  2:43352/1‑100 (MQ=255)
                                                     tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                           <  1:714966/100‑1 (MQ=255)
                                                     tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                           <  1:1174926/100‑1 (MQ=255)
                                                      caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                          >  2:1388825/1‑100 (MQ=255)
                                                      caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                          >  1:890702/1‑100 (MQ=255)
                                                      caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                          >  1:314933/1‑100 (MQ=255)
                                                       aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                         <  2:1934990/100‑1 (MQ=255)
                                                        atcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt                                        >  1:1138804/1‑100 (MQ=255)
                                                        atcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt                                        >  1:295868/1‑100 (MQ=255)
                                                          caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                      <  1:14702/100‑1 (MQ=255)
                                                           aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                     >  2:81586/1‑100 (MQ=255)
                                                           aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                     >  1:2055085/1‑100 (MQ=255)
                                                            aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                    <  2:248973/100‑1 (MQ=255)
                                                                  ggCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgt                              >  1:652066/1‑100 (MQ=255)
                                                                   gCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACtgtg                             <  1:115186/100‑1 (MQ=255)
                                                                    cAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGa                            <  2:1469465/100‑1 (MQ=255)
                                                                      gTTCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                          >  1:431347/1‑100 (MQ=255)
                                                                       ttCCCAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                         <  1:295304/100‑1 (MQ=255)
                                                                       ttCCCAACAACAAAATAAAGATGGGAGGT‑AAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagga                        <  1:1243586/100‑1 (MQ=255)
                                                                         cccAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                       <  1:885582/100‑1 (MQ=255)
                                                                         cccAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                       >  1:2032117/1‑100 (MQ=255)
                                                                         cccAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggag                       <  1:724804/100‑1 (MQ=255)
                                                                          ccAACAACAAAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAaggagg                      >  2:1994196/1‑100 (MQ=255)
                                                                           caacaacaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                     >  1:933395/1‑100 (MQ=255)
                                                                           caacaacaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                     >  2:271608/1‑100 (MQ=255)
                                                                           caacaacaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGAAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                     >  1:1583193/1‑100 (MQ=255)
                                                                             acaacaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACg                   <  1:43354/100‑1 (MQ=255)
                                                                               aacaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGcc                 >  1:1278959/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                >  2:108112/1‑100 (MQ=255)
                                                                                acaaAATAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGAGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                >  2:1278975/1‑100 (MQ=255)
                                                                                  aaaaTAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATg              <  1:122895/100‑1 (MQ=255)
                                                                                      tAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc          >  1:2050541/1‑100 (MQ=255)
                                                                                      tAAAGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc          <  1:1388828/100‑1 (MQ=255)
                                                                                       aaaGATGGGAGGAAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTCTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt         <  2:1247177/100‑1 (MQ=255)
                                                                                         aGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt       <  2:1138802/100‑1 (MQ=255)
                                                                                         aGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt       <  2:295866/100‑1 (MQ=255)
                                                                                         aGATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt       >  1:562847/1‑100 (MQ=255)
                                                                                          gATGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt      <  1:2054579/100‑1 (MQ=255)
                                                                                           aTGGGAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTa     >  2:2163884/1‑100 (MQ=255)
                                                                                              ggAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATc  >  1:495932/1‑100 (MQ=255)
                                                                                              ggAGGTAAAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTCGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATc  >  1:986692/1‑100 (MQ=255)
                                                                                                    |                                                                                              
AGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGT‑AAA‑AATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATC  >  NC_000913/2563414‑2563606

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