Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,649,511 G→A E197E (GAG→GAA yhiS → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,649,5110GA99.6% 949.1 / NA 254E197E (GAG→GAAyhiShypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (137/116);  minor base T (0/1);  total (137/117)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGT  >  NC_000913/3649421‑3649600
                                                                                          |                                                                                         
tAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAAca                                                                                         <  1:3804728/93‑1 (MQ=255)
tAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAAca                                                                                         >  2:5423226/1‑93 (MQ=255)
tAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAAca                                                                                         <  1:2524062/93‑1 (MQ=255)
 aaTTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAacac                                                                                        >  1:875486/1‑93 (MQ=255)
 aaTTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAacac                                                                                        >  2:2800297/1‑93 (MQ=255)
    tttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGt                                                                                     >  1:5017058/1‑93 (MQ=255)
    tttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGt                                                                                     <  2:5316859/93‑1 (MQ=255)
     ttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTg                                                                                    <  1:4705112/93‑1 (MQ=255)
     ttATGGCTGAACAACTATCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTAACAAAAAATGACGAAGTTATATATTGAACACAGTg                                                                                    <  1:1584724/93‑1 (MQ=255)
      tATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGa                                                                                   <  1:4254173/93‑1 (MQ=255)
        ttGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaa                                                                                 <  1:2711346/93‑1 (MQ=255)
         tGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaaa                                                                                <  1:1855350/93‑1 (MQ=255)
         tGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaaa                                                                                <  2:1173196/93‑1 (MQ=255)
           cTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATa                                                                              <  1:2255576/93‑1 (MQ=255)
            tGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATaa                                                                             >  1:2165170/1‑93 (MQ=255)
            tGAACAACGAGCTGACTTACCAACTAATATAAAAAATATTCAAAACCTGTATCAAAAAATGACGGAGTTATATATTGAAAACAGTGAAAATaa                                                                             >  2:2660527/1‑93 (MQ=255)
               acaacGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATaaaaa                                                                          <  2:3053916/93‑1 (MQ=255)
               acaacGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATaaaaa                                                                          >  2:5941792/1‑93 (MQ=255)
                caacGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAAc                                                                         <  1:6243039/93‑1 (MQ=255)
                 aacGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACa                                                                        >  1:5092414/1‑93 (MQ=255)
                  acGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAg                                                                       >  2:1012669/1‑93 (MQ=255)
                  acGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAg                                                                       >  1:2067551/1‑93 (MQ=255)
                   cGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAgg                                                                      <  2:4868993/93‑1 (MQ=255)
                     aGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGAt                                                                    >  2:3454127/1‑93 (MQ=255)
                        tGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaa                                                                 <  2:5870326/93‑1 (MQ=255)
                         gACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaaa                                                                >  1:2016892/1‑93 (MQ=255)
                          aCTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAg                                                               >  1:3369322/1‑93 (MQ=255)
                           cTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGt                                                              <  1:3577144/93‑1 (MQ=255)
                           cTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGt                                                              >  2:2295903/1‑93 (MQ=255)
                            ttACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTc                                                             >  1:1623025/1‑93 (MQ=255)
                            ttACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTc                                                             >  2:2999229/1‑93 (MQ=255)
                             tACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCt                                                            >  1:95168/1‑93 (MQ=255)
                             tACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCt                                                            >  2:2727425/1‑93 (MQ=255)
                              aCCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           <  1:3291350/93‑1 (MQ=255)
                              aCCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           <  2:4604103/93‑1 (MQ=255)
                              aCCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           <  2:1138980/93‑1 (MQ=255)
                               ccAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCttt                                                          >  1:3644311/1‑93 (MQ=255)
                                cAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTg                                                         >  2:4698758/1‑93 (MQ=255)
                                 aaCTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGc                                                        >  1:3429607/1‑93 (MQ=255)
                                    tAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCgg                                                     >  2:2132212/1‑93 (MQ=255)
                                    tAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCgg                                                     <  2:3208982/93‑1 (MQ=255)
                                    tAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCgg                                                     >  1:3002865/1‑93 (MQ=255)
                                      atatAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGaa                                                   <  2:543850/93‑1 (MQ=255)
                                       tataAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAAc                                                  >  2:338234/1‑93 (MQ=255)
                                        ataAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGaatc                                                 >  1:2199563/1‑91 (MQ=255)
                                        ataAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAAcc                                                 <  1:3980003/93‑1 (MQ=255)
                                          aaaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCaa                                               >  2:3756488/1‑93 (MQ=255)
                                          aaaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCaa                                               >  2:1545417/1‑93 (MQ=255)
                                          aaaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCaa                                               >  1:2095941/1‑93 (MQ=255)
                                           aaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAAc                                              >  2:1354482/1‑93 (MQ=255)
                                            aaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACt                                             <  1:3960026/93‑1 (MQ=255)
                                             aGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACtt                                            <  2:5022017/93‑1 (MQ=255)
                                               atatTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTa                                          >  2:4946475/1‑93 (MQ=255)
                                               atatTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTa                                          <  1:526611/93‑1 (MQ=255)
                                                tatTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTAt                                         <  1:608180/93‑1 (MQ=255)
                                                 atTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTAtt                                        >  1:14276/1‑93 (MQ=255)
                                                  ttCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTg                                       <  2:3825907/93‑1 (MQ=255)
                                                  ttCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTg                                       >  1:5593106/1‑93 (MQ=255)
                                                   tCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGa                                      >  2:2442297/1‑93 (MQ=255)
                                                   tCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGa                                      <  2:628206/93‑1 (MQ=255)
                                                    cAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGAt                                     <  2:5776036/93‑1 (MQ=255)
                                                       aGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTc                                  >  2:513776/1‑93 (MQ=255)
                                                        gCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCa                                 >  2:339879/1‑93 (MQ=255)
                                                        gCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCa                                 <  1:1136204/93‑1 (MQ=255)
                                                         ccTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCaa                                <  1:3196015/93‑1 (MQ=255)
                                                         ccTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCaa                                >  2:54983/1‑93 (MQ=255)
                                                           tGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAata                              >  2:2333453/1‑93 (MQ=255)
                                                            gTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAatat                             <  2:2621421/93‑1 (MQ=255)
                                                            gTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAatat                             >  1:48314/1‑93 (MQ=255)
                                                             tATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATg                            >  1:634761/1‑93 (MQ=255)
                                                             tATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATg                            <  2:1653612/93‑1 (MQ=255)
                                                              aTCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATga                           >  1:4431920/1‑93 (MQ=255)
                                                              aTCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATga                           >  1:4757692/1‑93 (MQ=255)
                                                               tCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgac                          >  1:1000205/1‑93 (MQ=255)
                                                               tCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgac                          <  1:1214759/93‑1 (MQ=255)
                                                               tCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgac                          >  1:2478168/1‑93 (MQ=255)
                                                                 aaaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacag                        <  2:70348/93‑1 (MQ=255)
                                                                 aaaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacag                        <  2:4263291/93‑1 (MQ=255)
                                                                  aaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagg                       <  1:1248568/93‑1 (MQ=255)
                                                                  aaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTAGCCGGAACCAACTTTATTGATTTCAATATgacagg                       <  2:214382/93‑1 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      <  2:3491983/93‑1 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      <  2:4467243/93‑1 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      >  2:2559172/1‑93 (MQ=255)
                                                                    aaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacaggac                     <  2:433461/93‑1 (MQ=255)
                                                                     aaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacaggaca                    >  2:951331/1‑93 (MQ=255)
                                                                      aTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacaggacag                   <  1:5846905/93‑1 (MQ=255)
                                                                       tGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGa                  >  1:2874068/1‑93 (MQ=255)
                                                                       tGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGa                  <  1:5978761/93‑1 (MQ=255)
                                                                           gAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGTAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACATAAcca              >  2:136634/1‑92 (MQ=255)
                                                                           gAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCt              >  2:2744547/1‑93 (MQ=255)
                                                                            aaGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTg             >  1:1808308/1‑93 (MQ=255)
                                                                             aGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGt            >  1:5138084/1‑93 (MQ=255)
                                                                               ttATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTcc          >  2:1803125/1‑93 (MQ=255)
                                                                                tatataTTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCg         >  2:4591746/1‑93 (MQ=255)
                                                                                tatataTTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCg         >  1:2830038/1‑93 (MQ=255)
                                                                                tatataTTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCg         >  2:5639686/1‑93 (MQ=255)
                                                                                  tatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCggg       >  1:4295685/1‑93 (MQ=255)
                                                                                   atatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGgt      >  2:44027/1‑93 (MQ=255)
                                                                                      ttGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGgttcg   >  2:6103212/1‑93 (MQ=255)
                                                                                       tGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGgttcgt  <  1:610059/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGT  >  NC_000913/3649421‑3649600

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