Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,383,404 A→G T95A (ACC→GCC)  yhcN → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,383,4040AG99.7% 1090.8 / ‑5.5 348T95A (ACC→GCC) yhcNhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (168/178);  total (168/180)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 357 total aligned reads displayed.

TCGCCAATGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCT  >  NC_000913/3383314‑3383493
                                                                                          |                                                                                         
tCGCCAATGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGcc                                                                                         <  2:5440267/93‑1 (MQ=255)
tCGCCAATGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGcc                                                                                         >  2:757748/1‑93 (MQ=255)
  gCCAATGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTg                                                                                       <  1:958890/93‑1 (MQ=255)
   ccAATGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTgg                                                                                      <  1:2869869/93‑1 (MQ=255)
       tGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACg                                                                                  >  1:836424/1‑93 (MQ=255)
        ggATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGc                                                                                 <  1:1090903/93‑1 (MQ=255)
         gATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCt                                                                                >  1:5660236/1‑93 (MQ=255)
          atatGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTa                                                                               <  1:5194275/93‑1 (MQ=255)
           tatGCGTGAAATGCTGAACGAAAAGGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTaa                                                                              <  1:7023597/93‑2 (MQ=255)
           tatGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTAc                                                                              <  2:5883434/93‑1 (MQ=255)
           tatGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTAc                                                                              >  2:627165/1‑93 (MQ=255)
           tatGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTAc                                                                              <  1:6653477/93‑1 (MQ=255)
            atGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACg                                                                             <  2:7686937/93‑1 (MQ=255)
             tGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACgg                                                                            <  1:1960390/93‑1 (MQ=255)
             tGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACgg                                                                            >  1:1882491/1‑93 (MQ=255)
              gCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGc                                                                           >  2:2806128/1‑93 (MQ=255)
               cGTGAAATGCTGAACAAAAAAGCGGACGAGAACGGCGCGATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGACACGCTACGGCt                                                                          <  2:7037990/93‑1 (MQ=255)
               cGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCt                                                                          <  2:5158137/93‑1 (MQ=255)
                gTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTg                                                                         >  1:2955676/1‑93 (MQ=255)
                gTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTg                                                                         >  2:6929392/1‑93 (MQ=255)
                 tGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGa                                                                        >  1:1642779/1‑93 (MQ=255)
                  gAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGaa                                                                       <  1:1590917/93‑1 (MQ=255)
                  gAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGaa                                                                       <  2:1143047/93‑1 (MQ=255)
                     aTGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACTCTAAGGCTGAACTg                                                                    >  1:6209065/1‑93 (MQ=255)
                      tGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGt                                                                   >  1:4007920/1‑93 (MQ=255)
                          gAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACaa                                                               <  2:4665210/93‑1 (MQ=255)
                          gAACAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACaa                                                               >  2:1756918/1‑93 (MQ=255)
                           aaCAAAAAAGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACaaa                                                              <  1:3771497/93‑1 (MQ=255)
                              aaaaaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATaa                                                           >  2:2118954/1‑93 (MQ=255)
                              aaaaaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATaa                                                           >  2:1692169/1‑93 (MQ=255)
                               aaaaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATaaa                                                          >  1:351937/1‑93 (MQ=255)
                               aaaaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATaaa                                                          <  1:3476621/93‑1 (MQ=255)
                               aaaaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGGACAAATaaa                                                          >  2:6711832/1‑93 (MQ=255)
                                aaaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAAc                                                         <  1:1785361/93‑1 (MQ=255)
                                 aaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAAcc                                                        <  2:264357/93‑1 (MQ=255)
                                 aaaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAAcc                                                        <  2:7364468/93‑1 (MQ=255)
                                  aaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAAccc                                                       <  1:818565/93‑1 (MQ=255)
                                  aaGCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAAccc                                                       >  1:3336682/1‑93 (MQ=255)
                                    gCGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTc                                                     <  1:6027785/93‑1 (MQ=255)
                                     cGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCa                                                    >  1:3859117/1‑93 (MQ=255)
                                     cGGAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCa                                                    <  2:2546217/93‑1 (MQ=255)
                                      ggAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCAt                                                   <  2:2893839/93‑1 (MQ=255)
                                       gAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATc                                                  >  1:7406162/1‑93 (MQ=255)
                                       gAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATc                                                  <  1:2540327/93‑1 (MQ=255)
                                       gAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATc                                                  <  2:5025798/93‑1 (MQ=255)
                                       gAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATc                                                  <  1:1335509/93‑1 (MQ=255)
                                       gAAGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATc                                                  >  1:265533/1‑93 (MQ=255)
                                        aaGAGAAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCg                                                 >  1:383088/1‑93 (MQ=255)
                                           agaAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCt                                              >  2:1467483/1‑93 (MQ=255)
                                            gaAAGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCtt                                             >  2:4856281/1‑93 (MQ=255)
                                             aaaGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCCTCGTCTTg                                            >  1:4409422/1‑93 (MQ=255)
                                               aGGCGCAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTc                                          >  1:4156093/1‑93 (MQ=255)
                                                 gcgcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcg                                        <  1:7356277/93‑1 (MQ=255)
                                                 gcgcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcg                                        <  1:3683990/93‑1 (MQ=255)
                                                  cgcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcga                                       <  2:7419229/93‑1 (MQ=255)
                                                  cgcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcga                                       <  1:2031749/93‑1 (MQ=255)
                                                  cgcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcga                                       <  2:4228761/93‑1 (MQ=255)
                                                  cgcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcga                                       <  1:233766/93‑1 (MQ=255)
                                                   gcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcgac                                      >  2:1182421/1‑93 (MQ=255)
                                                   gcAATGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcgac                                      <  2:370668/93‑1 (MQ=255)
                                                     aatggCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCcgacga                                    <  2:5360336/90‑1 (MQ=255)
                                                      atggCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGat                                   >  1:2586652/3‑93 (MQ=255)
                                                      atggCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGat                                   >  2:6092924/3‑93 (MQ=255)
                                                        ggCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGatat                                 <  2:3572914/93‑1 (MQ=255)
                                                         gCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATAtt                                <  2:6656086/93‑1 (MQ=255)
                                                         gCCTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATAtt                                >  1:1728339/1‑93 (MQ=255)
                                                           cTACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGc                              <  1:3174372/93‑1 (MQ=255)
                                                            tACCAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGcc                             >  1:7074861/1‑93 (MQ=255)
                                                              ccAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGcccc                           <  2:6404205/93‑1 (MQ=255)
                                                               cAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGccccc                          >  2:3221221/1‑93 (MQ=255)
                                                               cAGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGccccc                          <  2:3831609/93‑1 (MQ=255)
                                                                aGATTACTGAAGCTCGTAGCGGTGACTCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCg                         <  1:5342218/93‑1 (MQ=255)
                                                                aGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCg                         <  2:3495128/93‑1 (MQ=255)
                                                                aGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCg                         >  2:3972867/1‑93 (MQ=255)
                                                                aGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCg                         >  1:2266056/1‑93 (MQ=255)
                                                                aGATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCg                         <  1:7587308/93‑1 (MQ=255)
                                                                 gATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCgg                        >  1:1473757/1‑93 (MQ=255)
                                                                 gATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCgg                        >  1:5063185/1‑93 (MQ=255)
                                                                 gATTACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCgg                        <  2:5940705/93‑1 (MQ=255)
                                                                   ttACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGTCCCCGGtt                      >  1:447459/1‑93 (MQ=255)
                                                                   ttACTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGtt                      >  1:7502558/1‑93 (MQ=255)
                                                                     tcTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCg                    <  2:7747003/92‑1 (MQ=255)
                                                                     aCTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGTCCCCGGTTCg                    >  2:950134/1‑93 (MQ=255)
                                                                      cTGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCgg                   <  1:68829/93‑1 (MQ=255)
                                                                       tGAAGCTCGTAGCGTTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                  <  2:5992666/93‑1 (MQ=255)
                                                                       tGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                  <  1:5585561/93‑1 (MQ=255)
                                                                       tGAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                  >  2:2456092/1‑93 (MQ=255)
                                                                       ggAAGCGCGTAGCGGGGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggg                  <  1:6409218/92‑1 (MQ=255)
                                                                        gAAGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCgggg                 <  1:4850039/93‑1 (MQ=255)
                                                                         aaGCTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCggggg                <  2:7860711/93‑1 (MQ=255)
                                                                            cTCGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtt             <  2:3161379/93‑1 (MQ=255)
                                                                              cGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttt           <  1:790513/93‑1 (MQ=255)
                                                                              cGTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttt           <  2:1983671/93‑1 (MQ=255)
                                                                               gTAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCttttt          >  1:1014969/1‑93 (MQ=255)
                                                                                tAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttttt         >  1:3113019/1‑93 (MQ=255)
                                                                                tAGCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCcttttt         >  2:5776985/1‑93 (MQ=255)
                                                                                  gCGGTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCtttttttt       <  1:4622764/93‑1 (MQ=255)
                                                                                     gTGACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTgcg    >  1:4333699/1‑93 (MQ=255)
                                                                                       gACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCt  >  1:7701725/1‑93 (MQ=255)
                                                                                       gACGCCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCt  <  2:7205281/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCGCCAATGGATATGCGTGAAATGCTGAACAAAAAAGCGGAAGAGAAAGGCGCAACGGCCTACCAGATTACTGAAGCTCGTAGCGGTGACACCTGGCACGCTACGGCTGAACTGTACAAATAAACCCTCATCGTCTTGTCCGACGATATTGCCCCCGGTTCGGGGGCTTTTTTTTGCGCT  >  NC_000913/3383314‑3383493

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: