Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 657,402 G→A Q37Q (CAG→CAA cspE → transcription antiterminator and regulator of RNA stability

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,4020GA79.1% 785.8 / 197.1 306Q37Q (CAG→CAAcspEtranscription antiterminator and regulator of RNA stability
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (106/136);  total (138/168)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 307 total aligned reads displayed.

TTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAA  >  NC_000913/657302‑657501
                                                                                                    |                                                                                                   
  aaaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAAc                                                                                                   >  2:1084853/1‑101 (MQ=255)
  aaaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAAc                                                                                                   <  1:1084451/101‑1 (MQ=255)
  aaaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAAc                                                                                                   <  2:1084284/101‑1 (MQ=255)
   aaGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACt                                                                                                  >  1:1084856/1‑101 (MQ=255)
       tAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATg                                                                                              >  2:1084859/1‑101 (MQ=255)
       tAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATg                                                                                              <  1:1084293/101‑1 (MQ=255)
        aaCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATgg                                                                                             >  2:1084864/1‑101 (MQ=255)
          cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATg                                                                                              <  1:1084538/98‑1 (MQ=255)
          cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGtt                                                                                           <  1:1084448/101‑1 (MQ=255)
           gTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                                                                          <  1:1084394/101‑1 (MQ=255)
           gTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGttt                                                                                          <  1:1084416/101‑1 (MQ=255)
           gTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGttt                                                                                          >  2:1084867/1‑101 (MQ=255)
            ttAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGtttt                                                                                         <  1:1084513/101‑1 (MQ=255)
             tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTa                                                                                        >  2:1084870/1‑101 (MQ=255)
                   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAActc                                                                                  <  1:1084482/101‑1 (MQ=255)
                   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAActc                                                                                  <  1:1084348/101‑1 (MQ=255)
                   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAActc                                                                                  <  2:1084422/101‑1 (MQ=255)
                    tttAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActct                                                                                 >  2:1084878/1‑101 (MQ=255)
                    tttAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAActct                                                                                 <  1:1083890/101‑1 (MQ=255)
                    tttAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAActct                                                                                 >  1:1084876/1‑101 (MQ=255)
                    tttAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAActct                                                                                 <  2:1084701/101‑1 (MQ=255)
                     ttAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCtt                                                                                <  1:1084428/101‑1 (MQ=255)
                       aaTGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGc                                                                              >  2:1084881/1‑101 (MQ=255)
                       aaTGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGc                                                                              <  1:1084249/101‑1 (MQ=255)
                         tGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTg                                                                            <  2:1084329/101‑1 (MQ=255)
                         tGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTg                                                                            >  1:1084884/1‑101 (MQ=255)
                          gAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                                                                           <  1:1084403/101‑1 (MQ=255)
                          gAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGa                                                                           <  1:1084198/101‑1 (MQ=255)
                            gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAg                                                                         <  1:1084500/101‑1 (MQ=255)
                            gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAg                                                                         >  1:1084888/1‑101 (MQ=255)
                               cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTg                                                                               >  2:1084896/1‑92 (MQ=255)
                               cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                      <  1:1084565/101‑1 (MQ=255)
                               cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                      >  1:1084895/1‑101 (MQ=255)
                               cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTc                                                                      >  2:1084892/1‑101 (MQ=255)
                                aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCa                                                                     <  2:1084455/101‑1 (MQ=255)
                                aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCa                                                                     <  1:1084386/101‑1 (MQ=255)
                                  aGGATTCGGTTTCATTTCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCtt                                                                                <  1:1084436/88‑1 (MQ=255)
                                  aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc                                                                   <  1:1084367/101‑1 (MQ=255)
                                  aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc                                                                   <  1:1084868/101‑1 (MQ=255)
                                    gATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcgc                                                                 >  1:1084901/1‑101 (MQ=255)
                                     aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGcgcg                                                                >  1:1084905/1‑101 (MQ=255)
                                     aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGcgcg                                                                >  2:1084904/1‑101 (MQ=255)
                                       tCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTa                                                              <  1:1084705/101‑1 (MQ=255)
                                         ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTaga                                                            <  1:1084407/101‑1 (MQ=255)
                                          gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTagag                                                           <  1:1084624/101‑1 (MQ=255)
                                          gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTagag                                                           <  1:1084465/101‑1 (MQ=255)
                                           tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGt                                                          >  1:1084911/1‑101 (MQ=255)
                                           tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGt                                                          >  2:1084909/1‑101 (MQ=255)
                                            ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGtt                                                         <  1:1084543/101‑1 (MQ=255)
                                             tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTc                                                        <  1:1084279/101‑1 (MQ=255)
                                             tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTc                                                        <  2:1084526/101‑1 (MQ=255)
                                              cATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCg                                                       >  2:1084915/1‑101 (MQ=255)
                                               aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa                                                      <  1:1084304/101‑1 (MQ=255)
                                                ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaa                                                     >  1:1084917/1‑101 (MQ=255)
                                                ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaa                                                     <  2:1084668/101‑1 (MQ=255)
                                                  aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAAt                                                   >  2:1084923/1‑101 (MQ=255)
                                                   cTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATc                                                  >  1:1084924/1‑101 (MQ=255)
                                                    tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCa                                                 <  1:1084590/101‑1 (MQ=255)
                                                    tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCa                                                 <  2:1084444/101‑1 (MQ=255)
                                                     ccGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAc                                                >  2:1084928/1‑101 (MQ=255)
                                                         aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAAc                                            >  1:1084932/1‑101 (MQ=255)
                                                          aGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACg                                           >  2:1084934/1‑101 (MQ=255)
                                                            aCGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATggct                                                                                           >  1:1084937/1‑49 (MQ=255)
                                                            aCGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGt                                         <  1:1084462/101‑1 (MQ=255)
                                                              ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTg                                        <  1:1084441/100‑1 (MQ=255)
                                                              ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGc                                       >  2:1084940/1‑101 (MQ=255)
                                                              ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGc                                       >  1:1084943/1‑101 (MQ=255)
                                                               gcagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGcc                                      >  2:1084945/1‑101 (MQ=255)
                                                                 agcaAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCaa                                    >  2:1084948/1‑101 (MQ=255)
                                                                   caAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAg                                  >  1:1084953/1‑101 (MQ=255)
                                                                   caAAGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAg                                  <  2:1084552/101‑1 (MQ=255)
                                                                     aaGACGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGc                                >  2:1084956/1‑101 (MQ=255)
                                                                       gACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAg                                  >  2:1084959/1‑97 (MQ=255)
                                                                        aCGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCt                             <  2:1084160/101‑1 (MQ=255)
                                                                        aCGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCt                             >  1:1084960/1‑101 (MQ=255)
                                                                         cGTGTTCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCtt                            <  2:1084515/101‑1 (MQ=255)
                                                                          gtgtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTc                           <  1:1084663/101‑1 (MQ=255)
                                                                           tgtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTct                          <  2:1084402/101‑1 (MQ=255)
                                                                              tCGTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTctgct                       >  2:1084965/1‑101 (MQ=255)
                                                                                gTACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCtgctgc                     >  1:1084968/1‑101 (MQ=255)
                                                                                 tACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCa                    <  1:1084652/101‑1 (MQ=255)
                                                                                 tACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCa                    <  1:1084631/101‑1 (MQ=255)
                                                                                 tACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCa                    <  2:1084571/101‑1 (MQ=255)
                                                                                 tACACTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCa                    >  2:1084971/1‑101 (MQ=255)
                                                                                   cacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCaaa                  >  1:1084972/1‑101 (MQ=255)
                                                                                     cTTCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACg                <  2:1084483/101‑1 (MQ=255)
                                                                                      ttCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGt               <  1:1084424/101‑1 (MQ=255)
                                                                                      ttCTCTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGt               >  2:1084976/1‑101 (MQ=255)
                                                                                       tctcTGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTa              <  2:1084595/101‑1 (MQ=255)
                                                                                         tctGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAAt            <  1:1084520/101‑1 (MQ=255)
                                                                                           tGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAg                                                                         >  2:1084981/1‑38 (MQ=255)
                                                                                           tGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCg          <  1:1084690/101‑1 (MQ=255)
                                                                                           tGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCg          >  1:1084978/1‑101 (MQ=255)
                                                                                           tGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCg          <  2:1084592/101‑1 (MQ=255)
                                                                                           tGCAATCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCg          <  2:1084505/101‑1 (MQ=255)
                                                                                             cAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGct        >  2:1084984/1‑101 (MQ=255)
                                                                                               aTCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGctct      >  2:1084987/1‑101 (MQ=255)
                                                                                               aTCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGctat      <  1:1084478/101‑3 (MQ=255)
                                                                                                tCCAAACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTg     <  1:1084474/101‑1 (MQ=255)
                                                                                                  cAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTa   <  1:1084836/101‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTAAAGGTAACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCACTAACGGTGCCAAAGGCCCTTCTGCTGCAAACGTAATCGCTCTGTAA  >  NC_000913/657302‑657501

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: