Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,892,514 T→C intergenic (+58/+62) b1810 → / ← b1811 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,892,5140TC97.4% 463.6 / ‑0.8 152intergenic (+58/+62)b1810/b1811hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base C (80/68);  total (84/68)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 158 total aligned reads displayed.

TATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCA  >  NC_000913/1892424‑1892604
                                                                                          |                                                                                          
tatatGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGccc                                                                                          >  1:1744875/1‑93 (MQ=255)
 atatGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCt                                                                                         >  2:2282080/1‑93 (MQ=255)
   atGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCttt                                                                                       >  2:4191789/1‑93 (MQ=255)
   atGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCttt                                                                                       <  1:2898123/93‑1 (MQ=255)
   atGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCttt                                                                                       <  2:4663850/93‑1 (MQ=255)
    tGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCtttt                                                                                      >  1:2052584/1‑93 (MQ=255)
      cAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtc                                                                                    >  2:2649374/1‑93 (MQ=255)
       aGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtct                                                                                   >  1:1694158/1‑93 (MQ=255)
       aGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtct                                                                                   >  2:2491008/1‑93 (MQ=255)
       aGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtct                                                                                   >  1:2255434/1‑93 (MQ=255)
       aGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtct                                                                                   >  2:161305/1‑93 (MQ=255)
        ggAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtctc                                                                                  <  2:4524236/93‑1 (MQ=255)
        ggAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtctc                                                                                  >  2:308426/1‑93 (MQ=255)
        ggAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTtctc                                                                                  >  2:4094799/1‑93 (MQ=255)
         gaagaaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTca                                                                                 >  1:4005904/1‑93 (MQ=255)
         gaagaaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTca                                                                                 >  1:3038230/1‑93 (MQ=255)
         gaagaaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTca                                                                                 <  1:4090303/93‑1 (MQ=255)
            gaaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACaa                                                                              >  1:231873/1‑93 (MQ=255)
            gaaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACaa                                                                              <  1:1078094/93‑1 (MQ=255)
                tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAACAAGCGCCCTTTTTCTCACAATCaa                                                                          >  1:2950681/1‑93 (MQ=255)
                tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCaa                                                                          >  2:2388915/1‑93 (MQ=255)
                tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCaa                                                                          >  2:1503901/1‑93 (MQ=255)
                 aaaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCaaa                                                                         <  1:421030/93‑1 (MQ=255)
                 aaaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCaaa                                                                         <  1:3667855/93‑1 (MQ=255)
                  aaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAAt                                                                        <  1:800701/93‑1 (MQ=255)
                  aaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAAt                                                                        <  1:2998561/93‑1 (MQ=255)
                   aaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATg                                                                       >  2:4696005/1‑93 (MQ=255)
                   aaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATg                                                                       >  1:3241329/1‑93 (MQ=255)
                   aaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATg                                                                       <  2:2217782/93‑1 (MQ=255)
                    aGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGa                                                                      <  2:2955672/93‑1 (MQ=255)
                    aGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGa                                                                      <  2:3759669/93‑1 (MQ=255)
                     gCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGat                                                                     >  1:4140667/1‑93 (MQ=255)
                          ggTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCg                                                                <  2:3052840/93‑1 (MQ=255)
                          ggTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCg                                                                <  2:3427628/93‑1 (MQ=255)
                            tCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGc                                                              <  1:2736884/93‑1 (MQ=255)
                              tGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCtt                                                            >  2:22986/1‑93 (MQ=255)
                              tGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCtt                                                            >  2:3837653/1‑93 (MQ=255)
                              tGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCtt                                                            >  1:2639222/1‑93 (MQ=255)
                              tGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCtt                                                            >  1:2685345/1‑93 (MQ=255)
                                 gTTAAATTTATATCAGCTTAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCta                                                         <  1:4222765/93‑1 (MQ=255)
                                   tAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCtata                                                       <  1:131742/93‑1 (MQ=255)
                                        ttATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGaa                                                  >  1:4353461/1‑93 (MQ=255)
                                          atatCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAAt                                                >  2:3353841/1‑93 (MQ=255)
                                          atatCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAAt                                                <  1:4687223/93‑1 (MQ=255)
                                            atCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCg                                              >  2:1228539/1‑93 (MQ=255)
                                             tCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGa                                             <  2:3653803/93‑1 (MQ=255)
                                             tCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGa                                             <  1:2178245/93‑1 (MQ=255)
                                              cAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGAt                                            <  2:518749/93‑1 (MQ=255)
                                               aGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATa                                           >  1:384341/1‑93 (MQ=255)
                                               aGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATa                                           <  1:3469234/93‑1 (MQ=255)
                                                gCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACATTCAAATGATCATCGGCTTCTTTAACGAAATCGATAc                                          >  2:3829713/1‑93 (MQ=255)
                                                gCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATAc                                          <  2:2210439/93‑1 (MQ=255)
                                                gCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATAc                                          <  1:715725/93‑1 (MQ=255)
                                                   tAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGa                                       >  1:558831/1‑93 (MQ=255)
                                                    aaaTGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGAt                                      <  1:3493395/93‑1 (MQ=255)
                                                     aaTGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGAtt                                     >  2:4588564/1‑93 (MQ=255)
                                                     aaTGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGAtt                                     >  2:2797408/1‑93 (MQ=255)
                                                      aTGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTa                                    <  2:266056/93‑1 (MQ=255)
                                                       tGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAg                                   >  2:3738835/1‑93 (MQ=255)
                                                       tGGGTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAg                                   <  2:3309997/93‑1 (MQ=255)
                                                        gggTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgt                                  <  1:2526660/93‑1 (MQ=255)
                                                          gTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgtgt                                <  1:642179/93‑1 (MQ=255)
                                                          gTCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgtgt                                >  1:2035832/1‑93 (MQ=255)
                                                           tCAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgtgtg                               <  2:1402924/93‑1 (MQ=255)
                                                            cAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgtgtgt                              <  1:2616447/93‑1 (MQ=255)
                                                            cAGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgtgtgt                              <  1:4414642/93‑1 (MQ=255)
                                                             aGGACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGtt                             >  1:2072543/1‑93 (MQ=255)
                                                              ggACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTa                            >  1:2220685/1‑93 (MQ=255)
                                                               gACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTAc                           <  2:777749/93‑1 (MQ=255)
                                                               gACGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTAc                           >  1:3486967/1‑93 (MQ=255)
                                                                aCGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACt                          <  1:755446/93‑1 (MQ=255)
                                                                 ctcTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTc                         <  2:675756/91‑1 (MQ=255)
                                                                 cGCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTc                         <  1:4497886/93‑1 (MQ=255)
                                                                  gCTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCg                        <  2:3251390/93‑1 (MQ=255)
                                                                   cTTTTAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGt                       >  2:1314368/1‑93 (MQ=255)
                                                                    ttttAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTc                      >  1:2614744/1‑93 (MQ=255)
                                                                       tAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTc                   <  2:1894390/93‑1 (MQ=255)
                                                                       tAATCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTc                   <  1:4279148/93‑1 (MQ=255)
                                                                        aaTCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTct                  >  2:2545630/1‑92 (MQ=255)
                                                                        aaTCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCa                  <  1:1987406/93‑1 (MQ=255)
                                                                        aaTCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCa                  <  1:3833595/93‑1 (MQ=255)
                                                                        aaTCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCa                  >  2:27240/1‑93 (MQ=255)
                                                                          tCACATATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATc                <  1:3424392/93‑1 (MQ=255)
                                                                           cacaTATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATcc               >  1:2392032/1‑93 (MQ=255)
                                                                            acaTATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCt              >  2:4429934/1‑93 (MQ=255)
                                                                             caTATAAAAAGCGTACTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACTAAATCGATACGGATTAGAGTGTTACTCATCTTCATCCta             >  2:2088930/1‑92 (MQ=255)
                                                                             caTATAAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTc             <  1:3671846/93‑1 (MQ=255)
                                                                               tataAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaa           >  1:4281605/1‑93 (MQ=255)
                                                                               tataAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaa           >  2:1813278/1‑93 (MQ=255)
                                                                               tataAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaa           <  1:289606/93‑1 (MQ=255)
                                                                               tataAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaa           <  1:289605/93‑1 (MQ=255)
                                                                                 taAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaaaa         <  1:440635/93‑1 (MQ=255)
                                                                                 taAAAAGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaaaa         >  2:449764/1‑93 (MQ=255)
                                                                                  aaaaaGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTTATCCTCAAAAc        >  1:609838/1‑93 (MQ=255)
                                                                                  aaaaaGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAAc        >  2:4157708/1‑93 (MQ=255)
                                                                                   aaaaGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTCCTCGTCTTCATCCTCAAAAcg       <  1:2454769/93‑1 (MQ=255)
                                                                                    aaaGCGCCCTTGTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGCTTAGAGTGTTACTCGTCTTCATCCTCAAAActc      >  2:4429383/1‑91 (MQ=255)
                                                                                     aaGCGCCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAAcgcg     >  1:2388073/1‑93 (MQ=255)
                                                                                       gcgcCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGcc   >  1:3625527/1‑93 (MQ=255)
                                                                                       gcgcCCTTTTTCTCACAATCAAATGATCATCGGCTGCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATTCTCAAAACGCGcc   >  1:840013/1‑93 (MQ=255)
                                                                                          |                                                                                          
TATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCA  >  NC_000913/1892424‑1892604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: