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.

GTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAA  >  NC_000913/2758829‑2759026
                                                                                                  |                                                                                                    
 tcttCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGc                                                                                                    <  2:1499530/97‑1 (MQ=255)
  tATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCg                                                                                                   <  1:2214827/100‑1 (MQ=255)
    tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑a                                                                                                >  2:2436533/1‑100 (MQ=255)
    tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑a                                                                                                >  1:1644881/1‑100 (MQ=255)
        aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAg                                                                                            <  2:3002022/100‑1 (MQ=255)
         aCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTA‑TTCCGCTATTAATGGCGG‑aaatagc                                                                                          <  1:2998279/100‑6 (MQ=38)
         aCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTA‑TTCCGCTATTAATGGCGG‑aaatagc                                                                                          <  1:3393413/100‑6 (MQ=38)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGG‑aaatagca                                                                                         <  1:3263390/100‑7 (MQ=39)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCa                                                                                         >  2:648873/1‑100 (MQ=255)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCa                                                                                         <  1:3780411/100‑1 (MQ=255)
           gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTA‑TTCCGCTATTAATGGCGG‑aaatagcag                                                                                        <  1:1861870/100‑8 (MQ=37)
            cTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAAAATCCTGCTA‑TTCCGCTATTAATGGCGG‑aaatagcagg                                                                                       <  1:2352054/100‑9 (MQ=25)
             tGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCata                                                                                       <  1:2257083/100‑1 (MQ=255)
               tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATAtt                                                                                     <  2:863669/100‑1 (MQ=255)
                   cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCt                                                                                 <  2:1441275/100‑1 (MQ=255)
                   cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCt                                                                                 <  1:1435487/100‑1 (MQ=255)
                   cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCt                                                                                 >  2:581221/1‑100 (MQ=255)
                         tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACaa                                                                           >  2:4393818/1‑100 (MQ=255)
                         tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACaa                                                                           >  2:2116350/1‑100 (MQ=255)
                           aGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAAtt                                                                         >  2:2646909/1‑100 (MQ=255)
                            gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAAttt                                                                        <  2:4309291/100‑1 (MQ=255)
                            gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAAttt                                                                        >  2:2359566/1‑100 (MQ=255)
                              gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCt                                                                      >  2:674344/1‑100 (MQ=255)
                              gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCt                                                                      <  2:2171486/100‑1 (MQ=255)
                                 ttACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaa                                                                   <  1:1133342/100‑1 (MQ=255)
                                    taTGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaattt                                                                <  1:926520/99‑1 (MQ=255)
                                    cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATCGCGG‑AATAGCCATACTTGCTAAACAATTTCttaa                                                                   >  1:4090397/1‑97 (MQ=255)
                                    cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaattt                                                                <  1:310859/100‑1 (MQ=255)
                                    cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaattt                                                                >  2:2009236/1‑100 (MQ=255)
                                      taaaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  2:2982407/3‑100 (MQ=255)
                                      tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  1:4219872/1‑100 (MQ=255)
                                      tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  2:3485424/1‑100 (MQ=255)
                                      tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                              >  2:933030/1‑100 (MQ=255)
                                       gAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaat                                                             <  1:2132376/100‑1 (MQ=255)
                                        aaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTtaatttaatt                                                            <  2:2689269/100‑1 (MQ=255)
                                            cATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                        >  2:2083176/1‑100 (MQ=255)
                                            cATACTCAATTGCTGTCATAGGAGTGTAAACAGCCTGCTATTTCCGCCATTAATGGCGG‑AAAAGCCATATGTGCTAAACAACTTCTTAATTTAATTAAAt                                                        <  1:3195518/100‑1 (MQ=255)
                                             aTACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATa                                                       <  1:3749776/100‑1 (MQ=255)
                                              tACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATATCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                      >  1:3983549/1‑100 (MQ=255)
                                              tACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                      <  1:3323168/100‑1 (MQ=255)
                                                cTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAaga                                                    >  2:3684930/1‑100 (MQ=255)
                                                cTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAaga                                                    >  2:1394173/1‑100 (MQ=255)
                                                cTCAATTGCTGCCATCGGAGTGTAACCATCCTTCTATTTCCGCCATTAATAGCGT‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                                      <  1:4353180/98‑1 (MQ=255)
                                                 tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                   <  2:3411143/100‑1 (MQ=255)
                                                 tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                                   <  2:1173115/100‑1 (MQ=255)
                                                  cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAgg                                                  <  2:3549936/100‑1 (MQ=255)
                                                   aaTTGCTGTCATCGGAGTCTATACATCCTTCTATTTCCGCTATTAATAGCGA‑ATTTGCCATATTTGCTAACCAATTTCTTAATTGAATTAAATAagag                                                   <  1:2371161/98‑1 (MQ=255)
                                                       gCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTg                                             >  2:1854175/1‑100 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  1:3179025/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  1:2041659/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  1:1474301/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  2:3148672/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  2:1647568/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                           <  2:1693863/100‑1 (MQ=255)
                                                         tGTCATCGGAGTGCAAACATCCTGGTATTTCCCTCATTAATAACTG‑AATAGACATATTTGCTAAACAATTTCTTAATTTAATTAAAGAAGAGGCATTGCa                                           <  1:1933297/100‑1 (MQ=255)
                                                          gTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCAt                                          <  2:1544906/100‑1 (MQ=255)
                                                            cATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATg                                         <  2:3559859/99‑1 (MQ=255)
                                                               cGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                                     >  1:1286521/1‑100 (MQ=255)
                                                               cGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                                     <  2:198170/100‑1 (MQ=255)
                                                                ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                    >  2:2619199/1‑100 (MQ=255)
                                                                ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                    >  2:175729/1‑100 (MQ=255)
                                                                ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                                    >  2:806157/1‑100 (MQ=255)
                                                                 gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                                   <  1:1188249/100‑1 (MQ=255)
                                                                  aGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                                  <  1:1009895/100‑1 (MQ=255)
                                                                   gtgtAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                                 >  1:3310646/1‑100 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAGTTTAAGTAAATAAGAGGCATTGCATGTTCAATGgtcg                              >  1:2970625/1‑100 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                              <  2:1683187/100‑1 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                              >  1:733212/1‑100 (MQ=255)
                                                                      tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                              >  2:2821230/1‑100 (MQ=255)
                                                                        aaCATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                                  >  1:1059621/1‑94 (MQ=255)
                                                                           aatatggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCt                          >  1:1234457/11‑100 (MQ=25)
                                                                             ccTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                       <  2:3755598/100‑1 (MQ=255)
                                                                              atggctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTcc                      >  1:1825170/8‑100 (MQ=38)
                                                                               tggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                       >  1:1888462/7‑99 (MQ=37)
                                                                                ggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTcc                      >  1:1782406/6‑99 (MQ=37)
                                                                                ggctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCt                    >  1:3579736/6‑100 (MQ=255)
                                                                                gCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCACt                    >  2:1492203/1‑100 (MQ=255)
                                                                                 gctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCt                    >  1:564263/5‑100 (MQ=38)
                                                                                 gctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg                   >  2:3307118/5‑100 (MQ=255)
                                                                                 gctattCCGCTATTAATGGCGG‑AATAGACATATTTGCTAAACCATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCACTTTTCGTCAtt                         >  1:3761867/5‑94 (MQ=25)
                                                                                  tATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt                  <  2:4170731/100‑1 (MQ=255)
                                                                                  ctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg                   >  2:3271332/4‑100 (MQ=38)
                                                                                  ctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt                  >  1:2945033/4‑100 (MQ=255)
                                                                                   tattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCATGTa                 >  1:2537615/3‑100 (MQ=255)
                                                                                   aTTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTa                 <  1:4357526/100‑1 (MQ=255)
                                                                                    tttCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  2:4036313/1‑100 (MQ=255)
                                                                                    attCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGTGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  2:4202593/2‑100 (MQ=255)
                                                                                    attCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  2:2901065/2‑100 (MQ=255)
                                                                                     ttCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                >  1:4242337/1‑100 (MQ=255)
                                                                                      tCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg                <  1:1689008/98‑1 (MQ=255)
                                                                                      tCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt              >  2:1059206/1‑100 (MQ=255)
                                                                                      tCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt              <  2:707728/100‑1 (MQ=255)
                                                                                        cGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc            >  2:3919639/1‑100 (MQ=255)
                                                                                        cGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc            <  2:1470873/100‑1 (MQ=255)
                                                                                           cATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCAtt         >  2:3829500/1‑100 (MQ=255)
                                                                                           cATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCACTGTAGATGCAtt         >  2:430500/1‑100 (MQ=255)
                                                                                             ttAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATTTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTc       >  1:3922075/1‑100 (MQ=255)
                                                                                              tAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTcc      <  1:278473/100‑1 (MQ=255)
                                                                                              tAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTcc      <  1:1322688/100‑1 (MQ=255)
                                                                                                aTAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTa    >  1:536038/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
GTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATTTCCTAAA  >  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: ‑