Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,181,376 G→A A852T (GCT→ACT)  rpoB → RNA polymerase beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,181,3760GA85.6% 784.9 / 92.1 278A852T (GCT→ACT) rpoBRNA polymerase beta subunit
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base A (106/132);  total (124/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-01

Only 100 of 292 total aligned reads displayed.

GTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATT  >  NC_000913/4181284‑4181468
                                                                                            |                                                                                            
   ttgttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTg                                                                                           >  1:1112921/1‑93 (MQ=255)
     gttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGAc                                                                                         <  1:1239680/93‑1 (MQ=255)
        cAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATc                                                                                      <  1:586986/93‑1 (MQ=255)
        cAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATc                                                                                      >  2:6257224/1‑93 (MQ=255)
         aGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGTCACCAAGCTGGGTTCGGAAGAGATCACGCCTGACATcc                                                                                     >  1:5380300/1‑93 (MQ=255)
         aGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATcc                                                                                     >  2:6563833/1‑93 (MQ=255)
         aGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATcc                                                                                     <  2:6113599/93‑1 (MQ=255)
         aGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATcc                                                                                     >  2:3455251/1‑93 (MQ=255)
            aaGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGa                                                                                  <  2:4835943/93‑1 (MQ=255)
             acaCCGTTTCACAACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGaa                                                                                 <  2:646345/91‑1 (MQ=255)
              gACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACGAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAAc                                                                                <  2:4965955/93‑1 (MQ=255)
              gACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAAc                                                                                >  1:5918081/1‑93 (MQ=255)
              aaCCGTTACACCACCATCCACATTCAGGAACGGGCGTTTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAAc                                                                                <  1:4463320/92‑1 (MQ=255)
               aCCTTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACg                                                                               <  2:4736307/93‑1 (MQ=255)
               aCCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACg                                                                               <  1:922488/93‑1 (MQ=255)
                ccGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGt                                                                              >  1:2819831/1‑93 (MQ=255)
                 cGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTg                                                                             >  2:4638404/1‑93 (MQ=255)
                 cGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTg                                                                             <  2:4054385/93‑1 (MQ=255)
                  gTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTgg                                                                            >  2:3869753/1‑93 (MQ=255)
                   tttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTTCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTggg                                                                           >  1:2040930/1‑93 (MQ=255)
                   tttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTggg                                                                           <  1:5080702/93‑1 (MQ=255)
                    ttCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGt                                                                          <  2:6078485/93‑1 (MQ=255)
                     tCACCCCCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTg                                                                         <  1:1654325/93‑1 (MQ=255)
                     tCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTg                                                                         <  2:2958622/93‑1 (MQ=255)
                     tCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTg                                                                         <  2:1947089/93‑1 (MQ=255)
                       accaccATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGaa                                                                       >  2:2591442/1‑93 (MQ=255)
                       accaccATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGaa                                                                       <  2:4236514/93‑1 (MQ=255)
                       accaccATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGACCGGAATAGATGACCACTGACATCCCGAACGTGGGTgga                                                                       >  1:2764461/1‑91 (MQ=255)
                        ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAg                                                                      >  2:1234635/1‑93 (MQ=255)
                        ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAg                                                                      >  2:517391/1‑93 (MQ=255)
                        ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAg                                                                      <  2:5109881/93‑1 (MQ=255)
                        ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAg                                                                      <  2:805800/93‑1 (MQ=255)
                         caccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGc                                                                     <  1:3171156/93‑1 (MQ=255)
                          cccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCt                                                                    <  1:2355525/92‑1 (MQ=255)
                          accaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCt                                                                    <  1:716772/93‑1 (MQ=255)
                           ccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTg                                                                   >  1:2702931/1‑93 (MQ=255)
                            caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGTAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTgc                                                                  <  2:4792916/93‑1 (MQ=255)
                            caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCCCTGACATCCCGCACGTGGGTGACGCTgg                                                                  >  1:552756/1‑92 (MQ=255)
                            caTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTgc                                                                  <  1:5738434/93‑1 (MQ=255)
                               ccACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGct                                                               <  1:3368952/93‑1 (MQ=255)
                                 acaTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGctct                                                             <  2:6348097/93‑1 (MQ=255)
                                  caTTCAGGAACTGGAGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGCAGATATCACCACTGAGATAGCGAACGTGGGGGGAGCCGCGCgctc                                                            >  1:2957954/1‑93 (MQ=255)
                                   aTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTcc                                                           <  2:18233/93‑1 (MQ=255)
                                   aTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTcc                                                           >  1:3743861/1‑93 (MQ=255)
                                    ttCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCa                                                          <  1:2201405/93‑1 (MQ=255)
                                    ttCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCa                                                          >  2:3628327/1‑93 (MQ=255)
                                    ttCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCa                                                          >  2:2894333/1‑93 (MQ=255)
                                     tCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCaa                                                         >  1:3600321/1‑93 (MQ=255)
                                      cAGGAACTGGCGTGTGTGTTCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCaaa                                                        >  1:6496045/1‑93 (MQ=255)
                                      cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCaaa                                                        >  2:3754238/1‑93 (MQ=255)
                                       aGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                       <  2:5785948/93‑1 (MQ=255)
                                       aGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAAc                                                       >  1:4804491/1‑93 (MQ=255)
                                           aCTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGa                                                   <  2:1686438/93‑1 (MQ=255)
                                           aCTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGa                                                   <  2:1780862/93‑1 (MQ=255)
                                            cTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAt                                                  <  2:2445394/93‑1 (MQ=255)
                                            cTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAt                                                  <  2:5758441/93‑1 (MQ=255)
                                            cTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAt                                                  >  2:5353822/1‑93 (MQ=255)
                                            cTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAt                                                  <  2:4147844/93‑1 (MQ=255)
                                               gCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGaa                                               >  1:3234385/1‑93 (MQ=255)
                                                cGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAAt                                              >  1:4917420/1‑93 (MQ=255)
                                                 gtgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATc                                             <  1:1771994/93‑1 (MQ=255)
                                                  tgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATcc                                            <  2:5275258/93‑1 (MQ=255)
                                                  tgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATcc                                            >  1:408736/1‑93 (MQ=255)
                                                  tgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATcc                                            <  1:814472/93‑1 (MQ=255)
                                                  tgtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATcc                                            >  1:3770512/1‑93 (MQ=255)
                                                     gtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                         >  2:1060062/1‑93 (MQ=255)
                                                     gtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGt                                         >  1:5283285/1‑93 (MQ=255)
                                                      tgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTa                                        <  1:2622644/93‑1 (MQ=255)
                                                        tCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATc                                      >  2:1035079/1‑93 (MQ=255)
                                                          ccGTGACACCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGt                                    <  1:4702895/93‑1 (MQ=255)
                                                               acacCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACa                               <  2:4494650/93‑1 (MQ=255)
                                                               acacCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACa                               >  2:3678397/1‑93 (MQ=255)
                                                                cacCAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAt                              >  1:3102264/1‑93 (MQ=255)
                                                                   cAAGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTgg                           <  1:5977281/93‑1 (MQ=255)
                                                                    aaGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGt                          >  1:6394005/1‑93 (MQ=255)
                                                                    aaGCGGGATCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGt                          <  1:1377803/93‑1 (MQ=255)
                                                                     aGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTg                         <  2:2288089/93‑1 (MQ=255)
                                                                     aGCTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTg                         <  1:2384858/93‑1 (MQ=255)
                                                                       cTGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCg                       <  2:4166173/93‑1 (MQ=255)
                                                                        tGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCgg                      >  1:5211925/1‑93 (MQ=255)
                                                                        tGGGTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCgg                      <  1:5609407/93‑1 (MQ=255)
                                                                          ggTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGaa                    >  1:1483706/1‑93 (MQ=255)
                                                                          ggTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGaa                    <  1:4239658/93‑1 (MQ=255)
                                                                           gTCCGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAg                   <  1:4310740/93‑1 (MQ=255)
                                                                              cGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGa                >  2:6449645/1‑93 (MQ=255)
                                                                              cGGAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAACCTGGATGAATCCGGTATCGTTTACATTGGTGCGGACGTGa                >  2:1718390/1‑93 (MQ=255)
                                                                               ggAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAc               >  2:1435844/1‑93 (MQ=255)
                                                                               ggAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAc               <  1:4103036/93‑1 (MQ=255)
                                                                                gAAGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGAcc              <  2:3098645/93‑1 (MQ=255)
                                                                                 aaGAGATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg             <  1:4868988/93‑1 (MQ=255)
                                                                                  agagATCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg            >  2:5548992/1‑93 (MQ=255)
                                                                                   gagaTCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCggg           >  2:2162847/1‑92 (MQ=255)
                                                                                   gagaTCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGt           <  1:5498237/93‑1 (MQ=255)
                                                                                   gagaTCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGt           <  2:1512456/93‑1 (MQ=255)
                                                                                       tCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg       <  1:6132538/93‑1 (MQ=255)
                                                                                       tCACCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCg       <  2:160920/93‑1 (MQ=255)
                                                                                        caccacTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGa      <  2:2260089/93‑1 (MQ=255)
                                                                                        cACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGa      >  1:6065873/1‑93 (MQ=255)
                                                                                         aCCACTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc     <  1:3914632/93‑1 (MQ=255)
                                                                                            gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt  <  2:4005903/93‑1 (MQ=255)
                                                                                            |                                                                                            
GTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATT  >  NC_000913/4181284‑4181468

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