Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 651,666 T→C V70A (GTC→GCC)  citA → putative sensor‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913651,6660TC99.5% 702.8 / ‑4.9 213V70A (GTC→GCC) citAputative sensor‑type protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (98/114);  total (98/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

Only 100 of 220 total aligned reads displayed.

TTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGA  >  NC_000913/651577‑651755
                                                                                         |                                                                                         
ttttACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCAt                                                                                        >  1:4592550/1‑93 (MQ=255)
 tttACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCAtc                                                                                       >  2:1575286/1‑93 (MQ=255)
  ttACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGTCAtct                                                                                      <  1:3912737/93‑1 (MQ=255)
   taaggCCAGTTTGTAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCAtctc                                                                                     <  1:185997/90‑1 (MQ=255)
   tACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCAtctc                                                                                     >  2:3822284/1‑93 (MQ=255)
   tACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCAtctc                                                                                     <  2:2757765/93‑1 (MQ=255)
    aCGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATctct                                                                                    >  2:3805733/1‑93 (MQ=255)
    aCGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATctct                                                                                    <  2:2866011/93‑1 (MQ=255)
      ggCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGc                                                                                  <  2:4184926/93‑1 (MQ=255)
      ggCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGc                                                                                  <  1:3555599/93‑1 (MQ=255)
       gCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCg                                                                                 >  1:1287885/1‑93 (MQ=255)
         cAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGt                                                                               <  1:3524266/93‑1 (MQ=255)
         cAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGt                                                                               <  1:2602877/93‑1 (MQ=255)
         cAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGt                                                                               >  2:2828175/1‑93 (MQ=255)
          aGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTg                                                                              >  1:4018819/1‑93 (MQ=255)
               tGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAAc                                                                         >  1:1024244/1‑93 (MQ=255)
               tGAGGACTATGTAACGCTTCATGTACGCGACATGTCAATGAATCAGTCGAAAATTATTGCCTCCAATGACAGTGCCACCTCTGCGGTGAAAAc                                                                         <  2:2255331/93‑1 (MQ=255)
                  ggACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAAcgcg                                                                      >  1:3588817/1‑93 (MQ=255)
                   gACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGt                                                                     <  2:1516078/93‑1 (MQ=255)
                   gACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGt                                                                     <  1:1777551/93‑1 (MQ=255)
                   gACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGt                                                                     <  2:3005777/93‑1 (MQ=255)
                   gACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGt                                                                     <  2:3649162/93‑1 (MQ=255)
                      tATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGAc                                                                  <  1:2978380/93‑1 (MQ=255)
                       aTTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACt                                                                 >  2:2440641/1‑93 (MQ=255)
                            aCGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACaaa                                                            <  1:2694071/93‑1 (MQ=255)
                            aCGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACaaa                                                            <  1:3192193/93‑1 (MQ=255)
                             cGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAAc                                                           <  1:651799/93‑1 (MQ=255)
                              gCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACg                                                          >  1:3310510/1‑93 (MQ=255)
                              gCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACg                                                          >  2:4412173/1‑93 (MQ=255)
                              gCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACg                                                          >  1:1056854/1‑93 (MQ=255)
                                ttCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACACACGGc                                                        >  1:1206019/1‑93 (MQ=255)
                                ttCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGc                                                        >  2:815891/1‑93 (MQ=255)
                                ttCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGc                                                        <  1:1268315/93‑1 (MQ=255)
                                 tCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGTCt                                                       >  1:4714079/1‑93 (MQ=255)
                                 tCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCt                                                       <  1:2338300/93‑1 (MQ=255)
                                  cATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGAGGTGAAAACGCGTGACTACAAACGGCTg                                                      <  2:1391522/93‑1 (MQ=255)
                                   aTGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTgg                                                     >  1:1168175/1‑93 (MQ=255)
                                    tGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGc                                                    >  1:4434252/1‑93 (MQ=255)
                                    tGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGc                                                    <  2:3318751/93‑1 (MQ=255)
                                    tGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGc                                                    <  2:3728661/93‑1 (MQ=255)
                                       aCGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGAc                                                 >  1:2451007/1‑93 (MQ=255)
                                       aCGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGAc                                                 <  1:2221762/93‑1 (MQ=255)
                                        cgcgACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGAcc                                                <  2:4061469/93‑1 (MQ=255)
                                         gcgACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCa                                               <  2:2563057/93‑1 (MQ=255)
                                          cgACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCAt                                              >  1:2251139/1‑93 (MQ=255)
                                          cgACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCAt                                              <  2:4454688/93‑1 (MQ=255)
                                          cgACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCAt                                              <  1:505714/93‑1 (MQ=255)
                                          cgACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCAt                                              <  2:2124/93‑1 (MQ=255)
                                             cATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGc                                           <  1:4242653/93‑1 (MQ=255)
                                               tGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTa                                         >  1:4559434/1‑93 (MQ=255)
                                                 gCAATGAAGCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAACCGGCTGGCGACCATCGCTAAc                                       <  1:851625/93‑1 (MQ=255)
                                                  cAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACa                                      <  1:2279091/93‑1 (MQ=255)
                                                  cAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACa                                      >  1:4510667/1‑93 (MQ=255)
                                                   aaTGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACaa                                     <  1:4146614/93‑1 (MQ=255)
                                                   aaTGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACaa                                     <  2:1460148/93‑1 (MQ=255)
                                                    aTGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACaaa                                    >  1:1642569/1‑93 (MQ=255)
                                                    aTGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACaaa                                    <  2:8830/93‑1 (MQ=255)
                                                     tGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAt                                   >  2:1937499/1‑93 (MQ=255)
                                                     tGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAt                                   >  2:1662073/1‑93 (MQ=255)
                                                      gAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAtt                                  <  2:3603553/93‑1 (MQ=255)
                                                      gAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAtt                                  <  1:262951/93‑1 (MQ=255)
                                                      gAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAtt                                  <  1:782903/93‑1 (MQ=255)
                                                      gAATCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAtt                                  >  1:3429055/1‑93 (MQ=255)
                                                        aTCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTAc                                >  2:4741417/1‑93 (MQ=255)
                                                         tCAGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGAGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACa                               >  1:3040009/1‑93 (MQ=255)
                                                           aGGCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACaaa                             >  1:3680260/1‑93 (MQ=255)
                                                            ggggAAAATTATTGCCTCCAATGTCAGTGCGATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAag                            <  2:3378437/90‑1 (MQ=255)
                                                            ggCGAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAag                            <  2:244569/93‑1 (MQ=255)
                                                               gAAAATTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAgaga                         <  2:2893533/93‑1 (MQ=255)
                                                                aaaaTTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGAt                        <  1:2119914/93‑1 (MQ=255)
                                                                  aaTTTTTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATAc                      <  1:1515615/93‑1 (MQ=255)
                                                                  aaTTATTGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATAc                      <  1:1515617/93‑1 (MQ=255)
                                                                   attattGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATAcc                     >  2:1980928/1‑93 (MQ=255)
                                                                   attattGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATAcc                     <  1:1836693/93‑1 (MQ=255)
                                                                      attGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAt                  >  1:877690/1‑93 (MQ=255)
                                                                      attGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAt                  <  1:1714285/93‑1 (MQ=255)
                                                                       ttGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAtt                 <  2:2157064/93‑1 (MQ=255)
                                                                        tGCCTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAttt                <  1:3818554/93‑1 (MQ=255)
                                                                          ccTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTg              <  2:1618367/93‑1 (MQ=255)
                                                                           cTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGa             >  2:704363/1‑93 (MQ=255)
                                                                           cTCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGAGACCGATTTTgt             >  2:123674/1‑92 (MQ=255)
                                                                            tCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACGAATTACAAAGAGATACCGATTTTGAt            <  2:3993182/93‑1 (MQ=255)
                                                                            tCCAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGAt            <  1:2481846/93‑1 (MQ=255)
                                                                             ccAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGAtt           >  2:2976854/1‑93 (MQ=255)
                                                                              cAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTa          >  1:4584111/1‑93 (MQ=255)
                                                                              cAATGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTa          <  1:3644918/93‑1 (MQ=255)
                                                                               aaTGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGTGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAt         >  1:694193/1‑93 (MQ=255)
                                                                               aaTGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAt         <  1:352753/93‑1 (MQ=255)
                                                                               aaTGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAt         >  1:1003539/1‑93 (MQ=255)
                                                                                aTGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAtg        <  1:4301325/93‑1 (MQ=255)
                                                                                 tGACAGTGCCATCTCTGCGGTGAGAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGCTTAtgt       >  2:626778/1‑93 (MQ=255)
                                                                                 tGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAtgt       <  2:13618/93‑1 (MQ=255)
                                                                                 tGACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAtgt       <  2:180065/93‑1 (MQ=255)
                                                                                  gACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAtgtg      >  2:4711778/1‑93 (MQ=255)
                                                                                  gACAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTAtgtg      >  2:566950/1‑93 (MQ=255)
                                                                                    cAGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATgtggt    <  1:248895/93‑1 (MQ=255)
                                                                                     aGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATgtggtg   >  2:3494961/1‑93 (MQ=255)
                                                                                     aGTGCCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATgtggtg   >  2:1304224/1‑93 (MQ=255)
                                                                                      gtgcCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGa  <  1:590763/93‑1 (MQ=255)
                                                                                      gtgcCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGa  >  2:1911183/1‑93 (MQ=255)
                                                                                         |                                                                                         
TTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGA  >  NC_000913/651577‑651755

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