Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,027,751 (C)6→7 coding (191/921 nt) yedA → putative transmembrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,027,7451.C98.2% 571.1 / 1.9 167coding (185/921 nt)yedAputative transmembrane subunit
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (86/78);  total (86/81)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

GGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACA  >  NC_000913/2027653‑2027832
                                                                                               |                                                                                        
ggAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGGTACTGCGCGGACACCAACt                                                                                           >  1:5571576/1‑93 (MQ=255)
  aaaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                        <  2:1249573/93‑1 (MQ=255)
  aaaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcc                                                                                        <  2:2980616/93‑1 (MQ=255)
   aaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                      <  2:894872/93‑1 (MQ=255)
   aaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                      >  1:5464306/1‑93 (MQ=255)
   aaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                      <  2:2915132/93‑1 (MQ=255)
    aGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccc                                                                                     >  1:116666/1‑93 (MQ=255)
     gCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccccc                                                                                    >  2:3586260/1‑93 (MQ=255)
      cTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccc                                                                                   <  2:589497/93‑1 (MQ=255)
        ggCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccccg                                                                                 >  2:1377169/1‑91 (MQ=255)
        ggCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccccg                                                                                 >  2:4107355/1‑91 (MQ=255)
        ggCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccccg                                                                                 >  1:2307086/1‑91 (MQ=255)
        ggCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccccg                                                                                 >  1:2193675/1‑91 (MQ=255)
         gCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccccgc                                                                                >  2:2838923/1‑90 (MQ=255)
         gCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccccgc                                                                                >  2:2292002/1‑90 (MQ=255)
            tCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCCTTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTAc                                                                              >  1:6190022/1‑93 (MQ=255)
            tCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTAc                                                                              >  1:5895546/1‑93 (MQ=255)
            tCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTAc                                                                              <  2:3271014/93‑1 (MQ=255)
             ccGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACg                                                                             >  2:56852/1‑93 (MQ=255)
             ccGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACg                                                                             >  1:4005038/1‑93 (MQ=255)
               gTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTc                                                                           <  2:4207821/93‑1 (MQ=255)
               gTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTc                                                                           >  2:6352715/1‑93 (MQ=255)
               gTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTc                                                                           >  2:623653/1‑93 (MQ=255)
                ttAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTcc                                                                          >  2:2761400/1‑93 (MQ=255)
                ttAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTcc                                                                          <  1:2479298/93‑1 (MQ=255)
                ttAATGATGGCGGGAGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTcc                                                                          <  1:2112589/93‑1 (MQ=255)
                   atgatgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCgct                                                                       <  1:2511481/93‑1 (MQ=255)
                      atgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCgctgct                                                                    <  2:364125/93‑1 (MQ=255)
                       tgGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTc                                                                   <  2:877230/93‑1 (MQ=255)
                        ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCa                                                                  >  1:673717/1‑93 (MQ=255)
                        ggcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCa                                                                  >  1:6119948/1‑93 (MQ=255)
                         gcgggcgTTGTATTCCTGGCAGCCGGTATTTTAGTGCTGGCATTTTTGCTACTGCGCGGACACAAACGCCCCCCCGCTACGTCCGCTGCTCaa                                                                 <  2:816130/93‑1 (MQ=255)
                         gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCaa                                                                 <  2:3048830/93‑1 (MQ=255)
                          ggggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAACCTCCCCCCCGCTACGTCCGCTGCTCAAt                                                                <  1:6174152/92‑1 (MQ=255)
                            ggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAACCTCCCCCCCGCTACGTCCGCTGCTCAATGc                                                              <  1:4605704/93‑1 (MQ=255)
                            ggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGc                                                              >  2:2709841/1‑93 (MQ=255)
                            ggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGc                                                              >  2:1321532/1‑93 (MQ=255)
                             gcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGcc                                                             >  2:2882360/1‑93 (MQ=255)
                               gTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCcgc                                                           >  2:6249817/1‑93 (MQ=255)
                                  cGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCt                                                        >  1:542598/1‑93 (MQ=255)
                                      tCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGAtt                                                    <  1:1707404/93‑1 (MQ=255)
                                      tCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGAtt                                                    <  1:1634012/93‑1 (MQ=255)
                                          ggCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGcc                                                <  2:2008621/93‑1 (MQ=255)
                                           gCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCt                                               <  1:2420084/93‑1 (MQ=255)
                                           gCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCt                                               <  2:6321208/93‑1 (MQ=255)
                                           gCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCt                                               <  1:3177908/93‑1 (MQ=255)
                                             cgCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACCCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGt                                             <  1:4499634/92‑1 (MQ=255)
                                             aGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGt                                             >  2:1788705/1‑93 (MQ=255)
                                               ccGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTa                                           >  1:1996036/1‑93 (MQ=255)
                                                 ggTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTAtt                                         <  2:3997844/93‑1 (MQ=255)
                                                 ggTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTAtt                                         <  2:4900769/93‑1 (MQ=255)
                                                 ggTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTAtt                                         >  2:4984526/1‑93 (MQ=255)
                                                  gTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTg                                        >  1:5220856/1‑93 (MQ=255)
                                                    aTTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgct                                      <  1:1649970/93‑1 (MQ=255)
                                                    aTTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgct                                      <  1:2315523/93‑1 (MQ=255)
                                                     ttttATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctg                                     <  2:6299355/93‑1 (MQ=255)
                                                     ttttATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctg                                     <  1:2903860/93‑1 (MQ=255)
                                                     ttttATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctg                                     >  2:722477/1‑93 (MQ=255)
                                                         aTTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTgctggctg                                 >  2:2108752/1‑93 (MQ=255)
                                                          ttGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                <  2:2543774/93‑1 (MQ=255)
                                                          ttGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                <  1:4927139/93‑1 (MQ=255)
                                                          ttGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGt                                <  2:2550226/93‑1 (MQ=255)
                                                             cTGGCATTTTTGCGACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCgg                             <  1:3752177/93‑1 (MQ=255)
                                                              tGGTATTTTTGATACTGCGCGGACACAATCTCCACCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTAGGt                            <  1:5909902/93‑1 (MQ=255)
                                                              tGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGt                            >  2:3384596/1‑93 (MQ=255)
                                                                  aTTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATg                        <  2:5155651/93‑1 (MQ=255)
                                                                  aTTTTTGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATg                        <  1:4061761/93‑1 (MQ=255)
                                                                   tttttGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATgg                       >  2:4844314/1‑93 (MQ=255)
                                                                    ttttGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGc                      <  2:4811660/93‑1 (MQ=255)
                                                                     tttGCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCa                     <  2:69719/93‑1 (MQ=255)
                                                                        gCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                  >  2:4367294/1‑93 (MQ=255)
                                                                        gCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                  <  2:1571347/93‑1 (MQ=255)
                                                                        gCTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATgg                  <  2:3881684/93‑1 (MQ=255)
                                                                         cTACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGt                 >  2:1521683/1‑93 (MQ=255)
                                                                          tACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTg                <  1:3242764/93‑1 (MQ=255)
                                                                          tACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTg                >  1:4095065/1‑93 (MQ=255)
                                                                          tACTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTg                >  1:6078312/1‑93 (MQ=255)
                                                                           agtGCGCGGACACAAACTCCCCCCCTCTACGTCCGCTGCTCAATGCCGCGCGTATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGa               <  1:2653783/91‑1 (MQ=255)
                                                                            cTGCGCGGACACAATCTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGAc              <  2:6034300/93‑1 (MQ=255)
                                                                            cTGCGCGGACACAACCTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGAc              <  2:5706316/93‑1 (MQ=255)
                                                                            cTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGAc              >  1:1774970/1‑93 (MQ=255)
                                                                            cTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGAc              <  1:3372401/93‑1 (MQ=255)
                                                                            cTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTACTGCTGGCTGTCGGTAATGGCATGGTGAc              <  1:5713281/93‑1 (MQ=255)
                                                                            cTGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAAGCCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATggggcg              >  2:3189949/1‑89 (MQ=255)
                                                                             tGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACg             >  2:1168451/1‑93 (MQ=255)
                                                                             tGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACg             >  2:5748010/1‑93 (MQ=255)
                                                                             tGCGCGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACg             >  2:1915693/1‑93 (MQ=255)
                                                                              gcgcGGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACgg            <  2:3580070/93‑1 (MQ=255)
                                                                               cgcgGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGt           <  1:5395512/93‑1 (MQ=255)
                                                                                gcgGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGtt          >  1:6241452/1‑93 (MQ=255)
                                                                                gcgGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGtt          >  2:4117638/1‑93 (MQ=255)
                                                                                gcgGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGtt          <  1:3931154/93‑1 (MQ=255)
                                                                                ccgGACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGtt          <  1:1430057/92‑1 (MQ=255)
                                                                                  ggACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGc        >  2:5653157/1‑93 (MQ=255)
                                                                                   gACACAAACTCCCCCCCGCTACGTTCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGCGACGGTTGcc       >  2:3786656/1‑93 (MQ=255)
                                                                                   gACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGcc       >  2:1102377/1‑93 (MQ=255)
                                                                                   gACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGcc       >  1:6432005/1‑93 (MQ=255)
                                                                                   gACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGcc       <  1:1677762/93‑1 (MQ=255)
                                                                                   gACACAAACTCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGGTGcc       >  2:2232637/1‑93 (MQ=255)
                                                                                     cacaAACCCCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGa     <  2:4755142/93‑1 (MQ=255)
                                                                                               |                                                                                        
GGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACA  >  NC_000913/2027653‑2027832

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