Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,758,920 C→T 16.6% intergenic (+64/‑65) 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,9200CT16.6% 477.6 / 74.2 203intergenic (+64/‑65)alpA/yfjICP4‑57 prophage; DNA‑binding transcriptional activator/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (86/83);  new base T (21/13);  total (107/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 = 1.77e-02

Only 100 of 209 total aligned reads displayed.

GATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATT  >  NC_000913/2758821‑2759019
                                                                                                   |                                                                                                    
gatgaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGcc                                                                                                      <  1:615066/100‑1 (MQ=255)
  tgaGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCAt                                                                                                    <  1:4267053/100‑1 (MQ=255)
    aGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTa                                                                                                  >  2:1483631/1‑100 (MQ=255)
    aGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTa                                                                                                  <  2:2501990/100‑1 (MQ=255)
     gTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                 >  1:2808356/1‑100 (MQ=255)
     gTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTaa                                                                                                 <  1:1359824/100‑1 (MQ=255)
        gTTATTACAACGCTGACAACCTGCTTGAGGCCGTTACTTGAATTCAGACGCAGTTGCTGGCATCGGAGTGTAAACATCCTGCTATTTCCGACGTtagtag                                                                                              <  1:114654/100‑1 (MQ=255)
         ttATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGc                                                                                             >  2:3371801/1‑100 (MQ=255)
         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)
            tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑a                                                                                         >  2:560083/1‑100 (MQ=255)
               caACGCTGTAAACTTATTTGAGGGTGTTATATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGG‑aaat                                                                                      <  2:4194857/100‑3 (MQ=255)
               caACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGG‑aaat                                                                                      <  1:562539/100‑3 (MQ=255)
               caACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATa                                                                                      >  1:2147391/1‑100 (MQ=255)
                aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGG‑aaata                                                                                     <  2:633106/100‑4 (MQ=255)
                aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAg                                                                                     <  2:3002022/100‑1 (MQ=255)
                 aCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGc                                                                                    <  1:1944260/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)
                   gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTA‑TTCCGCTATTAATGGCGG‑aaatagcag                                                                                 <  2:4211635/100‑8 (MQ=37)
                   gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTA‑TTCCGCTATTAATGGCGG‑aaatagcag                                                                                 <  2:4367288/100‑8 (MQ=37)
                    cTGTAAACTTATTTTAAGGTGTTACATAAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGG‑aaatagcag                                                                                 <  1:3828044/100‑8 (MQ=37)
                     tGTAAACTTAGTTGAGTGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGT‑AATAGCCata                                                                                <  1:3571736/100‑1 (MQ=255)
                      gTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCatat                                                                               >  2:179936/1‑100 (MQ=255)
                       tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATAtt                                                                              <  2:863669/100‑1 (MQ=255)
                         aaCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTg                                                                            >  1:1903925/1‑100 (MQ=255)
                           cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCt                                                                          <  2:1441275/100‑1 (MQ=255)
                           cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCt                                                                          <  1:1034361/100‑1 (MQ=255)
                           cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCt                                                                          <  1:1435487/100‑1 (MQ=255)
                               tttGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAAc                                                                      <  1:2476026/100‑1 (MQ=255)
                                 tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACaa                                                                    >  2:2116350/1‑100 (MQ=255)
                                 tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACaa                                                                    >  1:88965/1‑100 (MQ=255)
                                   aGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAAtt                                                                  >  2:2646909/1‑100 (MQ=255)
                                    gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAAttt                                                                 >  2:2359566/1‑100 (MQ=255)
                                     ggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AGTAGCCATATTTGCTAAACAATTTc                                                                <  1:3763404/100‑1 (MQ=255)
                                      gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCt                                                               >  2:4306212/1‑100 (MQ=255)
                                        gtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCtta                                                             >  1:2296071/1‑100 (MQ=255)
                                            cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATCGCGG‑AATAGCCATACTTGCTAAACAATTTCttaa                                                            >  1:4090397/1‑97 (MQ=255)
                                            cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaattt                                                         >  2:2009236/1‑100 (MQ=255)
                                            cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaattt                                                         <  1:310859/100‑1 (MQ=255)
                                              taaaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                       >  2:2982407/3‑100 (MQ=255)
                                              tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                       >  2:933030/1‑100 (MQ=255)
                                              tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCttaatttaa                                                       >  1:4219872/1‑100 (MQ=255)
                                                aaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTtaatttaatt                                                     <  2:2689269/100‑1 (MQ=255)
                                                 aTTCATACTCAATTGCTGTCATCAGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATttaatta                                                    >  1:949356/1‑100 (MQ=255)
                                                    cATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                 >  2:2083176/1‑100 (MQ=255)
                                                    cATACTCAATTGCTGTCATAGGAGTGTAAACAGCCTGCTATTTCCGCCATTAATGGCGG‑AAAAGCCATATGTGCTAAACAACTTCTTAATTTAATTAAAt                                                 <  1:3195518/100‑1 (MQ=255)
                                                    cATACTCAATAGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATCGACATCTTTTCTAAACAATTTCttaatttaatcacat                                                 >  1:1291807/1‑95 (MQ=255)
                                                     ntACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATa                                                >  1:3717806/2‑100 (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:1394173/1‑100 (MQ=255)
                                                        cTCAATTGCTGCCATCGGAGTGTAACCATCCTTCTATTTCCGCCATTAATAGCGT‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                               <  1:4353180/98‑1 (MQ=255)
                                                         tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                            <  2:1173115/100‑1 (MQ=255)
                                                          cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAgg                                           >  1:1278290/1‑100 (MQ=255)
                                                           aaTTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCAGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGc                                          >  1:2364364/1‑100 (MQ=255)
                                                              tGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCAtt                                       >  2:3472258/1‑100 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  2:1647568/100‑1 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  1:1474301/100‑1 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  1:3179025/100‑1 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  2:3148672/100‑1 (MQ=255)
                                                                  gTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCAt                                   <  2:1544906/100‑1 (MQ=255)
                                                                     aTCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGtt                                >  1:1850136/1‑100 (MQ=255)
                                                                      tCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTc                               <  2:3636287/100‑1 (MQ=255)
                                                                       cGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                              <  2:198170/100‑1 (MQ=255)
                                                                        ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                             >  2:806157/1‑100 (MQ=255)
                                                                         gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                            >  2:3708832/1‑100 (MQ=255)
                                                                         gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                            <  1:1188812/100‑1 (MQ=255)
                                                                          aGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                           <  1:1009895/100‑1 (MQ=255)
                                                                           gtgtAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                          >  2:1116652/1‑100 (MQ=255)
                                                                           gtgtAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                          >  1:3310646/1‑100 (MQ=255)
                                                                            tgtAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgt                         >  2:1887117/1‑100 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAGTTTAAGTAAATAAGAGGCATTGCATGTTCAATGgtcg                       >  1:2970625/1‑100 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                       <  1:1785450/100‑1 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                       >  2:2821230/1‑100 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                       >  1:733212/1‑100 (MQ=255)
                                                                                    tCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCttt                 <  2:1687534/100‑1 (MQ=255)
                                                                                     ccTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                <  2:3755598/100‑1 (MQ=255)
                                                                                     ccTGCTATTTCCCCCATTAATAGCGG‑AATAGTCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                <  1:1102348/100‑1 (MQ=255)
                                                                                      atggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                >  2:657447/8‑100 (MQ=37)
                                                                                       tggctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTccc              >  2:3556755/7‑100 (MQ=39)
                                                                                         gctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg            >  2:1717225/5‑100 (MQ=255)
                                                                                         cTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg            <  1:2046822/100‑1 (MQ=255)
                                                                                          ctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg            >  2:3279779/4‑100 (MQ=38)
                                                                                          ctattCCGCTATTAATGGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt           >  1:2945033/4‑100 (MQ=255)
                                                                                            tttCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg         >  2:4036313/1‑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:707728/100‑1 (MQ=255)
                                                                                              tCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt       >  2:352774/1‑100 (MQ=255)
                                                                                                cGCCATTACTAGCGG‑AATAGCCATATTTGCTACACAACTTCTTAATTTAATTAAATAAGAGGCATTGCATGATCAATGGTCGTCCTTTCCCTGTAGATGc     >  1:187995/1‑100 (MQ=255)
                                                                                                cGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc     <  2:1470873/100‑1 (MQ=255)
                                                                                                cGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc     >  2:3919639/1‑100 (MQ=255)
                                                                                                cGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc     <  2:2247867/100‑1 (MQ=255)
                                                                                                   cATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCAtt  >  2:3829500/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
GATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGG‑AATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATT  >  NC_000913/2758821‑2759019

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