Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,969 C→A 7.3% M5I (ATG→ATT)  yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,9690CA 7.3% 1081.1 / 9.2 347M5I (ATG→ATT) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (183/126);  new base A (3/34);  total (187/160)
Fisher's exact test for biased strand distribution p-value = 1.54e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 349 total aligned reads displayed.

CCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGA  >  NC_000913/281820‑282118
                                                                                                                                                     |                                                                                                                                                     
  ttCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCAt                                                                                                                                                     >  1:5684084/1‑150 (MQ=255)
    ccTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATaa                                                                                                                                                   <  1:386550/150‑1 (MQ=255)
     cTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAg                                                                                                                                                  <  2:5098218/150‑1 (MQ=255)
     cTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAg                                                                                                                                                  <  2:733707/150‑1 (MQ=255)
     cTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAg                                                                                                                                                  >  1:3133517/1‑150 (MQ=255)
      tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGc                                                                                                                                                 >  1:2778515/1‑150 (MQ=255)
      tGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGc                                                                                                                                                 >  2:3657350/1‑150 (MQ=255)
                 cTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATaaa                                                                                                                                      >  1:841764/1‑150 (MQ=255)
                    gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATCAACgg                                                                                                                                   >  1:4782478/1‑150 (MQ=255)
                    gAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACgg                                                                                                                                   >  2:6403855/1‑150 (MQ=255)
                       ccACTTGTTGCGGGTGGCAGGTGAAATGCCGAAGCGAGGGCTGTGGGAAGGGTTGTTCGCCCCGTCCTGCTTGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGtt                                                                                                                                <  2:5275829/150‑1 (MQ=255)
                       ccACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGtt                                                                                                                                >  1:4477060/1‑150 (MQ=255)
                            tGTAGCGGGTGGCTGTTGAAATGCCTAAGCTAGGGCATAGGGAACGGATGTTCACCCGGTCCTGCTTGGAGAACAAAACAAACTCAGTACGTAATGACTTGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTAt                                                                                                                           <  2:5952046/150‑1 (MQ=255)
                            tGTAGCCGGTGGCTGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCACCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTAt                                                                                                                           <  2:1079934/150‑1 (MQ=255)
                              ttgcgggTGGCTGGTGATTGGCCGACGCGTGGGCTGTGGTAAGGGTTGTTCGCCCGGTCGTGGGTGGCGAACAAAACAACCTCAGTCGGTATTGACTGGGTATCTCTCGCCGCCCGGGGCATAAGCGACTCCATAAACGGGTTCTTATGc                                                                                                                         <  2:3001811/145‑1 (MQ=255)
                                 ccGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCtt                                                                                                                      >  2:2380194/1‑150 (MQ=255)
                                  cGGTGGCAGGTGAAATGCCGAAGCGAGGGCAGAGGGAACGTATGTTCGCCCGGTCCTGCGAGGCGAACAAAAAAAACTCATTACTTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTa                                                                                                                     <  2:2082169/150‑1 (MQ=255)
                                  cGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTa                                                                                                                     >  2:4500281/1‑150 (MQ=255)
                                   tgTGGCAGGTGAAATGCGGAAGCGACGGCTGTGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAg                                                                                                                    <  2:352391/149‑1 (MQ=255)
                                   ggTGGCTGGTTAACGGCCGAAGCGACGGCAGTGGGAAGGGATGTTCGCCCCGTCCTGCTAGGCTAAAAACACAAACCCCGTACGTACTGACCTGGTATCCCTCGAATCCCAGGTCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAg                                                                                                                    <  2:5702033/150‑1 (MQ=255)
                                    gTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGt                                                                                                                   >  2:1933128/1‑150 (MQ=255)
                                    gTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCGTAGt                                                                                                                   >  1:4907614/1‑150 (MQ=255)
                                       gCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGATCTTATGCCTTAGTTGt                                                                                                                >  1:4188562/1‑150 (MQ=255)
                                         aGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAACACAAACTCAGTACGTAATGACATGGTCTCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCGTATGCCGTAGGTGTaa                                                                                                              >  1:6462547/1‑150 (MQ=255)
                                           gTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgt                                                                                                            >  1:5103997/1‑150 (MQ=255)
                                               aaTGCCGAAGCCACGGCAGAGGCAACGGAAGTACGCCCCGTCCTCCGAGGCGAACAAAACAAACTCAGTACGTAAAGACAACGTATCTCTCGCATCCCAGGGCATAAGCCACACCAAAAACGCGTTCTTATGCCTTAGTTGTAACTGACt                                                                                                        >  1:4762953/1‑150 (MQ=255)
                                                aTGCCGACGCGACGGCAGAGGCAACGGCTGTGCGCCCGGCCCTGCGAGGCTAACAAAAAAAACTCCGTACGTAATGACATGGTATCCCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTa                                                                                                       <  1:4296977/150‑1 (MQ=255)
                                                aTGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTa                                                                                                       >  1:4029489/1‑150 (MQ=255)
                                                  gCCGACGCGACGGCTGTGGGAAGGGTTGTTCGCCCCGTCCTGCGTGGCTAAAAAAACAACCCCAGTACTTAATGACATGGTACCTCTCGCACCCCCGGGAAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAcc                                                                                                     <  1:4117339/150‑1 (MQ=255)
                                                   ccGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCa                                                                                                    >  2:4290334/1‑150 (MQ=255)
                                                    cGACGCGCCGGCTGTGGGACCGGTTGTTCGCCCCGTCCTGCGTGGCCAACAAAAAAAACTCAGTACGTAATGACTTGGTATCTCTCGCATCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAt                                                                                                   <  1:4731280/150‑1 (MQ=255)
                                                    cGAAGCGACGGCAGAGGTAACGGATGTTCGCCCCGTCCTGCGTGGCGAACAAAACAAACTCAGTACGTATTGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAt                                                                                                   <  2:3436860/150‑1 (MQ=255)
                                                          gccGGCAGTGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAAAAAAACAAACCCAGTACTTAATGCCAGGGTATCTCGCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTcccc                                                                                             <  1:2045581/148‑1 (MQ=255)
                                                            cGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAAAAAAACAAACTCAGTACGTAATGACTTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGa                                                                                           <  2:4197402/150‑1 (MQ=255)
                                                              gCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTAAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTGTTATGCCTTAGTTGTAAGTGTCTACCATGTCACCGAAc                                                                                         >  1:1694992/1‑150 (MQ=255)
                                                               cAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACa                                                                                        >  2:1478969/1‑150 (MQ=255)
                                                                  agGGAACGTATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACAGGGTATCTCCCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgt                                                                                     <  2:1204684/150‑1 (MQ=255)
                                                                   ggtaaCGTATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGTTATCTCTCGCATCCCCGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgtg                                                                                    <  2:2299085/147‑1 (MQ=255)
                                                                     gAACGGAGGTGCGCCCCGTCCTGCGAGGCGAACAAAACAACCTCAGTACGTAATGACATGGTATCTCTCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGtt                                                                                  <  2:5762044/150‑1 (MQ=255)
                                                                      aaCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTGCTTCTGCCTGAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTc                                                                                 >  1:3651267/1‑150 (MQ=255)
                                                                      aaCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGGTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTc                                                                                 >  2:5446273/1‑150 (MQ=255)
                                                                         ggATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACt                                                                              >  2:2611282/1‑150 (MQ=255)
                                                                          taTGTTCGCCCCGTCCTGCGTGGCTAACAAAAAAAACACAGTACTTAATGACATGTTATCTCTCGAATCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTa                                                                             <  2:1334345/149‑1 (MQ=255)
                                                                           aTGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATACACGGGTTCTTAGGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTAt                                                                            >  1:2143873/1‑150 (MQ=255)
                                                                              ttCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTc                                                                         >  2:6286638/1‑150 (MQ=255)
                                                                               tCGCCCCTTCCTGCGAGGAGAAAAAAAAAAACTCAGTACTTAATGACATGGTATCTCTCGAATCCCCGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTcc                                                                        <  2:5212746/150‑1 (MQ=255)
                                                                                cGCCCCGTCCTGCGAGGCTAACAAAAAAAACTCAGTACTTAATGACTGGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTccc                                                                       <  1:5397856/150‑1 (MQ=255)
                                                                                 gCCCCGTCCTGCGTGGCTAACAAAACAAACCCAGTACGTAATGACTTGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTcccc                                                                      <  1:2652669/150‑1 (MQ=255)
                                                                                  ccccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCg                                                                     >  1:6291897/1‑150 (MQ=255)
                                                                                  ccccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCg                                                                     >  2:1171425/1‑150 (MQ=255)
                                                                                  ccccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCGTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCg                                                                     >  1:3278425/1‑150 (MQ=255)
                                                                                   cccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCgg                                                                    >  1:1013834/1‑150 (MQ=255)
                                                                                   cccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCgg                                                                    >  1:5256399/1‑150 (MQ=255)
                                                                                    ccGTCCTGCGCGGCGAACAAAAAAAACAAAGTACTTAATGACATGGTATCCCTCGCACCCCTGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGa                                                                   <  1:2309746/150‑1 (MQ=255)
                                                                                    ccGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCGTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGa                                                                   >  2:929373/1‑150 (MQ=255)
                                                                                        ccTGCGCGGCGAACAAAACAAACTCTGTCCGTACGGACCTGGTATCTCTCGCATCCCCGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGt                                                               <  2:3599227/150‑1 (MQ=255)
                                                                                          tGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTac                                                             >  1:3841696/1‑150 (MQ=255)
                                                                                            cGAGGCGAACAAACCAAACCCAGTACGTAATGACATGGTATCTCTCGCAGCCCAGGGCATACGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACAC                                                           <  1:742979‑M1/150‑1 (MQ=255)
                                                                                               ggCGAAAAAAACAAACTCCGTACGTAATGACTTGTTATCTCTCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCC                                                        <  1:4811966‑M1/150‑1 (MQ=255)
                                                                                                gCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCA                                                       >  2:58649‑M1/1‑150 (MQ=255)
                                                                                                  gAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAA                                                     >  1:1491979‑M1/1‑150 (MQ=255)
                                                                                                   aacaaaacaaACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAG                                                    <  1:5599142‑M1/150‑1 (MQ=255)
                                                                                                     caaaacaacCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATACACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGG                                                  <  2:4395623‑M1/150‑1 (MQ=255)
                                                                                                     caaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCAGCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGG                                                  <  2:3312205‑M1/150‑1 (MQ=255)
                                                                                                     aaaaacaaacTCCGTACTTAATGACATGGTATCTCTCGCAGCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGG                                                  <  1:2420383‑M1/149‑1 (MQ=255)
                                                                                                      aaaacaaaCTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGG                                                 >  1:1029511‑M1/1‑150 (MQ=255)
                                                                                                       aaacaaccTCAGTACGTAAGGACTTGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGga                                                <  2:4138729/150‑1 (MQ=255)
                                                                                                        aacaaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGgag                                               >  2:2745708/1‑150 (MQ=255)
                                                                                                          caaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGagag                                             <  2:4588111/150‑1 (MQ=255)
                                                                                                           aaacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAgg                                            >  1:1559724/1‑150 (MQ=255)
                                                                                                            aacTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAggg                                           >  1:662500/1‑150 (MQ=255)
                                                                                                            aacTCACTACTAAATGACATGGTATATCTCGAAACCCAGGCCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAggg                                           <  2:3211652/150‑1 (MQ=255)
                                                                                                                 aGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTGTGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTg                                      >  2:2321480/1‑150 (MQ=255)
                                                                                                                    aCGTAATGACTTGGTATCTCTCGCAGCCCCGGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCAc                                   <  2:820388/150‑1 (MQ=255)
                                                                                                                     cGTAATGACAGGGTATCTCTCGCATCCCCGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCAcc                                  <  1:3437469/150‑1 (MQ=255)
                                                                                                                      gTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCg                                 >  2:4753576/1‑150 (MQ=255)
                                                                                                                        aaTGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGGTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAg                               >  2:373857/1‑150 (MQ=255)
                                                                                                                          tGACATGGTATCTCTCGCATCCCTGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGcc                             <  2:3263690/150‑1 (MQ=255)
                                                                                                                             aaTGGTTTCTTTCGAATCGCTGGGAATAAGCTTCTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGTGACTGCACCGAGCGATc                          <  2:6625388/149‑1 (MQ=255)
                                                                                                                                  tATCTCTCGCAACCCCGGGCAAAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTc                     <  1:1275876/150‑1 (MQ=255)
                                                                                                                                   aTCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTcc                    <  2:2136838/150‑1 (MQ=255)
                                                                                                                                     ctctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTcccc                  >  2:3073258/1‑150 (MQ=255)
                                                                                                                                     ctctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTcccc                  <  1:3531940/150‑1 (MQ=255)
                                                                                                                                      tctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccccc                 >  1:4561541/1‑150 (MQ=255)
                                                                                                                                      tctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccccc                 >  2:401767/1‑150 (MQ=255)
                                                                                                                                      tctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTccccc                 >  1:5031287/1‑150 (MQ=255)
                                                                                                                                       ctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTAcc                                                                                                     <  2:482418/65‑1 (MQ=255)
                                                                                                                                       ctcGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCt                >  2:2867957/1‑150 (MQ=255)
                                                                                                                                          gCATCCCCGGGCATAAGCGACTCCATAACCGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGc             >  2:3965336‑M2/1‑149 (MQ=255)
                                                                                                                                            aacccagGGAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccc           >  2:1694992‑M2/3‑147 (MQ=255)
                                                                                                                                                cAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGccccttt       <  1:2365695‑M2/150‑8 (MQ=255)
                                                                                                                                                 aGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttg      >  1:5995942‑M2/1‑142 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATACACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  2:435769‑M2/1‑141 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  1:1761267‑M2/1‑141 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  2:1012320‑M2/1‑141 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  1:6286638‑M2/1‑141 (MQ=255)
                                                                                                                                                  gggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  2:4875986‑M2/1‑141 (MQ=255)
                                                                                                                                                  gggAATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttgg     >  1:3979767‑M2/1‑141 (MQ=255)
                                                                                                                                                   ggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggg    <  1:901994‑M2/150‑11 (MQ=255)
                                                                                                                                                   ggCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGcccctttggg    >  1:2461154‑M2/1‑140 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACCCCAAAGGGGAGAGGGGACTGCACCGAGCCATCTTTTCCCCCTCGCCCCTTTGGGGA  >  NC_000913/281820‑282118

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