Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,396,030 G→A intergenic (‑108/‑46) ynaI ← / → insH‑4 small conductance mechanosensitive channel YnaI/IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,396,0300GA92.1% 1246.9 / 69.4 340intergenic (‑108/‑46)ynaI/insH‑4small conductance mechanosensitive channel YnaI/IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (11/16);  new base A (143/170);  total (154/186)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

Only 100 of 345 total aligned reads displayed.

TCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTA  >  NC_000913/1395930‑1396127
                                                                                                    |                                                                                                 
tCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAg                                                                                                   <  2:2552561/101‑1 (MQ=255)
   ggCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACaa                                                                                                <  1:2552347/101‑1 (MQ=255)
     ctcctcCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACa                                                                                              >  1:2553232/1‑101 (MQ=255)
      tcctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACAt                                                                                             <  2:2552482/101‑1 (MQ=255)
       cctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACAtt                                                                                            <  1:2552418/101‑1 (MQ=255)
       cctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACAtt                                                                                            <  2:2552687/101‑1 (MQ=255)
        ctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTa                                                                                           <  1:2552759/101‑1 (MQ=255)
        ctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTa                                                                                           >  2:2553240/1‑101 (MQ=255)
        ctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTa                                                                                           <  1:2552590/101‑1 (MQ=255)
         tccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAg                                                                                          <  2:2552743/101‑1 (MQ=255)
          ccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAg                                                                                          <  2:2552636/100‑1 (MQ=255)
          ccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGc                                                                                         <  2:2552723/101‑1 (MQ=255)
          ccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGc                                                                                         <  2:2552652/101‑1 (MQ=255)
            aGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTa                                                                                       >  1:2553242/1‑101 (MQ=255)
            aGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTa                                                                                       <  1:2552354/101‑1 (MQ=255)
            aGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTa                                                                                       <  1:2552986/101‑1 (MQ=255)
            aGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTa                                                                                       <  2:2552611/101‑1 (MQ=255)
             gATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaa                                                                                      <  1:2553049/101‑1 (MQ=255)
              aTAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaag                                                                                     >  1:2553247/1‑101 (MQ=255)
                aaTTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaagga                                                                                   >  1:2553252/1‑101 (MQ=255)
                aaTTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaagga                                                                                   >  2:2553249/1‑101 (MQ=255)
                 aTTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaaggaa                                                                                  <  1:2552507/101‑1 (MQ=255)
                  ttGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaaggaag                                                                                 <  2:2552668/101‑1 (MQ=255)
                   tGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTaaggaagg                                                                                <  1:2552839/101‑1 (MQ=255)
                    gtcgtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGt                                                                               <  1:2552639/101‑1 (MQ=255)
                    gtcgtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGt                                                                               <  1:2552458/101‑1 (MQ=255)
                    gtcgtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGt                                                                               >  1:2553255/1‑101 (MQ=255)
                      cgtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGc                                                                             <  2:2552818/101‑1 (MQ=255)
                          aTGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAAc                                                                         <  2:2552288/101‑1 (MQ=255)
                          aTGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAAc                                                                         >  2:2553259/1‑101 (MQ=255)
                            gCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACaa                                                                       <  2:2552951/101‑1 (MQ=255)
                               ggAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTc                                                                    >  1:2553263/1‑101 (MQ=255)
                                gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGt                                                                               <  2:2552674/89‑1 (MQ=255)
                                gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTcc                                                                   <  1:2552328/101‑1 (MQ=255)
                                gAAACCCCTGGCAGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTcc                                                                   >  2:2553265/1‑101 (MQ=255)
                                  aaCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCt                                                                 <  2:2552756/101‑1 (MQ=255)
                                  aaCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCt                                                                 <  1:2552575/101‑1 (MQ=255)
                                  aaCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCt                                                                 <  2:2552909/101‑1 (MQ=255)
                                   aCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTg                                                                <  2:2552888/101‑1 (MQ=255)
                                   aCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTg                                                                <  1:2552678/101‑1 (MQ=255)
                                   aCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTg                                                                <  2:2552444/101‑1 (MQ=255)
                                      ccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGata                                                             >  1:2553271/1‑101 (MQ=255)
                                      ccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGata                                                             <  2:2552828/101‑1 (MQ=255)
                                       cTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGatat                                                            >  2:2553273/1‑101 (MQ=255)
                                       cTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGatat                                                            <  2:2552791/101‑1 (MQ=255)
                                       cTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGatat                                                            <  2:2552515/101‑1 (MQ=255)
                                          gCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATgag                                                         >  1:2553276/1‑101 (MQ=255)
                                           cGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATgaga                                                        <  1:2552618/101‑1 (MQ=255)
                                           cGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATgaga                                                        <  2:2552991/101‑1 (MQ=255)
                                            ggggCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGAt                                                       <  1:2552707/101‑1 (MQ=255)
                                            ggggCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGAt                                                       >  2:2553282/1‑101 (MQ=255)
                                             gggCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATc                                                      >  1:2553284/1‑101 (MQ=255)
                                             gggCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATc                                                      <  1:2552716/101‑1 (MQ=255)
                                               gCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCAt                                                    >  2:2553286/1‑101 (MQ=255)
                                                 tATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGt                                                  <  1:2552651/101‑1 (MQ=255)
                                                  aTTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGtt                                                 >  1:2553291/1‑101 (MQ=255)
                                                    tttACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTg                                               >  1:2553297/1‑101 (MQ=255)
                                                    tttACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTg                                               >  2:2553296/1‑101 (MQ=255)
                                                    tttACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTg                                               >  1:2553301/1‑101 (MQ=255)
                                                      tACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTc                                             <  2:2553054/101‑1 (MQ=255)
                                                      tACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTc                                             >  2:2553303/1‑101 (MQ=255)
                                                       aCCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCa                                            >  2:2553311/1‑101 (MQ=255)
                                                       aCCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCa                                            >  1:2553309/1‑101 (MQ=255)
                                                        ccGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAt                                           <  1:2552907/101‑1 (MQ=255)
                                                          gcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCt                                         <  1:2552950/101‑1 (MQ=255)
                                                          gcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCt                                         <  2:2553019/101‑1 (MQ=255)
                                                          gcgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCt                                         >  2:2553315/1‑101 (MQ=255)
                                                           cgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTg                                        <  1:2553010/101‑1 (MQ=255)
                                                           cgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTg                                        >  1:2553318/1‑101 (MQ=255)
                                                            gACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTgg                                       <  2:2552735/101‑1 (MQ=255)
                                                            gACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTgg                                       >  2:2553321/1‑101 (MQ=255)
                                                             aCAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGa                                      <  1:2552815/101‑1 (MQ=255)
                                                               aaTTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGc                                    >  2:2553328/1‑101 (MQ=255)
                                                               aaTTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGc                                    >  1:2553327/1‑101 (MQ=255)
                                                                  tcattcaGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCAt                                 >  1:2553332/1‑101 (MQ=255)
                                                                      tcaGATCATCAATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGaa                             >  2:2553335/1‑101 (MQ=255)
                                                                         gATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAg                          >  1:2553338/1‑101 (MQ=255)
                                                                            catcaATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGt                       <  2:2552865/101‑1 (MQ=255)
                                                                             atcaATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                      >  2:2553345/1‑101 (MQ=255)
                                                                             atcaATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                      >  1:2553344/1‑101 (MQ=255)
                                                                              tcaATAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGGAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtc                     >  2:2553348/1‑101 (MQ=255)
                                                                                aaTAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                   >  1:2553351/1‑101 (MQ=255)
                                                                                  tAGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcatca                 <  1:2552858/101‑1 (MQ=255)
                                                                                   aGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                >  1:2553356/1‑101 (MQ=255)
                                                                                   aGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                <  1:2552749/101‑1 (MQ=255)
                                                                                   aGTCAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                <  1:2552873/101‑1 (MQ=255)
                                                                                    gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATg               <  2:2552810/101‑1 (MQ=255)
                                                                                    gTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATg               >  1:2553360/1‑101 (MQ=255)
                                                                                      cAGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAg             >  2:2553364/1‑101 (MQ=255)
                                                                                       aGGGAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGt            <  2:2552764/101‑1 (MQ=255)
                                                                                        gggAAGGAAGTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtc           >  1:2553367/1‑101 (MQ=255)
                                                                                          gaaggaagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcat         >  2:2553370/1‑101 (MQ=255)
                                                                                           aaggaagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatc        <  1:2552786/101‑1 (MQ=255)
                                                                                            aggaagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTccc                                                                  >  1:2553374/1‑42 (MQ=255)
                                                                                            aggaagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatca       <  2:2552776/101‑1 (MQ=255)
                                                                                             ggaagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCaa      >  2:2553376/1‑101 (MQ=255)
                                                                                              gaagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAAc     <  2:2552393/101‑1 (MQ=255)
                                                                                               aagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACt    <  2:2553213/101‑1 (MQ=255)
                                                                                               aagTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACt    >  1:2553381/1‑101 (MQ=255)
                                                                                                agTAACAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACtt   >  2:2553383/1‑101 (MQ=255)
                                                                                                    |                                                                                                 
TCAGGCTCCTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTA  >  NC_000913/1395930‑1396127

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