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.A98.5% 472.5 / ‑0.3 136intergenic (+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 (72/62);  total (73/63)
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.86e-01

Only 100 of 155 total aligned reads displayed.

GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563415‑2563607
                                                                                                   |                                                                                              
ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                                                               >  2:75273/1‑100 (MQ=255)
ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTa                                                                                               >  2:2230696/1‑100 (MQ=255)
 gTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaa                                                                                              <  1:2283937/100‑1 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                             >  1:1047730/1‑100 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                             <  1:4137821/100‑1 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                             >  2:1749340/1‑100 (MQ=255)
  tGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                                                             >  2:242454/1‑100 (MQ=255)
   gAATTTTGTGATGCAGATCGCTTTTTTTTTTTTTTTTTTTTTGTTTAATTCAATCAATGAAAGGCAGTTTTCTACTACAAAATAAAGATTGGTGGTaaa                                                                                             >  1:3804262/1‑99 (MQ=255)
    aaTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaa                                                                                           <  2:3591487/100‑1 (MQ=255)
    aaTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaa                                                                                           <  2:534993/100‑1 (MQ=255)
     aTTTTGTGATTCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTACCTAACAACAAAATAAAGATGGGAGGTaaaaaa                                                                                          <  1:2807767/100‑2 (MQ=255)
     aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaa                                                                                          >  1:3464897/1‑99 (MQ=255)
         tgtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTaaaaaatgaa                                                                                      <  2:3893151/100‑6 (MQ=255)
           tgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAAc                                                                                     <  2:1123558/100‑1 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                    >  2:621974/1‑100 (MQ=255)
            gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACg                                                                                    >  2:3415714/1‑100 (MQ=255)
               gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTg                                                                                 <  2:1349089/100‑1 (MQ=255)
                 aGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTggg                                                                               >  1:2145678/1‑100 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTggg                                                                               >  2:2586/1‑99 (MQ=255)
                  gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGa                                                                              >  2:3560533/1‑100 (MQ=255)
                   aTCGCTTTTTCCCTTCAGGCATTTTTTTTTTTTTTTTTCAACGAAAAGCAGTTCCCAACAATAAATTAATGTTGGGAGGTAAAAAATGAAACGTTGTGa                                                                              >  1:1500407/1‑99 (MQ=255)
                   aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGaa                                                                             >  1:2225324/1‑100 (MQ=255)
                    tCGCTTTTTCCCTGCATGCAGTTTTGTTTAATTCCATCAACGAAATGCAGTTTCCAACAACATACTAAAGATGGGAGGTAAAAAATGAAACGGTgg                                                                                <  1:2522183/96‑1 (MQ=255)
                        tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGc                                                                        <  2:84834/100‑1 (MQ=255)
                         ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCt                                                                       <  2:1352101/100‑1 (MQ=255)
                          tttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCtt                                                                      <  1:3567230/100‑1 (MQ=255)
                           ttCCCTTCATGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCttt                                                                     <  2:1943540/100‑1 (MQ=255)
                            tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTa                                                                    >  2:433415/1‑100 (MQ=255)
                             cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTAc                                                                   <  2:632537/100‑1 (MQ=255)
                              ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACt                                                                  <  2:4036974/100‑1 (MQ=255)
                              ccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACt                                                                  <  2:2830536/100‑1 (MQ=255)
                               cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTg                                                                 >  2:1729080/1‑100 (MQ=255)
                               cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTg                                                                 >  2:3563179/1‑100 (MQ=255)
                                ttcctgCATTTTTATTTAATTCACTCAACGAAAGGCAGCTTCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  1:3824813/95‑1 (MQ=255)
                                ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                                                                <  2:3491071/100‑1 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  1:1567676/1‑100 (MQ=255)
                                  cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              <  2:2203686/100‑1 (MQ=255)
                                   aGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                                                              >  1:1557147/1‑99 (MQ=255)
                                    ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCGAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            >  1:2623295/1‑100 (MQ=255)
                                    ggCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            >  1:1345056/1‑100 (MQ=255)
                                      cATTTTTGTTTAATTCAATCAACTAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGAt                                                            <  1:3742911/98‑1 (MQ=255)
                                       aTTTTTGTGTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                                                         <  1:2300063/100‑1 (MQ=255)
                                           ttGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                                                       >  1:3942445/1‑98 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   >  2:512079/1‑100 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGtt                                                   <  2:1316148/100‑1 (MQ=255)
                                             gTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                                                    >  2:2511528/1‑99 (MQ=255)
                                               ttAATTCAATCAACCAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAACCGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                                                    <  1:3709083/97‑1 (MQ=255)
                                                tAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTgg                                                 >  2:2870254/1‑99 (MQ=255)
                                                 aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               >  2:2527008/1‑100 (MQ=255)
                                                 aaTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt                                               >  2:2835006/1‑100 (MQ=255)
                                                  aTTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt                                              >  2:1985362/1‑100 (MQ=255)
                                                    tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  1:2512396/100‑1 (MQ=255)
                                                    tcaatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGc                                            <  1:802188/100‑1 (MQ=255)
                                                     caatcaaCGAAGGGCAGTTCCCAACAACAAGATAAATAAGGGAGGTAAAAAATGAAACGGTAGGAAGTGGCTTTACTGGTGATGGTATGTTGGTTTGCTg                                          >  1:1758891/1‑100 (MQ=255)
                                                     caatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt                                           <  2:2263034/100‑1 (MQ=255)
                                                      aatcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg                                          <  1:1767682/100‑1 (MQ=255)
                                                        tcaaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa                                        <  2:271779/100‑1 (MQ=255)
                                                         caaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAt                                       >  2:3876306/1‑100 (MQ=255)
                                                          aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                      >  2:2756775/1‑100 (MQ=255)
                                                          aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                      <  1:68616/100‑1 (MQ=255)
                                                          aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                      >  1:2811072/1‑100 (MQ=255)
                                                          aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                      >  1:4071777/1‑100 (MQ=255)
                                                          aaCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTAtt                                      <  2:695494/100‑1 (MQ=255)
                                                           aCGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     >  1:1132838/1‑100 (MQ=255)
                                                           aCGAAAGGCAGTTCCCAACAACAAAATAAAGAGGGGAGGTAAAAAATGAAACGGTGGGAAGTGTCTTTAGTGGTTATGGTAATGTTGGTTTGCTGTATTg                                     >  1:1848694/1‑100 (MQ=255)
                                                             gAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTg                                     <  2:64574/98‑1 (MQ=255)
                                                              aaaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  >  1:797308/1‑100 (MQ=255)
                                                              aaaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAAc                                  >  1:1004413/1‑100 (MQ=255)
                                                               aaGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACt                                 <  2:3024603/100‑1 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGaa                            >  1:3127074/1‑100 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGaa                            >  2:2088637/1‑100 (MQ=255)
                                                                    aGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGAGGGAAGTGGCTTTACTGGTGATGATAATGTTGGTTTGCTGTATTGAACTGTGaa                            >  1:3019925/1‑100 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           <  2:4034007/100‑1 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           >  2:1501893/1‑100 (MQ=255)
                                                                     gTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAag                           >  1:208673/1‑100 (MQ=255)
                                                                      ttCCCAACAACAAAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAagg                          <  1:1176599/100‑1 (MQ=255)
                                                                          caacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                      <  2:3691594/100‑1 (MQ=255)
                                                                          caacaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAggagga                      >  2:1856289/1‑100 (MQ=255)
                                                                            acaacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACg                    >  2:1245471/1‑100 (MQ=255)
                                                                             caacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  2:1571741/1‑100 (MQ=255)
                                                                             caacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGc                   >  2:1208262/1‑100 (MQ=255)
                                                                              aacaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGcc                  <  2:2993141/100‑1 (MQ=255)
                                                                               acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 <  2:2991916/100‑1 (MQ=255)
                                                                               acaaAATAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCa                 >  2:3802935/1‑100 (MQ=255)
                                                                                  aaaTAAAGATGTGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGa              >  2:3398081/1‑100 (MQ=255)
                                                                                     tAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  2:146395/100‑1 (MQ=255)
                                                                                     tAAAGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAAc           <  2:1773553/100‑1 (MQ=255)
                                                                                      aaaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt          <  1:433413/100‑1 (MQ=255)
                                                                                      aaaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt          >  2:2410558/1‑100 (MQ=255)
                                                                                      aaaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACt          <  2:3270471/100‑1 (MQ=255)
                                                                                       aaGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtt         >  2:1319291/1‑100 (MQ=255)
                                                                                        aGATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACttt        <  1:2088633/100‑1 (MQ=255)
                                                                                         gATGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACtttt       >  1:4250830/1‑100 (MQ=255)
                                                                                           tGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  2:1574545/100‑1 (MQ=255)
                                                                                           tGGGAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTaa     <  2:731888/100‑1 (MQ=255)
                                                                                            gggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    <  2:3003395/100‑1 (MQ=255)
                                                                                            gggAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAAt    <  2:3245348/100‑1 (MQ=255)
                                                                                              gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  >  2:378352/1‑100 (MQ=255)
                                                                                              gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  <  1:474694/100‑1 (MQ=255)
                                                                                              gAGGTAAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCg  >  2:2634724/1‑100 (MQ=255)
                                                                                                   |                                                                                              
GGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATTGAACTGTGAAGGAGGACGCCATGAAACTTTTAATCG  >  NC_000913/2563415‑2563607

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: