Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,151 T→C 15.7% intergenic (+18/‑120) ygdQ → / → ygdR UPF0053 family inner membrane protein/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,971,1510TC15.7% 512.0 / 81.3 217intergenic (+18/‑120)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (80/103);  new base C (13/21);  total (93/124)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 221 total aligned reads displayed.

AGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTA  >  NC_000913/2971052‑2971250
                                                                                                   |                                                                                                   
aGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTCGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCt                                                                                                     >  2:728508/1‑100 (MQ=255)
   ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCg                                                                                                  <  2:1267479/100‑1 (MQ=255)
   ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGca                                                                                                  <  2:1864922/100‑2 (MQ=255)
    tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCgg                                                                                                 <  2:3814003/100‑1 (MQ=255)
     aCATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                <  2:2454988/100‑1 (MQ=255)
      cATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGa                                                                                               <  1:1247117/100‑1 (MQ=255)
      cATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGa                                                                                               >  2:1573802/1‑100 (MQ=255)
        tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                             <  2:1941849/100‑1 (MQ=255)
        tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                             >  1:2433709/1‑100 (MQ=255)
        tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                             >  2:2362651/1‑100 (MQ=255)
         cTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                            <  1:3136387/100‑1 (MQ=255)
          tACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGc                                                                                           <  2:3188026/100‑1 (MQ=255)
           aCTTCGCTATGATCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGGTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                          >  1:2314625/1‑100 (MQ=255)
           aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                          >  1:1576990/1‑100 (MQ=255)
            cTTCGCGATGTTCTTCTCTATTGCGGTTGAAATCCTCAACTTGAGTCGCAACATAAAGAATCCGCTCTAATACTCCGTATGCTCTCCTGCGGGAggg                                                                                            <  1:1881060/97‑1 (MQ=255)
            cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGt                                                                                         <  2:3836661/100‑1 (MQ=255)
            cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGAACGCCCTCCCGCAGGAGAGCGt                                                                                         <  1:2342672/100‑1 (MQ=255)
              tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                       <  1:3690068/100‑1 (MQ=255)
              tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                       <  2:2584827/100‑1 (MQ=255)
              tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                       <  1:2479098/100‑1 (MQ=38)
              tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                       <  1:657606/100‑1 (MQ=39)
               cgcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGAGACTCCGTACGCCCTCCCGCAGGAGAGCGTacg                                                                                      <  2:4182303/100‑2 (MQ=38)
                gcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGGTTCGCAAATAATAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGTGAGTactc                                                                                     >  1:4293902/1‑100 (MQ=255)
                gcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactc                                                                                     <  1:3608728/100‑1 (MQ=255)
                gcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactc                                                                                     <  2:1816328/100‑1 (MQ=255)
                 cgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                    <  1:317200/100‑4 (MQ=38)
                  gATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcac                                                                                   <  1:3720391/100‑1 (MQ=255)
                  gATGTTCTTCGCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacggag                                                                                   <  2:1707767/100‑5 (MQ=37)
                    tGTTCTTCTCTATTGCGGTTGAAAACCTCAACTTGATTCGCAACGAAAAGAATCCGCTCTGATACGCCGTACGCCCTCCCGCAGGAGAGCGTacggag                                                                                   <  1:4169755/98‑5 (MQ=25)
                     gTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                <  2:483348/100‑1 (MQ=255)
                     gTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                <  1:3490350/100‑1 (MQ=255)
                       tcttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtt                                                                              <  2:1027747/100‑1 (MQ=255)
                        cttcttTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                             >  1:3089477/1‑100 (MQ=255)
                        cttctCTATTGCGGTTGAAAGCCTCAACTTGATTATCAACAAAAAGAATACGCTCTGATACTACGTACGCTCTCCTGCGGGAGGGCGTATTCACTACttt                                                                             >  1:1807121/1‑100 (MQ=255)
                          tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTa                                                                           <  2:2045268/100‑1 (MQ=255)
                          tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTa                                                                           >  2:786884/1‑100 (MQ=255)
                           ctctATTGCGGTTGAAAGCCTCAACTTGATTGGCAACAAAGAGAATCCGCTTTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                          <  1:1986293/100‑1 (MQ=255)
                           ctctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                          <  2:3276445/100‑1 (MQ=255)
                             ctATGGCGGTTGTAAGCCTCGACGTGATTCGCAACAAAAAGAATCCGCTCTGATACGCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtttt                                                                            <  1:90438/96‑1 (MQ=255)
                              tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCt                                                                       >  1:1355322/1‑100 (MQ=255)
                               aTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGGATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTc                                                                      <  2:1360979/100‑1 (MQ=255)
                                 tGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                    <  2:305722/100‑1 (MQ=255)
                                 tGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                    >  2:54430/1‑100 (MQ=255)
                                 tGCGGTTGAAAGCCGCAACTTGATTCGCAACATAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                    <  1:1513016/100‑1 (MQ=255)
                                  gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcaca                                                                   >  2:3465399/1‑100 (MQ=255)
                                  gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTAATACTCCGTACGCTCTCCTGCGGGCGGGCGTACTCACTCCTTTTACCCTcac                                                                    >  1:2449175/1‑99 (MQ=255)
                                     gTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGtt                                                                >  1:499345/1‑100 (MQ=255)
                                       tGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCATCACAGTTTa                                                              <  2:2935439/100‑1 (MQ=255)
                                       tGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCATTTTACCCTCACAGTTTa                                                              >  1:3283386/1‑100 (MQ=255)
                                        gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                             <  1:1023442/100‑1 (MQ=255)
                                         aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                            >  2:3309463/1‑100 (MQ=255)
                                         aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                            >  1:4082515/1‑100 (MQ=255)
                                         aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                            <  1:218348/100‑1 (MQ=255)
                                         aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                            <  1:1658901/100‑1 (MQ=255)
                                              cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                       >  2:239964/1‑100 (MQ=255)
                                                cAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTcc                                                     <  2:1274732/100‑1 (MQ=255)
                                                   cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTg                                                   >  2:1604265/1‑99 (MQ=255)
                                                     tGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCtt                                                <  2:3749434/100‑1 (MQ=255)
                                                       aTTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                              >  1:4102904/1‑100 (MQ=255)
                                                         tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                            <  1:696617/100‑1 (MQ=255)
                                                         tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                            >  1:3908794/1‑100 (MQ=255)
                                                          cGCAACAAAAAGAATCCTCTCTGTTACCCCGTACGTTCACCTCCGTGACGGCGTACTCACACCGTTTGCCCTCACAGCTTAGGATTTTCCtt                                                   <  1:3058895/92‑2 (MQ=255)
                                                          cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                           <  2:1094259/100‑1 (MQ=255)
                                                          cGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                           <  2:4209054/100‑1 (MQ=255)
                                                           gatctaaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                          <  1:2811936/95‑1 (MQ=255)
                                                             aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                        <  1:785106/100‑1 (MQ=255)
                                                             aacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                        <  2:3117418/100‑1 (MQ=255)
                                                              acaaAGAGAATCCTGTCTGATGCGCCGTACGCTCCCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatat                                       <  1:3525621/100‑1 (MQ=255)
                                                               caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                      >  2:3851601/1‑100 (MQ=255)
                                                               caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                      <  1:1602101/100‑1 (MQ=255)
                                                               caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatata                                      <  1:2612855/100‑1 (MQ=255)
                                                               caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGataaa                                      <  2:2987281/100‑3 (MQ=255)
                                                                    aGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                 >  1:1942294/1‑100 (MQ=255)
                                                                    aGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                 >  1:3031126/1‑100 (MQ=255)
                                                                       aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCCTTCAAGATATATAACGTCg                              >  2:3907866/1‑100 (MQ=255)
                                                                          cGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                           <  1:63190/100‑1 (MQ=255)
                                                                           gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGttt                          <  2:2314635/100‑1 (MQ=255)
                                                                            ctctGATACTCCGTACGCTCCCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                         <  2:3492920/100‑1 (MQ=255)
                                                                             tctGATACTCCGTACGCTCTCCTGCGGGTGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtat                        <  2:3650216/100‑1 (MQ=255)
                                                                                  tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                   >  1:1436682/1‑100 (MQ=255)
                                                                                   aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                  <  2:3081092/100‑1 (MQ=255)
                                                                                    ctgcGTACGCTCTCCTGCGGGACGTCTTACTGACTCCTTTTACCCTCACAGTTTAAGATTTTCCTTCTTTCAAGATATATACCGTCGGTTTATAAACa                   <  1:3365035/95‑1 (MQ=255)
                                                                                    cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                 >  2:3219338/1‑100 (MQ=255)
                                                                                    cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                 >  2:1487467/1‑100 (MQ=255)
                                                                                      gagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt               >  2:2692226/3‑100 (MQ=38)
                                                                                      gagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt               >  2:1212813/3‑100 (MQ=38)
                                                                                       agTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta              >  2:2727470/2‑100 (MQ=38)
                                                                                        gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat             >  2:435937/1‑100 (MQ=38)
                                                                                          aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta           <  2:825721/100‑1 (MQ=255)
                                                                                          aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta           <  2:1588058/100‑1 (MQ=255)
                                                                                          aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatt            >  2:2506552/1‑99 (MQ=255)
                                                                                          aCGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta           >  2:3992155/1‑100 (MQ=255)
                                                                                            gCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc         >  2:3340501/1‑100 (MQ=255)
                                                                                             ctctCCTGCGGGAGGGCGTACTCACTCCTTTTACCCACACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa        >  2:3100233/1‑100 (MQ=255)
                                                                                             ccctcCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTACTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTAACa        >  1:3758863/3‑100 (MQ=255)
                                                                                               ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta      >  2:436838/1‑100 (MQ=255)
                                                                                               ctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta      <  1:1050426/100‑1 (MQ=255)
                                                                                                 ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATACACAGACTATTATCAtagg    >  1:529408/1‑100 (MQ=255)
                                                                                                 ccTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg    >  2:3942734/1‑100 (MQ=255)
                                                                                                  cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggt   <  1:1569090/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTA  >  NC_000913/2971052‑2971250

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