Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,371,235 T→C intergenic (‑24/+98) nanA ← / ← b3226 N‑acetylneuraminate lyase/putative FADA‑type transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,371,2350TC95.0% 799.5 / 19.6 277intergenic (‑24/+98)nanA/b3226N‑acetylneuraminate lyase/putative FADA‑type transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (10/4);  new base C (131/132);  total (141/136)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

TCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGTCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCACTA  >  NC_000913/3371144‑3371324
                                                                                           |                                                                                         
 ccAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGTCt                                                                                         <  1:4363607/93‑1 (MQ=255)
   aGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGa                                                                                       <  1:3587015/93‑1 (MQ=255)
   aGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGa                                                                                       >  2:3613229/1‑93 (MQ=255)
   aGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGa                                                                                       <  1:1608776/93‑1 (MQ=255)
    gTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGAt                                                                                      <  1:3230954/93‑1 (MQ=255)
    gTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGAt                                                                                      >  1:1657943/1‑93 (MQ=255)
    gTGCTTGTTGTTCGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGAt                                                                                      <  1:1074369/93‑1 (MQ=255)
    gTGCTTGTTGTTAGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGAt                                                                                      <  1:83319/93‑1 (MQ=255)
     tGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATa                                                                                     <  2:3874996/93‑1 (MQ=255)
     tGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATa                                                                                     >  2:3345379/1‑93 (MQ=255)
      gCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATaa                                                                                    >  2:1689912/1‑93 (MQ=255)
       cTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATaaa                                                                                   <  1:2553534/93‑1 (MQ=255)
        ttgttgttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAAc                                                                                  <  1:478847/93‑1 (MQ=255)
        ttgttgttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAAc                                                                                  >  2:3801325/1‑93 (MQ=255)
        ttgttgttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAAc                                                                                  <  1:2060288/93‑1 (MQ=255)
           ttgttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGat                                                                               <  2:4322888/93‑1 (MQ=255)
           ttgttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGat                                                                               <  2:1713787/93‑1 (MQ=255)
             gttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGatat                                                                             <  1:23059/93‑1 (MQ=255)
             gttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGatat                                                                             >  2:65129/1‑93 (MQ=255)
             gttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGatat                                                                             <  1:579333/93‑1 (MQ=255)
              ttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGAtata                                                                            >  1:4035535/1‑93 (MQ=255)
              ttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGAtata                                                                            <  2:3977184/93‑1 (MQ=255)
              ttgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGAtata                                                                            <  2:2522125/93‑1 (MQ=255)
               tgGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATAc                                                                           <  1:2031873/93‑1 (MQ=255)
                 gTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCt                                                                         <  1:131592/93‑1 (MQ=255)
                  tCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCtt                                                                        <  2:2352739/93‑1 (MQ=255)
                   cAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCttt                                                                       <  2:1994767/93‑1 (MQ=255)
                    aaaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTta                                                                      <  2:3041010/93‑1 (MQ=255)
                    aaaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTta                                                                      <  1:884039/93‑1 (MQ=255)
                    aaaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTta                                                                      <  1:1766034/93‑1 (MQ=255)
                     aaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTtat                                                                     >  1:3671741/1‑93 (MQ=255)
                     aaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTtat                                                                     >  1:1375786/1‑93 (MQ=255)
                     aaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTtat                                                                     >  1:2521575/1‑93 (MQ=255)
                     aaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTtat                                                                     >  2:2818584/1‑93 (MQ=255)
                     aaaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTtat                                                                     >  2:2676438/1‑93 (MQ=255)
                      aaGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTCGCCATAAATACCTCTGAAGTGATGCTTGTCTGATAAACGATATACCTTtata                                                                    >  1:182878/1‑93 (MQ=255)
                        ggAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATAcc                                                                  <  2:3282845/93‑1 (MQ=255)
                        ggAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATAcc                                                                  <  1:2707388/93‑1 (MQ=255)
                         gAGTCAGGAGTGCGGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCt                                                                 >  1:2214224/1‑93 (MQ=255)
                           gTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGt                                                               >  1:3840913/1‑93 (MQ=255)
                             cAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTta                                                             <  2:3187516/93‑1 (MQ=255)
                               ggAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTtata                                                           >  1:757936/1‑93 (MQ=255)
                               ggAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTtata                                                           <  2:4344539/93‑1 (MQ=255)
                                gAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATAc                                                          <  2:2733038/93‑1 (MQ=255)
                                gAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATAc                                                          >  2:1833697/1‑93 (MQ=255)
                                gAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATAc                                                          <  1:3880857/93‑1 (MQ=255)
                                    gCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGa                                                      <  2:1105939/93‑1 (MQ=255)
                                    gCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGa                                                      <  2:4107684/93‑1 (MQ=255)
                                       gCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCa                                                   <  1:1241434/93‑1 (MQ=255)
                                        ccATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCaa                                                  <  1:3214615/93‑1 (MQ=255)
                                         cATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAAt                                                 >  1:1172084/1‑93 (MQ=255)
                                           ttACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTa                                               <  1:4302157/93‑1 (MQ=255)
                                           ttACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTa                                               <  2:3925938/93‑1 (MQ=255)
                                             aCGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAg                                             >  1:4462979/1‑93 (MQ=255)
                                               gCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGTCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCa                                           >  2:1644057/1‑93 (MQ=255)
                                               gCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCa                                           >  1:1720656/1‑93 (MQ=255)
                                                ccACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCaa                                          >  2:2061113/1‑93 (MQ=255)
                                                 cACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAac                                         >  2:1475191/1‑93 (MQ=255)
                                                 cACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAac                                         >  1:4528747/1‑93 (MQ=255)
                                                  aCGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGTCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAaca                                        >  2:370592/1‑93 (MQ=255)
                                                  aCGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAaca                                        >  1:2863587/1‑93 (MQ=255)
                                                   cGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAacac                                       >  1:3732102/1‑93 (MQ=255)
                                                   cGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAacac                                       >  1:2355990/1‑93 (MQ=255)
                                                    gTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACAcc                                      >  2:1276664/1‑93 (MQ=255)
                                                    gTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACAcc                                      >  1:4055968/1‑93 (MQ=255)
                                                      aaaTTCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACAcccc                                    >  2:2993763/1‑93 (MQ=255)
                                                         ttCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATa                                 >  1:3522958/1‑93 (MQ=255)
                                                          tCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATAc                                <  1:2953287/93‑1 (MQ=255)
                                                          tCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATAc                                <  2:4625404/93‑1 (MQ=255)
                                                          tCGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATAc                                <  2:1210793/93‑1 (MQ=255)
                                                           cGTTGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACa                               <  1:3430830/93‑1 (MQ=255)
                                                              tGCCATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGaa                            >  2:862613/1‑93 (MQ=255)
                                                                 cATAAATACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGc                         <  1:997136/93‑1 (MQ=255)
                                                                    aaaTACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTta                      <  2:479794/93‑1 (MQ=255)
                                                                     aaTAGCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTtat                     >  2:4539769/1‑93 (MQ=255)
                                                                      aTACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTtata                    <  1:1862469/93‑1 (MQ=255)
                                                                      aTACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTtata                    <  2:3517195/93‑1 (MQ=255)
                                                                      aTACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTtata                    <  2:3555780/93‑1 (MQ=255)
                                                                       tACCTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATaa                   >  1:2174083/1‑93 (MQ=255)
                                                                         ccTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATg                 >  2:4661477/1‑93 (MQ=255)
                                                                         ccTCTGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATg                 >  2:2911917/1‑93 (MQ=255)
                                                                             tGAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGAt             >  2:3151104/1‑93 (MQ=255)
                                                                              gAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATc            >  1:96040/1‑93 (MQ=255)
                                                                              gAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATc            <  1:3311837/93‑1 (MQ=255)
                                                                              gAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATc            <  1:413827/93‑1 (MQ=255)
                                                                              gAAGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATc            <  2:4020842/93‑1 (MQ=255)
                                                                               aaGTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTAAACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCt           >  2:3427583/1‑93 (MQ=255)
                                                                                 gTGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGc         >  1:3544383/1‑93 (MQ=255)
                                                                                  tGATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCt        >  1:1531335/1‑93 (MQ=255)
                                                                                   gATGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCtt       <  1:2760235/93‑1 (MQ=255)
                                                                                    aTGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTc      <  1:3102318/93‑1 (MQ=255)
                                                                                     tGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCa     >  1:418383/1‑93 (MQ=255)
                                                                                     tGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCa     >  2:4468098/1‑93 (MQ=255)
                                                                                     tGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCa     >  2:1608619/1‑93 (MQ=255)
                                                                                     tGCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCa     >  2:357751/1‑93 (MQ=255)
                                                                                      gCTTGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCAc    <  2:2416407/93‑1 (MQ=255)
                                                                                        ttGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCACTa  <  2:4581798/93‑1 (MQ=255)
                                                                                        ttGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCACTa  >  2:1779344/1‑93 (MQ=255)
                                                                                        ttGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCACTa  >  2:590822/1‑93 (MQ=255)
                                                                                        ttGCCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCACTa  <  2:2611977/93‑1 (MQ=255)
                                                                                           |                                                                                         
TCCAGTGCTTGTTGTTGGTCAAAAGGAGTCAGGAGTGCAGCCATTACGCCACGTAAATTCGTTGCCATAAATACCTCTGAAGTGATGCTTGTCTGATAAACGATATACCTTTATACCTGTTATACCAGATCAATTAAGCAACACCCCATACAGAAAGCTTATAATGCGATCTGCTTCACTA  >  NC_000913/3371144‑3371324

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