Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,758,927 A→G 16.6% intergenic (+71/‑58) alpA → / → yfjI CP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,758,9270AG16.6% 458.8 / 79.3 204intergenic (+71/‑58)alpA/yfjICP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (86/83);  new base G (21/13);  total (108/96)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

Only 100 of 206 total aligned reads displayed.

GTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAA  >  NC_000913/2758829‑2759026
                                                                                                  |                                                                                                    
 tcttCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGc                                                                                                    <  2:1499530/97‑1 (MQ=255)
  tATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCg                                                                                                   <  1:2214827/100‑1 (MQ=255)
    tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGa                                                                                                >  2:2436533/1‑100 (MQ=255)
    tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGa                                                                                                >  1:1644881/1‑100 (MQ=255)
        aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAg                                                                                            <  2:3002022/100‑1 (MQ=255)
         aCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTCCGCTATTAATGGCGGaaatagc                                                                                          <  1:2998279/100‑6 (MQ=38)
         aCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTCCGCTATTAATGGCGGaaatagc                                                                                          <  1:3393413/100‑6 (MQ=38)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaatagca                                                                                         <  1:3263390/100‑7 (MQ=39)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCa                                                                                         >  2:648873/1‑100 (MQ=255)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCa                                                                                         <  1:3780411/100‑1 (MQ=255)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTCCGCTATTAATGGCGGaaatagcag                                                                                        <  1:1861870/100‑8 (MQ=37)
            cTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAAAATCCTGCTATTCCGCTATTAATGGCGGaaatagcagg                                                                                       <  1:2352054/100‑9 (MQ=25)
             tGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCata                                                                                       <  1:2257083/100‑1 (MQ=255)
               tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATAtt                                                                                     <  2:863669/100‑1 (MQ=255)
                   cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCt                                                                                 <  2:1441275/100‑1 (MQ=255)
                   cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCt                                                                                 <  1:1435487/100‑1 (MQ=255)
                   cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCt                                                                                 >  2:581221/1‑100 (MQ=255)
                         tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACaa                                                                           >  2:4393818/1‑100 (MQ=255)
                         tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACaa                                                                           >  2:2116350/1‑100 (MQ=255)
                           aGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAtt                                                                         >  2:2646909/1‑100 (MQ=255)
                            gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAttt                                                                        <  2:4309291/100‑1 (MQ=255)
                            gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAttt                                                                        >  2:2359566/1‑100 (MQ=255)
                              gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCt                                                                      >  2:674344/1‑100 (MQ=255)
                              gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCt                                                                      <  2:2171486/100‑1 (MQ=255)
                                 ttACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaa                                                                   <  1:1133342/100‑1 (MQ=255)
                                    taTGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaattt                                                                <  1:926520/99‑1 (MQ=255)
                                    cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATCGCGGAATAGCCATACTTGCTAAACAATTTCttaa                                                                   >  1:4090397/1‑97 (MQ=255)
                                    cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaattt                                                                <  1:310859/100‑1 (MQ=255)
                                    cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaattt                                                                >  2:2009236/1‑100 (MQ=255)
                                      taaaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  2:2982407/3‑100 (MQ=255)
                                      tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  1:4219872/1‑100 (MQ=255)
                                      tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  2:3485424/1‑100 (MQ=255)
                                      tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  2:933030/1‑100 (MQ=255)
                                       gAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaat                                                             <  1:2132376/100‑1 (MQ=255)
                                        aaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTtaatttaatt                                                            <  2:2689269/100‑1 (MQ=255)
                                            cATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                        >  2:2083176/1‑100 (MQ=255)
                                            cATACTCAATTGCTGTCATAGGAGTGTAAACAGCCTGCTATTTCCGCCATTAATGGCGGAAAAGCCATATGTGCTAAACAACTTCTTAATTTAATTAAAt                                                        <  1:3195518/100‑1 (MQ=255)
                                             aTACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATa                                                       <  1:3749776/100‑1 (MQ=255)
                                              tACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATATCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                      >  1:3983549/1‑100 (MQ=255)
                                              tACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                      <  1:3323168/100‑1 (MQ=255)
                                                cTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAaga                                                    >  2:3684930/1‑100 (MQ=255)
                                                cTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAaga                                                    >  2:1394173/1‑100 (MQ=255)
                                                cTCAATTGCTGCCATCGGAGTGTAACCATCCTTCTATTTCCGCCATTAATAGCGTAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                      <  1:4353180/98‑1 (MQ=255)
                                                 tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                   <  2:3411143/100‑1 (MQ=255)
                                                 tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                   <  2:1173115/100‑1 (MQ=255)
                                                  cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAgg                                                  <  2:3549936/100‑1 (MQ=255)
                                                   aaTTGCTGTCATCGGAGTCTATACATCCTTCTATTTCCGCTATTAATAGCGAATTTGCCATATTTGCTAACCAATTTCTTAATTGAATTAAATAagag                                                   <  1:2371161/98‑1 (MQ=255)
                                                       gCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTg                                             >  2:1854175/1‑100 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  1:3179025/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  1:2041659/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  1:1474301/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  2:3148672/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  2:1647568/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  2:1693863/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGCAAACATCCTGGTATTTCCCTCATTAATAACTGAATAGACATATTTGCTAAACAATTTCTTAATTTAATTAAAGAAGAGGCATTGCa                                           <  1:1933297/100‑1 (MQ=255)
                                                          gTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCAt                                          <  2:1544906/100‑1 (MQ=255)
                                                            cATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATg                                         <  2:3559859/99‑1 (MQ=255)
                                                               cGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                                     >  1:1286521/1‑100 (MQ=255)
                                                               cGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                                     <  2:198170/100‑1 (MQ=255)
                                                                ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                    >  2:2619199/1‑100 (MQ=255)
                                                                ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                    >  2:175729/1‑100 (MQ=255)
                                                                ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                    >  2:806157/1‑100 (MQ=255)
                                                                 gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                   <  1:1188249/100‑1 (MQ=255)
                                                                  aGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                                  <  1:1009895/100‑1 (MQ=255)
                                                                   gtgtAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                                 >  1:3310646/1‑100 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAGTTTAAGTAAATAAGAGGCATTGCATGTTCAATGgtcg                              >  1:2970625/1‑100 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                              <  2:1683187/100‑1 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                              >  1:733212/1‑100 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                              >  2:2821230/1‑100 (MQ=255)
                                                                        aaCATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                                  >  1:1059621/1‑94 (MQ=255)
                                                                           aatatggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCt                          >  1:1234457/11‑100 (MQ=25)
                                                                             ccTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                       <  2:3755598/100‑1 (MQ=255)
                                                                              atggctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTcc                      >  1:1825170/8‑100 (MQ=38)
                                                                               tggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                       >  1:1888462/7‑99 (MQ=37)
                                                                                ggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTcc                      >  1:1782406/6‑99 (MQ=37)
                                                                                ggctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCt                    >  1:3579736/6‑100 (MQ=255)
                                                                                gCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCACt                    >  2:1492203/1‑100 (MQ=255)
                                                                                 gctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCt                    >  1:564263/5‑100 (MQ=38)
                                                                                 gctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg                   >  2:3307118/5‑100 (MQ=255)
                                                                                 gctattCCGCTATTAATGGCGGAATAGACATATTTGCTAAACCATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCACTTTTCGTCAtt                         >  1:3761867/5‑94 (MQ=25)
                                                                                  tATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt                  <  2:4170731/100‑1 (MQ=255)
                                                                                  ctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg                   >  2:3271332/4‑100 (MQ=38)
                                                                                  ctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt                  >  1:2945033/4‑100 (MQ=255)
                                                                                   tattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCATGTa                 >  1:2537615/3‑100 (MQ=255)
                                                                                   aTTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTa                 <  1:4357526/100‑1 (MQ=255)
                                                                                    tttCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  2:4036313/1‑100 (MQ=255)
                                                                                    attCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGTGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  2:4202593/2‑100 (MQ=255)
                                                                                    attCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  2:2901065/2‑100 (MQ=255)
                                                                                     ttCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  1:4242337/1‑100 (MQ=255)
                                                                                      tCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                <  1:1689008/98‑1 (MQ=255)
                                                                                      tCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt              >  2:1059206/1‑100 (MQ=255)
                                                                                      tCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt              <  2:707728/100‑1 (MQ=255)
                                                                                        cGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc            >  2:3919639/1‑100 (MQ=255)
                                                                                        cGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc            <  2:1470873/100‑1 (MQ=255)
                                                                                           cATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCAtt         >  2:3829500/1‑100 (MQ=255)
                                                                                           cATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCACTGTAGATGCAtt         >  2:430500/1‑100 (MQ=255)
                                                                                             ttAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATTTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTc       >  1:3922075/1‑100 (MQ=255)
                                                                                              tAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTcc      <  1:278473/100‑1 (MQ=255)
                                                                                              tAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTcc      <  1:1322688/100‑1 (MQ=255)
                                                                                                aTAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTa    >  1:536038/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
GTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAA  >  NC_000913/2758829‑2759026

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: