Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 272,712 (G)5→6 intergenic (+523/‑135) insN → / → ykfC pseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/CP4‑6 prophage; conserved protein;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,7071.G98.2% 188.3 / ‑2.8 54intergenic (+518/‑140)insN/ykfCpseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/CP4‑6 prophage; conserved protein;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (24/29);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

AAACACCGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAAGACCTTTG  >  NC_000913/272561‑272849
                                                                                                                                                   |                                                                                                                                              
aaaCACCGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAgggg                                                                                                                                            >  2:10737/1‑151 (MQ=255)
 aaCACCGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggg                                                                                                                                           <  1:69236/151‑1 (MQ=255)
   cacCGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggggc                                                                                                                                         >  2:44091/1‑149 (MQ=255)
     ccGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggggcat                                                                                                                                       <  1:196254/151‑5 (MQ=255)
      cGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAggggggcatg                                                                                                                                      <  2:3029/151‑6 (MQ=255)
       ggggCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATgg                                                                                                                                      >  1:415675/1‑151 (MQ=255)
          gcagcaGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTc                                                                                                                                   >  2:525527/1‑151 (MQ=255)
           cagcagCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTcc                                                                                                                                  <  2:649898/151‑1 (MQ=255)
               agcGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATg                                                                                                                              <  2:781799/151‑1 (MQ=255)
                      aTCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACgcgc                                                                                                                        <  1:10716/150‑1 (MQ=255)
                      aTCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCgct                                                                                                                       <  2:472184/151‑1 (MQ=255)
                                        aCTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCgt                                                                                                     >  2:249051/1‑151 (MQ=255)
                                         cTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCgt                                                                                                     >  2:62996/1‑150 (MQ=255)
                                              gCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCgtagta                                                                                                 <  2:258976/149‑1 (MQ=255)
                                               ccAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTc                                                                                              <  1:169198/151‑1 (MQ=255)
                                                   cTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCcgcg                                                                                          >  2:120202/1‑151 (MQ=255)
                                                       cGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATgg                                                                                      >  2:518604/1‑151 (MQ=255)
                                                        ggCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATggg                                                                                     <  1:479000/151‑1 (MQ=255)
                                                         gcaagcaaGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATgggg                                                                                    >  1:237718/1‑151 (MQ=255)
                                                           aagcaagGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGa                                                                                   >  2:233955/1‑150 (MQ=255)
                                                             gcaagGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCt                                                                                                          <  1:280498/125‑1 (MQ=255)
                                                             gcaagGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCt                                                                                                          >  2:280523/1‑125 (MQ=255)
                                                                agGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGccc                                                                             <  1:44069/151‑1 (MQ=255)
                                                                   aTAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCAtt                                                                          <  1:677625/151‑1 (MQ=255)
                                                                    tAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCAtt                                                                          >  2:690235/1‑150 (MQ=255)
                                                                     aaaGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTa                                                                         <  1:710468/150‑1 (MQ=255)
                                                                            gCGACGGGCAGACTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGa                                                                 <  1:203417/151‑1 (MQ=255)
                                                                              gACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGa                                                                 >  2:91355/1‑149 (MQ=255)
                                                                               aCGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAggg                                                             <  2:428878/151‑1 (MQ=255)
                                                                                cGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAggg                                                             <  2:388364/151‑1 (MQ=255)
                                                                                 gggCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAggg                                                             >  1:686066/1‑150 (MQ=255)
                                                                                 gggCAGCCACCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGa                                                            <  2:311290/151‑1 (MQ=255)
                                                                                   gCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAg                                                               <  1:362374/146‑1 (MQ=255)
                                                                                   gCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAg                                                               >  2:362397/1‑146 (MQ=255)
                                                                                      gCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAgg                                                              <  2:85016/144‑1 (MQ=255)
                                                                                      gCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAgg                                                              >  1:84993/1‑144 (MQ=255)
                                                                                             cAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAAt                                                <  2:581668/151‑1 (MQ=255)
                                                                                             cAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAAt                                                <  1:766150/151‑1 (MQ=255)
                                                                                                   gCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGtt                                           <  1:545829/150‑1 (MQ=255)
                                                                                                    ccTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTg                                         <  2:406143/151‑1 (MQ=255)
                                                                                                         gTCCCGGCAGGTAAACCGGGGAAGGACAGCGACGGATGTTAAGGGGGGCATGTCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATCGCGAAGTGAAGCAGTTTGAATGTGTTTTCGACg                                    <  1:639023/151‑1 (MQ=255)
                                                                                                            caggcaggTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGt                                   <  1:800832/149‑1 (MQ=255)
                                                                                                            caggcaggTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGa                                 <  2:788746/151‑1 (MQ=255)
                                                                                                             aggcaggTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGa                                 >  2:614546/1‑150 (MQ=255)
                                                                                                                  ggTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAAc                           >  2:43520/1‑151 (MQ=255)
                                                                                                                  ggTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAAc                           >  2:180657/1‑151 (MQ=255)
                                                                                                                     aaaCCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTg                         >  1:343478/1‑150 (MQ=255)
                                                                                                                      aaCCGGGGAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGAc                       <  2:601705/151‑1 (MQ=255)
                                                                                                                          ggggAAGGTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTa                    >  2:347606/1‑150 (MQ=255)
                                                                                                                                ggTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGt             <  1:521842/151‑1 (MQ=255)
                                                                                                                                 gTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTg            <  1:270148/151‑1 (MQ=255)
                                                                                                                                 gTCAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTg            <  2:106862/151‑1 (MQ=255)
                                                                                                                                   cAGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGaa          <  1:271080/151‑1 (MQ=255)
                                                                                                                                    aGCGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGa           >  2:748731/1‑149 (MQ=255)
                                                                                                                                    aGCGACGGATGTTAAGGGGGGCATGGCGCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGaa          <  1:784644/150‑1 (MQ=255)
                                                                                                                                      cGACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAAGAc       >  1:763809/1‑151 (MQ=255)
                                                                                                                                       gACGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAAGAcc      >  1:369996/1‑151 (MQ=255)
                                                                                                                                         cGGATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAAGACCtt    >  1:363011/1‑151 (MQ=255)
                                                                                                                                           gATGTTAAGGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAAGACCTTTg  >  1:733092/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AAACACCGGGGCAGCAGCGTCTATCACGCTTGCGTTGCTGACTTCTGCCAACTTGCGGCAAGCAAGGATAAAGAGTGCGACGGGCAGCCTCCTCAGTATGCCTGAGTCCAGGCAGGTAAACCGGGGAAGGTCAGCGACGGATGTTAAGGGGGCATGGCTCCGATGACGCGCTGGCTGGCGGAGCTTCCGTAGTAGTCCGCGATGGGGAAAGCCCATTACATGGCGAAGGGAAGCAGTTTGAATGTGTTTGCGACGTGAATTAACTGACCTAACGAGGTGAAGACCTTTG  >  NC_000913/272561‑272849

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