Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,222,734 T→C intergenic (‑583/‑184) b1172 ← / → b1173 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,7340TC99.2% 400.5 / ‑4.8 129intergenic (‑583/‑184)b1172/b1173hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (67/61);  total (67/62)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 143 total aligned reads displayed.

TTATGCGGGATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGTCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGTCG  >  NC_000913/1222644‑1222823
                                                                                          |                                                                                         
ttATGCGGGATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCctggcct                                                                                         >  1:4065029/1‑93 (MQ=255)
 tATGCGGGATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCctggcctg                                                                                        <  2:86543/93‑1 (MQ=255)
 tATGCGGGATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCctggcctg                                                                                        >  1:2548327/1‑93 (MQ=255)
      gggATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGcc                                                                                   >  1:1462425/1‑93 (MQ=255)
       ggATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCt                                                                                  <  2:3732703/93‑1 (MQ=255)
       ggATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCt                                                                                  >  1:508133/1‑93 (MQ=255)
       ggATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCt                                                                                  <  2:2174881/93‑1 (MQ=255)
       ggATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCt                                                                                  >  2:3582904/1‑93 (MQ=255)
        gATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTg                                                                                 >  2:1769715/1‑93 (MQ=255)
          tatCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCg                                                                               >  2:4324137/1‑93 (MQ=255)
          tatCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCg                                                                               >  1:3678060/1‑93 (MQ=255)
           atcatcATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCgg                                                                              <  1:3928526/93‑1 (MQ=255)
             catcgtTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGc                                                                            >  2:1717232/1‑93 (MQ=255)
             catcatTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGc                                                                            <  2:3804600/93‑1 (MQ=255)
             catcatTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGc                                                                            <  2:424002/93‑1 (MQ=255)
             catcatTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGc                                                                            >  1:293761/1‑93 (MQ=255)
              atcatTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCt                                                                           >  2:4164073/1‑93 (MQ=255)
              atcatTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCt                                                                           >  1:1382375/1‑93 (MQ=255)
               tcatTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTg                                                                          >  2:4185878/1‑93 (MQ=255)
                catTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGa                                                                         <  2:4071354/93‑1 (MQ=255)
                 atTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGAt                                                                        >  2:520213/1‑93 (MQ=255)
                 atTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGAt                                                                        >  2:2274791/1‑93 (MQ=255)
                   ttCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGCCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGt                                                                      <  1:445973/93‑1 (MQ=255)
                    tcatcaTGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGtt                                                                     >  1:3702065/1‑93 (MQ=255)
                       tcatGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTgca                                                                  >  2:3362662/1‑93 (MQ=255)
                       tcatGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTgca                                                                  <  2:781140/93‑1 (MQ=255)
                       tcatGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTgca                                                                  >  1:2209540/1‑93 (MQ=255)
                         atgatgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTgcagc                                                                >  2:4564849/1‑93 (MQ=255)
                          tgatgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTgcagca                                                               <  2:3858972/93‑1 (MQ=255)
                           gatgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGcagcag                                                              <  1:516957/93‑1 (MQ=255)
                           gatgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGcagcag                                                              >  1:3405143/1‑93 (MQ=255)
                            atgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCaga                                                             >  1:1673359/1‑93 (MQ=255)
                            atgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCaga                                                             <  2:307171/93‑1 (MQ=255)
                            atgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCaga                                                             >  1:321404/1‑93 (MQ=255)
                            atgTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCaga                                                             <  2:2671153/93‑1 (MQ=255)
                              gTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGc                                                           <  2:1992444/93‑1 (MQ=255)
                              gTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGc                                                           <  1:4188343/93‑1 (MQ=255)
                              gTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGc                                                           >  1:3184516/1‑93 (MQ=255)
                              gTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGc                                                           >  2:4682080/1‑93 (MQ=255)
                               tCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCt                                                          <  2:4479532/93‑1 (MQ=255)
                               tCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCt                                                          >  1:1251784/1‑93 (MQ=255)
                                cTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCtt                                                         >  1:1561228/1‑93 (MQ=255)
                                 tttGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAGTCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTa                                                        <  1:1909016/93‑1 (MQ=255)
                                 tttGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTa                                                        >  1:3392340/1‑93 (MQ=255)
                                 tttGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTa                                                        >  1:751019/1‑93 (MQ=255)
                                  ttGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAg                                                       >  1:1091445/1‑93 (MQ=255)
                                    ggTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGt                                                     >  2:2590953/1‑93 (MQ=255)
                                    ggTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGt                                                     >  1:3766095/1‑93 (MQ=255)
                                    ggTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGt                                                     <  2:4595173/93‑1 (MQ=255)
                                     gTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTc                                                    <  2:1500216/93‑1 (MQ=255)
                                       gAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGGGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCat                                                  >  1:450206/1‑92 (MQ=255)
                                        aGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACt                                                 <  2:3399990/93‑1 (MQ=255)
                                          cGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTa                                               <  2:673434/93‑1 (MQ=255)
                                            gTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGa                                             <  1:1901972/93‑1 (MQ=255)
                                               aaCACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGc                                          >  2:3890838/1‑93 (MQ=255)
                                                   caATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGctact                                      <  2:1500736/93‑1 (MQ=255)
                                                   caATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGctact                                      >  1:401494/1‑93 (MQ=255)
                                                    aaTACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGctacta                                     >  2:3246163/1‑93 (MQ=255)
                                                     aTACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTAt                                    <  2:3613088/93‑1 (MQ=255)
                                                     aTACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTAt                                    <  2:3581381/93‑1 (MQ=255)
                                                       acacCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCa                                  <  1:2767235/93‑1 (MQ=255)
                                                       acacCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCa                                  <  1:4111843/93‑1 (MQ=255)
                                                        cacCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAt                                 >  1:1373425/1‑93 (MQ=255)
                                                        cacCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAt                                 >  1:1787338/1‑93 (MQ=255)
                                                        cacCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAt                                 <  2:3340624/93‑1 (MQ=255)
                                                         cccTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAtt                                <  1:2537478/92‑1 (MQ=255)
                                                         acCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAtt                                >  1:3224840/1‑93 (MQ=255)
                                                         acCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAtt                                <  2:2781627/93‑1 (MQ=255)
                                                         acCTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCAtt                                <  1:2033301/93‑1 (MQ=255)
                                                           cTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTa                              <  2:93161/93‑1 (MQ=255)
                                                           cTGCGCTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTa                              >  1:1688968/1‑93 (MQ=255)
                                                              cgcTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATcc                           <  2:3951221/93‑1 (MQ=255)
                                                               gcTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCa                          >  1:4245910/1‑93 (MQ=255)
                                                                cTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCaa                         <  2:443434/93‑1 (MQ=255)
                                                                cTGTCTCTTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCaa                         <  1:1013808/93‑1 (MQ=255)
                                                                    ctctTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGc                     <  1:2843066/93‑1 (MQ=255)
                                                                     tctTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCa                    >  1:1329980/1‑93 (MQ=255)
                                                                     tctTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCa                    >  1:1203233/1‑93 (MQ=255)
                                                                     tctTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCa                    >  2:673994/1‑93 (MQ=255)
                                                                     tctTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCa                    >  1:2227220/1‑93 (MQ=255)
                                                                      ctTCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCAt                   <  1:1246769/93‑1 (MQ=255)
                                                                       ttCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATg                  >  2:932927/1‑93 (MQ=255)
                                                                       ttCAGGAGGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATg                  <  2:2745677/93‑1 (MQ=255)
                                                                        tCAGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATgg                 >  2:3266263/1‑93 (MQ=255)
                                                                          aGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGaa               >  2:3664906/1‑93 (MQ=255)
                                                                          aGGATGAATCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGaa               <  1:1095376/93‑1 (MQ=255)
                                                                                aaTCCCCTGGTCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGaaa         <  1:4620100/93‑1 (MQ=255)
                                                                                aaTCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGaaa         <  1:3677009/93‑1 (MQ=255)
                                                                                aaTCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGaaa         >  2:659098/1‑93 (MQ=255)
                                                                                 aTCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAAt        >  2:3010556/1‑93 (MQ=255)
                                                                                  tCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATc       <  1:2006679/93‑1 (MQ=255)
                                                                                  tCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATc       <  1:2974541/93‑1 (MQ=255)
                                                                                  tCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATc       <  2:1859268/93‑1 (MQ=255)
                                                                                  tCCCCTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATc       <  1:1183167/93‑1 (MQ=255)
                                                                                    cccTGGCCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAg     <  2:1075345/93‑1 (MQ=255)
                                                                                     cctggcctggTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGt    >  1:2477731/1‑93 (MQ=255)
                                                                                      ctggcctggtGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGTc   >  1:431498/1‑93 (MQ=255)
                                                                                       tggcctggtGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGTCg  >  2:1086149/1‑93 (MQ=255)
                                                                                       tggcctggtGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGTCg  <  1:4663457/93‑1 (MQ=255)
                                                                                       tggcctggtGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGTCg  <  1:3777227/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTATGCGGGATATCATCATTTCATCATGATGTCTTTGGTGAGCGGTGAACACAATACACCTGCGCTGTCTCTTCAGGATGAATCCCCTGGTCTGGTGCCTGCGGGCTGATGTTGCAGCAGAGCTTAGGTCACTTAGACGCTACTATCATTTATCCAATGGCATGGAATCGAAATCAGTCG  >  NC_000913/1222644‑1222823

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