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.9% 2881.1 / NA 780E197E (GAG→GAAyhiShypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (361/418);  total (361/419)
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.00e+00

Only 100 of 796 total aligned reads displayed.

TAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGT  >  NC_000913/3649421‑3649600
                                                                                          |                                                                                         
tAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAAca                                                                                         <  2:3984801/93‑1 (MQ=255)
tAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAAca                                                                                         <  2:3015786/93‑1 (MQ=255)
 aaTTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAacac                                                                                        >  2:5353610/1‑93 (MQ=255)
  aTTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAAcaca                                                                                       <  2:7187491/93‑1 (MQ=255)
   ttttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAg                                                                                      <  1:8433065/93‑1 (MQ=255)
   ttttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAg                                                                                      <  1:8613599/93‑1 (MQ=255)
    tttATGGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGt                                                                                     <  2:6225841/93‑1 (MQ=255)
     ttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTg                                                                                    <  1:4650597/93‑1 (MQ=255)
     ttATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTg                                                                                    <  1:3492555/93‑1 (MQ=255)
      tATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGa                                                                                   >  1:5191612/1‑93 (MQ=255)
       aTTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaa                                                                                  >  2:2513330/1‑93 (MQ=255)
       aTTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaa                                                                                  <  2:897138/93‑1 (MQ=255)
       aTTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaa                                                                                  >  1:3947892/1‑93 (MQ=255)
       aTTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaa                                                                                  >  2:5060589/1‑93 (MQ=255)
        ttGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaa                                                                                 >  1:182971/1‑93 (MQ=255)
        ttGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaa                                                                                 <  1:2109718/93‑1 (MQ=255)
        ttGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaa                                                                                 <  2:7377030/93‑1 (MQ=255)
        ttGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGaaa                                                                                 <  1:5766453/93‑1 (MQ=255)
           cTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATa                                                                              <  2:604783/93‑1 (MQ=255)
             gAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATaaa                                                                            >  1:10009428/1‑93 (MQ=255)
                caacGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAAc                                                                         <  2:1406756/93‑1 (MQ=255)
                 aacGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACa                                                                        <  1:6570737/93‑1 (MQ=255)
                  acGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAg                                                                       >  2:1738681/1‑93 (MQ=255)
                  acGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAg                                                                       <  2:955325/93‑1 (MQ=255)
                    gAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGa                                                                     >  1:4309315/1‑93 (MQ=255)
                    gAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGa                                                                     >  1:9465211/1‑93 (MQ=255)
                     aGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGAt                                                                    >  1:7204474/1‑93 (MQ=255)
                       cTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGa                                                                  >  2:8807061/1‑93 (MQ=255)
                        tGACTTACCAACTAATATTAAAGATATTCAAAGCCTGCATCAACTAATGACGAAGTTGTATTTTTAAGATAGTGAATATATAAACATGATGaa                                                                 >  2:2015333/1‑93 (MQ=255)
                        tGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaa                                                                 >  1:7578243/1‑93 (MQ=255)
                        tGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaa                                                                 >  2:3374951/1‑93 (MQ=255)
                         gACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaaa                                                                <  1:6733382/93‑1 (MQ=255)
                         gACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaaa                                                                <  2:8226876/93‑1 (MQ=255)
                         gACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGaaa                                                                <  2:3748587/93‑1 (MQ=255)
                          aCTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAg                                                               >  1:7869083/1‑93 (MQ=255)
                           cTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGt                                                              <  2:4943621/93‑1 (MQ=255)
                           cTTAACAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGt                                                              <  1:5343527/93‑1 (MQ=255)
                            ttACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTc                                                             <  2:2271124/93‑1 (MQ=255)
                            ttACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTc                                                             >  2:4645079/1‑93 (MQ=255)
                             tACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCt                                                            >  2:4419527/1‑93 (MQ=255)
                             tACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCt                                                            <  1:8074608/93‑1 (MQ=255)
                              acaaaCTAATATAAAATATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           <  2:7494543/90‑1 (MQ=255)
                              aCCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           <  1:7741886/93‑1 (MQ=255)
                              aCCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           >  2:5836543/1‑93 (MQ=255)
                              aCCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCtt                                                           <  1:897023/93‑1 (MQ=255)
                               ccAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCttt                                                          >  2:4305802/1‑93 (MQ=255)
                                 aaCTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGc                                                        <  2:1160171/93‑1 (MQ=255)
                                  aCTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGcc                                                       <  1:1675243/93‑1 (MQ=255)
                                    tAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCgg                                                     <  2:3602655/93‑1 (MQ=255)
                                    tAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCgg                                                     >  1:1989032/1‑93 (MQ=255)
                                       tataAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAAc                                                  <  1:6977385/93‑1 (MQ=255)
                                       tacaaaaGATACTCAAAGCATGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAAc                                                  <  2:6918911/90‑1 (MQ=255)
                                        ataAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAAcc                                                 >  1:2898270/1‑93 (MQ=255)
                                         taAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGNTGAAAGTCTTTGCCGGAACCa                                                >  1:1113391/1‑93 (MQ=255)
                                          aaaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCaa                                               <  1:1442463/93‑1 (MQ=255)
                                           aaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATCTATTGAACACAGTGAAAATAAACACAGGATGAAAGTCTTTGCCGGAACCAAc                                              >  1:5513812/1‑93 (MQ=255)
                                           aaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAAc                                              <  1:8156998/93‑1 (MQ=255)
                                           aaaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAAc                                              <  2:8069306/93‑1 (MQ=255)
                                            aaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACt                                             <  2:5196732/93‑1 (MQ=255)
                                            aaGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACt                                             <  2:9384827/93‑1 (MQ=255)
                                             aGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACtt                                            >  1:9636032/1‑93 (MQ=255)
                                                tatTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGCACCAACTTTAt                                         >  1:6420862/1‑93 (MQ=255)
                                                  ttCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTg                                       >  2:6689793/1‑93 (MQ=255)
                                                  ttCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTg                                       <  2:2662093/93‑1 (MQ=255)
                                                   tCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGa                                      <  1:2261660/93‑1 (MQ=255)
                                                     aaaGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGAtt                                    >  1:3819200/1‑93 (MQ=255)
                                                     aaaGCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGAtt                                    >  2:7950921/1‑93 (MQ=255)
                                                        gCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCa                                 >  1:4153081/1‑93 (MQ=255)
                                                        gCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCa                                 >  2:4837672/1‑93 (MQ=255)
                                                        gCCTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCa                                 >  2:972601/1‑93 (MQ=255)
                                                          cTGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAat                               <  1:1317967/93‑1 (MQ=255)
                                                           tGTATCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAata                              <  1:8684772/93‑1 (MQ=255)
                                                             tATCAAAAAATGACGAAGTTACATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATg                            >  2:6599651/1‑93 (MQ=255)
                                                              aTCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATga                           <  1:9267978/93‑1 (MQ=255)
                                                              aTCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATga                           <  2:10000845/93‑1 (MQ=255)
                                                               tCAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgac                          >  2:7087385/1‑93 (MQ=255)
                                                                cAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgaca                         <  1:3342316/93‑1 (MQ=255)
                                                                cAAAAAATGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgaca                         <  1:4898575/93‑1 (MQ=255)
                                                                 aaaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacag                        >  1:3625674/1‑93 (MQ=255)
                                                                 aaaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacag                        >  2:5914313/1‑93 (MQ=255)
                                                                  aaaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagg                       >  1:5023369/1‑93 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      <  1:4544157/93‑1 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      >  1:6050685/1‑93 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      <  1:9864381/93‑1 (MQ=255)
                                                                   aaaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacagga                      <  2:4540535/93‑1 (MQ=255)
                                                                    aaaTGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATgacaggac                     <  1:8812054/93‑1 (MQ=255)
                                                                       tGACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGa                  >  1:9168504/1‑93 (MQ=255)
                                                                        gACGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAAGACAGaa                 >  2:3877550/1‑93 (MQ=255)
                                                                          cGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAAcc               <  2:5670538/93‑1 (MQ=255)
                                                                          cGAAGTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAAcc               <  1:2746676/93‑1 (MQ=255)
                                                                              gTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTc           <  2:7853125/93‑1 (MQ=255)
                                                                              gTTATATATTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTc           >  2:8472941/1‑93 (MQ=255)
                                                                                tatataTTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCg         >  1:2525855/1‑93 (MQ=255)
                                                                                 atatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCgg        >  2:1259692/1‑93 (MQ=255)
                                                                                 atatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCgg        <  2:9678287/93‑1 (MQ=255)
                                                                                 atatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCgg        >  2:1538721/1‑93 (MQ=255)
                                                                                  tatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCggg       >  2:9187639/1‑93 (MQ=255)
                                                                                  tatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCggg       >  1:7365642/1‑93 (MQ=255)
                                                                                  tatatTGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCggg       >  2:4781459/1‑93 (MQ=255)
                                                                                       tGAACACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGgttcgt  <  2:6421203/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAATTTTATTGCTGAACAACGAGCTGACTTACCAACTAATATAAAAGATATTCAAAGCCTGTATCAAAAAATGACGAAGTTATATATTGAGCACAGTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGT  >  NC_000913/3649421‑3649600

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