Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 581,394 T→C 6.6% pseudogene (41/309 nt) aaaD → pseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913581,3940TC 6.6% 160.9 / 3.5 62pseudogene (41/309 nt)aaaDpseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (33/23);  new base C (2/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

GCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAA  >  NC_000913/581247‑581538
                                                                                                                                                   |                                                                                                                                                
gcgTTTTCCGGGGCTGCAAAGCCGACATGTTGATTGCCGGAAACGGCATATGATCAGAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGCGCAGTCTGGTTAACAGGCTGCGGCCtttt                                                                                                                                                 <  2:576286/149‑1 (MQ=14)
  gTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTCATACCGGGGTTGCTGAGTGAATCTATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGt                                                                                                                                               <  2:649931/149‑1 (MQ=255)
        cGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTCTGTCCgcgc                                                                                                                                         >  1:864204/1‑149 (MQ=255)
        cGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTCTGTCCgcgc                                                                                                                                         >  1:393701/1‑149 (MQ=255)
         ggAACTGGAAAACGGCCATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCAGCGCTGGACGAACTGATACCGGGGCTGCTGAGTGGATATGTCGAACAGGCAGGTGTAAAGGCTGCGGCCTTGCGTCCGAGcc                                                                                                                                        >  2:752958/1‑149 (MQ=14)
           aaCTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCGGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAAGATATCGAACAGTCATGTTAACAGGCTGCGGCATTTCGTCCGCCCCgg                                                                                                                                      >  2:505532/1‑149 (MQ=18)
            aCTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCggg                                                                                                                                     <  1:316120/149‑1 (MQ=255)
                 aaaaCCGACATGTTGATTTCCTGAAACCGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCg                                                                                                                                <  1:809656/149‑1 (MQ=255)
                   aatcGACATGTTGACTGCCTGAAACGGGATATCAGCAAAGCCATGAACAAAGCAGCCTCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCt                                                                                                                              <  2:856751/146‑1 (MQ=21)
                    aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTc                                                                                                                             >  2:482873/1‑149 (MQ=255)
                    aCCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTc                                                                                                                             >  2:651320/1‑149 (MQ=255)
                     ccGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATAACGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCa                                                                                                                            <  1:503060/149‑1 (MQ=255)
                       gACATGTTGATTTCCTGAAGCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACt                                                                                                                          <  1:788241/149‑1 (MQ=255)
                       gACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACt                                                                                                                          <  1:257205/149‑1 (MQ=255)
                            gTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCCCTGTTCa                                                                                                                     >  1:413826/1‑149 (MQ=255)
                             ttGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCGTTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAg                                                                                                                    >  2:342622/1‑149 (MQ=255)
                              tGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAgg                                                                                                                   <  1:774130/149‑1 (MQ=255)
                                    ccTGAAACGGGATAGCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAg                                                                                                             <  1:15869/149‑1 (MQ=255)
                                      tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGcc                                                                                                           >  1:837625/1‑149 (MQ=255)
                                      tGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGcc                                                                                                           >  1:101222/1‑149 (MQ=255)
                                       gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTCTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCca                                                                                                          <  2:864204/149‑1 (MQ=255)
                                       gAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTCTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCca                                                                                                          <  2:393701/149‑1 (MQ=255)
                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCCCTGTTAAGGCCGGAGCcac                                                                                                         >  2:907817/1‑149 (MQ=255)
                                        aaaCGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCcac                                                                                                         >  2:98585/1‑149 (MQ=255)
                                         aaCGGGATATCATCAAAGCCAGGAACAAAGCAGCCGCGCCGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCcaca                                                                                                        <  1:253390/149‑1 (MQ=255)
                                              gATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATCCACTGCTCCCGGGGTTGCTGAGTGAATATCTCGAACAGTCAGGCCAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACCGTTCAGGCCGGAGCCACAGAccg                                                                                                   >  2:831063/1‑149 (MQ=25)
                                                  tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGt                                                                                               >  1:712478/1‑149 (MQ=255)
                                                  tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGt                                                                                               >  1:219258/1‑149 (MQ=255)
                                                  tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGt                                                                                               <  2:546267/149‑1 (MQ=255)
                                                  tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGAAAGGATACGCTAACAGTTCAAGCAGGAGCCAAAGAtggccgg                                                                                               >  1:191646/1‑148 (MQ=21)
                                                  tcatcaAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCAGGGGTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGt                                                                                               >  1:712440/1‑149 (MQ=255)
                                                       aaaGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAGCATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTCCGCGCCCTGTTCAGTCCCGAGCCACAGACCCCCGTTGAAt                                                                                          >  1:840772/1‑149 (MQ=37)
                                                       aaaGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAAt                                                                                          <  1:59728/149‑1 (MQ=255)
                                                                    aaGCACCCTCGCTGGATGTACTGCTACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGCCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTaa                                                                             <  1:616041/149‑1 (MQ=39)
                                                                       cagcCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGACCAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTa                                                                          <  2:357354/149‑1 (MQ=255)
                                                                       cagcCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAAttc                                                                          >  2:385411/1‑148 (MQ=255)
                                                                                       aCTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGaa                                                          >  1:896686/1‑149 (MQ=255)
                                                                                         tGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCCGGTTAACAGGCTGCGGCATTTTGTACGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATc                                                        >  1:624506/1‑149 (MQ=255)
                                                                                         tGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATc                                                        >  1:257683/1‑149 (MQ=255)
                                                                                          gacaCCGGGGTTGCTGAGTGAATATCTCGAACCGTCAGGTTAACAGGCTGCGGCATTTTGGCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATcc                                                       <  2:689299/146‑1 (MQ=255)
                                                                                           aTACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCg                                                      >  2:815730/1‑149 (MQ=255)
                                                                                             aCCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCa                                                    >  1:406796/1‑149 (MQ=255)
                                                                                               cGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATa                                                  <  2:219258/149‑1 (MQ=255)
                                                                                                ggggTTGCTGAGTGAATATATAGAACAGACAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATAc                                                 >  2:377600/1‑149 (MQ=255)
                                                                                                 gggTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCCCCCGAAAGAATCCGCCTAcc                                                >  2:745938/1‑149 (MQ=255)
                                                                                                       cTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAg                                          <  2:286385/149‑1 (MQ=255)
                                                                                                        tGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAgg                                         >  1:362830/1‑149 (MQ=255)
                                                                                                        tGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTACACTCCCTGTACAGGCCGGAGCCACAGACCGCCGCAGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAgc                                         >  2:148323/1‑148 (MQ=255)
                                                                                                         gggTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGCCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAggg                                        <  2:89672/147‑1 (MQ=255)
                                                                                                            tGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGgcgc                                     >  1:586360/1‑149 (MQ=255)
                                                                                                            tGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGccgc                                     >  1:743/1‑149 (MQ=255)
                                                                                                             gAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCt                                    >  1:564096/1‑149 (MQ=255)
                                                                                                                  tatCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGaa                               <  2:697998/149‑1 (MQ=255)
                                                                                                                  tatCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGaa                               <  1:793167/149‑1 (MQ=255)
                                                                                                                  tatCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTAACGGTTCAGGCCGGAGCCACAGACCGCCGTTCAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGCAGCTGGGaa                               >  2:378083/1‑149 (MQ=255)
                                                                                                                   atCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGCAAGAATCCGCATACCAGGAGGGGCGCTGGGaaa                              >  2:169702/1‑149 (MQ=255)
                                                                                                                      gAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTGTCTCCCGAAAGAATCCGCATACCAAGAAGGGCGCTGGGAAacgc                           >  2:75466/1‑147 (MQ=255)
                                                                                                                      gAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAacac                           >  2:792205/1‑149 (MQ=255)
                                                                                                                         cAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGc                        <  1:227366/149‑1 (MQ=255)
                                                                                                                              aGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCttt                   <  1:98585/149‑1 (MQ=255)
                                                                                                                               ggTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTc                  >  2:258134/1‑149 (MQ=255)
                                                                                                                                      aGGCTGCGGCATGGTGTCTGCGGCGCGCCTCGCTCACTGTTCAGGCCGGAGCCACCGACCGCCGTTGAATGGGCGGATGCTAATAACTACCCACCGAAAGAAGCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGc           <  2:29912/149‑1 (MQ=255)
                                                                                                                                               cATTTTGTCCGCGCCGGGCTTCGCAACCTCGTCATGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCCCCCGAAAGAATCCGCATACCAGGTAGGGCGCTGGGAAACAATGCCCTTTCAGCGGGCCATCTTGaa  <  2:191646/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAAGCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATCGAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAA  >  NC_000913/581247‑581538

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