Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,992 A→G intergenic (‑145/+13) ylaB ← / ← ylaC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,9920AG99.0% 635.1 / ‑2.9 195intergenic (‑145/+13)ylaB/ylaChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (95/98);  total (95/100)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 214 total aligned reads displayed.

TGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/477902‑478081
                                                                                          |                                                                                         
tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGtgtg                                                                                         >  1:1349131/1‑93 (MQ=255)
tggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGtgtg                                                                                         >  1:372103/1‑93 (MQ=255)
 ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGgtgtgt                                                                                        <  1:87120/93‑1 (MQ=255)
 ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGgtgtgt                                                                                        >  2:1214269/1‑93 (MQ=255)
 ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGgtgtgt                                                                                        >  1:1811143/1‑93 (MQ=255)
   agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTaa                                                                                      <  1:1881594/93‑1 (MQ=255)
      tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAAccc                                                                                   >  1:345784/1‑93 (MQ=255)
       aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCt                                                                                  <  2:600675/93‑1 (MQ=255)
       aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCt                                                                                  >  1:1832651/1‑93 (MQ=255)
         ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTaa                                                                                >  1:3072736/1‑93 (MQ=255)
          gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGGTTTACGCCTGGTGTGTAACCCTAAc                                                                               >  2:878968/1‑93 (MQ=255)
            ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACtt                                                                             >  1:3677679/1‑93 (MQ=255)
            ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACtt                                                                             >  1:3363904/1‑93 (MQ=255)
            ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACtt                                                                             <  1:3032761/93‑1 (MQ=255)
             gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTGACCCTAACTTa                                                                            >  2:3201859/1‑93 (MQ=255)
               agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAg                                                                          >  1:2932870/1‑93 (MQ=255)
               agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAg                                                                          >  2:2780776/1‑93 (MQ=255)
                gagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGa                                                                         <  1:566466/93‑1 (MQ=255)
                  gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCATAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATg                                                                       >  2:2209123/1‑93 (MQ=255)
                    tgctATTTCATGAAGTGATGGGTGGACACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     <  1:692480/90‑1 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     >  2:3466772/1‑93 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     >  2:4173698/1‑93 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     <  1:1659756/93‑1 (MQ=255)
                     ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTg                                                                    >  2:1946035/1‑93 (MQ=255)
                      gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGa                                                                   >  2:2577827/1‑93 (MQ=255)
                       tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAt                                                                  >  1:616567/1‑93 (MQ=255)
                       tATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAt                                                                  >  2:4306438/1‑93 (MQ=255)
                        aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAtt                                                                 >  2:826705/1‑93 (MQ=255)
                         tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTc                                                                >  2:3354024/1‑93 (MQ=255)
                         tttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTc                                                                <  1:2725960/93‑1 (MQ=255)
                           tCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCgg                                                              <  1:2482721/93‑1 (MQ=255)
                             aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcg                                                            <  1:2397761/93‑1 (MQ=255)
                              tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                           <  2:3943349/93‑1 (MQ=255)
                                aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcga                                                         <  1:1180764/93‑1 (MQ=255)
                                 aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGcgcgcg                                                        >  2:4064447/1‑91 (MQ=255)
                                 aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGcgcgcg                                                        >  1:2508344/1‑91 (MQ=255)
                                 aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                        <  2:2447644/93‑1 (MQ=255)
                                  gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                       <  2:2312698/93‑1 (MQ=255)
                                  gTGATGGGTGGAAACGATCATTGTCATGTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                       <  1:3662185/93‑1 (MQ=255)
                                    gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcgg                                                     >  2:3842296/1‑92 (MQ=255)
                                    gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                     <  2:2544374/93‑1 (MQ=255)
                                     aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                    >  1:3436965/1‑93 (MQ=255)
                                     aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                    >  1:2577488/1‑93 (MQ=255)
                                     aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                    <  1:4232975/93‑1 (MQ=255)
                                     aTGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                    <  2:4531890/93‑1 (MQ=255)
                                      tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   <  1:1873948/93‑1 (MQ=255)
                                      tGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   >  1:3829809/1‑93 (MQ=255)
                                        ggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                 <  1:3880578/93‑1 (MQ=255)
                                          tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaa                                               >  2:1989654/1‑93 (MQ=255)
                                           ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                              <  1:3846903/93‑1 (MQ=255)
                                             aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                            <  1:1062484/93‑1 (MQ=255)
                                               aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                          >  1:1910478/1‑93 (MQ=255)
                                                cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                         <  2:3037635/93‑1 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        <  1:1390507/93‑1 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        <  2:1444379/93‑1 (MQ=255)
                                                  aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                       >  2:2093299/1‑93 (MQ=255)
                                                  aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                       >  2:4436513/1‑93 (MQ=255)
                                                  aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                       >  2:1316928/1‑93 (MQ=255)
                                                   tCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGt                                      >  1:2696244/1‑93 (MQ=255)
                                                    cATGGTCATTTGACCAGGCGCAGAGCTTTACGCTTGGTGTGTACCCTTACCTTAAGATGTTGATTGGCGCGAGCGAGGGTAAAAATATCGTaa                                    <  2:3818219/93‑1 (MQ=255)
                                                     aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                    <  1:767776/93‑1 (MQ=255)
                                                        gTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                 <  2:3776872/93‑1 (MQ=255)
                                                         tCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                                <  1:945224/93‑1 (MQ=255)
                                                          cATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                               <  2:1224546/93‑1 (MQ=255)
                                                          cATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                               <  2:2048539/93‑1 (MQ=255)
                                                          cATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                               >  2:1826556/1‑93 (MQ=255)
                                                           aTTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                              <  2:3464514/93‑1 (MQ=255)
                                                            tttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             >  1:4577402/1‑93 (MQ=255)
                                                            tttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             >  1:4118354/1‑93 (MQ=255)
                                                             ttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                            <  1:1972334/93‑1 (MQ=255)
                                                             ttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                            <  1:494157/93‑1 (MQ=255)
                                                             ttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                            <  2:2670479/93‑1 (MQ=255)
                                                             ttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                            <  2:4278567/93‑1 (MQ=255)
                                                               gACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                          >  2:3519795/1‑93 (MQ=255)
                                                                  cAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       >  2:166210/1‑93 (MQ=255)
                                                                    ggCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTACAAATATCGTAAAAAGACAGTTCACCt                     >  1:98382/1‑93 (MQ=255)
                                                                     gcgcAGAACTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                    >  1:462037/1‑93 (MQ=255)
                                                                      cgcAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                   <  1:1461080/93‑1 (MQ=255)
                                                                       gcAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                  >  2:62365/1‑93 (MQ=255)
                                                                       gcAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                  >  1:4415854/1‑93 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 <  2:2976420/93‑1 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 <  2:3766719/93‑1 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 >  1:148616/1‑93 (MQ=255)
                                                                         agagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                >  2:1896417/1‑93 (MQ=255)
                                                                         agagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                >  2:4242986/1‑93 (MQ=255)
                                                                          gagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg               >  2:1956410/1‑93 (MQ=255)
                                                                           agCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa              <  2:3294731/93‑1 (MQ=255)
                                                                            gCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGAc             >  1:1423596/1‑93 (MQ=255)
                                                                              tttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa           >  1:4048013/1‑93 (MQ=255)
                                                                              tttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTACACCTTTACGGACGa           >  2:3998387/1‑93 (MQ=255)
                                                                                tACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc         <  1:695667/93‑1 (MQ=255)
                                                                                tACGCCGGGGGTGTAACCCTAACTTAAGATGTTGAGTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc         <  1:155311/93‑1 (MQ=255)
                                                                                 tcGCCTAGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa        <  2:1482043/92‑1 (MQ=255)
                                                                                   gCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtt      <  2:775731/93‑1 (MQ=255)
                                                                                    ccTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttt     <  2:4742438/93‑1 (MQ=255)
                                                                                     cTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt    >  1:3919367/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   <  1:1082646/93‑1 (MQ=255)
                                                                                       ggTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt  <  1:2020281/93‑1 (MQ=255)
                                                                                       ggTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt  >  2:1377995/1‑93 (MQ=255)
                                                                                       ggTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGCTAACCTTTACGTACGATCAtttttt  >  2:3677882/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/477902‑478081

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