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.

GATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATT  >  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)
            tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGa                                                                                         >  2:2436533/1‑100 (MQ=255)
            tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGa                                                                                         >  1:1644881/1‑100 (MQ=255)
            tCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGa                                                                                         >  2:560083/1‑100 (MQ=255)
               caACGCTGTAAACTTATTTGAGGGTGTTATATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaat                                                                                      <  2:4194857/100‑3 (MQ=255)
               caACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaat                                                                                      <  1:562539/100‑3 (MQ=255)
               caACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATa                                                                                      >  1:2147391/1‑100 (MQ=255)
                aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaata                                                                                     <  2:633106/100‑4 (MQ=255)
                aaCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAg                                                                                     <  2:3002022/100‑1 (MQ=255)
                 aCGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGc                                                                                    <  1:1944260/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)
                   gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTCCGCTATTAATGGCGGaaatagcag                                                                                 <  2:4211635/100‑8 (MQ=37)
                   gCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTCCGCTATTAATGGCGGaaatagcag                                                                                 <  2:4367288/100‑8 (MQ=37)
                    cTGTAAACTTATTTTAAGGTGTTACATAAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCTATTAATGGCGGaaatagcag                                                                                 <  1:3828044/100‑8 (MQ=37)
                     tGTAAACTTAGTTGAGTGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGTAATAGCCata                                                                                <  1:3571736/100‑1 (MQ=255)
                      gTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCatat                                                                               >  2:179936/1‑100 (MQ=255)
                       tAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATAtt                                                                              <  2:863669/100‑1 (MQ=255)
                         aaCTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTg                                                                            >  1:1903925/1‑100 (MQ=255)
                           cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCt                                                                          <  2:1441275/100‑1 (MQ=255)
                           cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCt                                                                          <  1:1034361/100‑1 (MQ=255)
                           cTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCt                                                                          <  1:1435487/100‑1 (MQ=255)
                               tttGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAAc                                                                      <  1:2476026/100‑1 (MQ=255)
                                 tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACaa                                                                    >  2:2116350/1‑100 (MQ=255)
                                 tGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACaa                                                                    >  1:88965/1‑100 (MQ=255)
                                   aGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAtt                                                                  >  2:2646909/1‑100 (MQ=255)
                                    gggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAAttt                                                                 >  2:2359566/1‑100 (MQ=255)
                                     ggTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAGTAGCCATATTTGCTAAACAATTTc                                                                <  1:3763404/100‑1 (MQ=255)
                                      gtgtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCt                                                               >  2:4306212/1‑100 (MQ=255)
                                        gtTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCtta                                                             >  1:2296071/1‑100 (MQ=255)
                                            cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATCGCGGAATAGCCATACTTGCTAAACAATTTCttaa                                                            >  1:4090397/1‑97 (MQ=255)
                                            cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaattt                                                         >  2:2009236/1‑100 (MQ=255)
                                            cATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaattt                                                         <  1:310859/100‑1 (MQ=255)
                                              taaaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                       >  2:2982407/3‑100 (MQ=255)
                                              tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                       >  2:933030/1‑100 (MQ=255)
                                              tGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCttaatttaa                                                       >  1:4219872/1‑100 (MQ=255)
                                                aaTTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTtaatttaatt                                                     <  2:2689269/100‑1 (MQ=255)
                                                 aTTCATACTCAATTGCTGTCATCAGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATttaatta                                                    >  1:949356/1‑100 (MQ=255)
                                                    cATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAAt                                                 >  2:2083176/1‑100 (MQ=255)
                                                    cATACTCAATTGCTGTCATAGGAGTGTAAACAGCCTGCTATTTCCGCCATTAATGGCGGAAAAGCCATATGTGCTAAACAACTTCTTAATTTAATTAAAt                                                 <  1:3195518/100‑1 (MQ=255)
                                                    cATACTCAATAGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATCGACATCTTTTCTAAACAATTTCttaatttaatcacat                                                 >  1:1291807/1‑95 (MQ=255)
                                                     ntACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATa                                                >  1:3717806/2‑100 (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:1394173/1‑100 (MQ=255)
                                                        cTCAATTGCTGCCATCGGAGTGTAACCATCCTTCTATTTCCGCCATTAATAGCGTAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATaa                                               <  1:4353180/98‑1 (MQ=255)
                                                         tCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAagag                                            <  2:1173115/100‑1 (MQ=255)
                                                          cAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAgg                                           >  1:1278290/1‑100 (MQ=255)
                                                           aaTTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCAGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGc                                          >  1:2364364/1‑100 (MQ=255)
                                                              tGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCAtt                                       >  2:3472258/1‑100 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  2:1647568/100‑1 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  1:1474301/100‑1 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  1:3179025/100‑1 (MQ=255)
                                                                 tGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCa                                    <  2:3148672/100‑1 (MQ=255)
                                                                  gTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCAt                                   <  2:1544906/100‑1 (MQ=255)
                                                                     aTCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGtt                                >  1:1850136/1‑100 (MQ=255)
                                                                      tCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTc                               <  2:3636287/100‑1 (MQ=255)
                                                                       cGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCa                              <  2:198170/100‑1 (MQ=255)
                                                                        ggAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCaa                             >  2:806157/1‑100 (MQ=255)
                                                                         gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                            >  2:3708832/1‑100 (MQ=255)
                                                                         gAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAAt                            <  1:1188812/100‑1 (MQ=255)
                                                                          aGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATg                           <  1:1009895/100‑1 (MQ=255)
                                                                           gtgtAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                          >  2:1116652/1‑100 (MQ=255)
                                                                           gtgtAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATgg                          >  1:3310646/1‑100 (MQ=255)
                                                                            tgtAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgt                         >  2:1887117/1‑100 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAGTTTAAGTAAATAAGAGGCATTGCATGTTCAATGgtcg                       >  1:2970625/1‑100 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                       <  1:1785450/100‑1 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                       >  2:2821230/1‑100 (MQ=255)
                                                                              tAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGgtcg                       >  1:733212/1‑100 (MQ=255)
                                                                                    tCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCttt                 <  2:1687534/100‑1 (MQ=255)
                                                                                     ccTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                <  2:3755598/100‑1 (MQ=255)
                                                                                     ccTGCTATTTCCCCCATTAATAGCGGAATAGTCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                <  1:1102348/100‑1 (MQ=255)
                                                                                      atggctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTc                >  2:657447/8‑100 (MQ=37)
                                                                                       tggctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTccc              >  2:3556755/7‑100 (MQ=39)
                                                                                         gctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg            >  2:1717225/5‑100 (MQ=255)
                                                                                         cTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg            <  1:2046822/100‑1 (MQ=255)
                                                                                          ctattCCGCTATTAATGGCGGAAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTg            >  2:3279779/4‑100 (MQ=38)
                                                                                          ctattCCGCTATTAATGGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGt           >  1:2945033/4‑100 (MQ=255)
                                                                                            tttCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAg         >  2:4036313/1‑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:707728/100‑1 (MQ=255)
                                                                                              tCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGAt       >  2:352774/1‑100 (MQ=255)
                                                                                                cGCCATTACTAGCGGAATAGCCATATTTGCTACACAACTTCTTAATTTAATTAAATAAGAGGCATTGCATGATCAATGGTCGTCCTTTCCCTGTAGATGc     >  1:187995/1‑100 (MQ=255)
                                                                                                cGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc     <  2:1470873/100‑1 (MQ=255)
                                                                                                cGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc     >  2:3919639/1‑100 (MQ=255)
                                                                                                cGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGc     <  2:2247867/100‑1 (MQ=255)
                                                                                                   cATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCAtt  >  2:3829500/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
GATGAGTGGTTATCACAACGCTGTAAACTTATTTGAGGGTGTTACATGAATTCATACTCAATTGCTGTCATCGGAGTGTAAACATCCTGCTATTTCCGCCATTAATAGCGGAATAGCCATATTTGCTAAACAATTTCTTAATTTAATTAAATAAGAGGCATTGCATGTTCAATGGTCGTCCTTTCCCTGTAGATGCATT  >  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: