Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,806,424 G→A T86T (ACG→ACA kdtA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,806,4240GA99.2% 1794.4 / ‑2.6 523T86T (ACG→ACAkdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (250/267);  total (252/271)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 542 total aligned reads displayed.

CTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCG  >  NC_000913/3806334‑3806513
                                                                                          |                                                                                         
  ccgtcGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCaa                                                                                       >  2:8408578/1‑93 (MQ=255)
       ggTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGt                                                                                  <  2:4084489/93‑1 (MQ=255)
       ggTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGt                                                                                  >  1:3995798/1‑93 (MQ=255)
        gTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGtt                                                                                 <  1:1470844/93‑1 (MQ=255)
        gTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGtt                                                                                 <  1:1752973/93‑1 (MQ=255)
         tGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTc                                                                                <  2:3898279/93‑1 (MQ=255)
          gAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCg                                                                               <  1:8284598/93‑1 (MQ=255)
            agctctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGa                                                                             <  1:7744112/91‑1 (MQ=255)
            aaCTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGa                                                                             <  1:6631122/93‑1 (MQ=255)
               tctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTTCCGATTTCCGTAACAACCTTGACACCAAACTGTTTGTAgcg                                                                          >  2:12029/1‑93 (MQ=255)
               tctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAgcg                                                                          <  1:10025306/93‑1 (MQ=255)
               tctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAgcg                                                                          <  2:2721916/93‑1 (MQ=255)
               tctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAgcg                                                                          >  1:4681736/1‑93 (MQ=255)
               tctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAgcg                                                                          >  2:5906765/1‑93 (MQ=255)
               tctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAgcg                                                                          >  1:8465072/1‑93 (MQ=255)
                ctGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAgcgc                                                                         <  2:8873066/93‑1 (MQ=255)
                 tGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGcgcg                                                                        >  2:1521301/1‑93 (MQ=255)
                  ggcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGt                                                                       >  2:23470/1‑93 (MQ=255)
                   gcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTa                                                                      >  1:3459936/1‑93 (MQ=255)
                   gcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTa                                                                      <  1:7850636/93‑1 (MQ=255)
                   gcggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTa                                                                      <  2:578589/93‑1 (MQ=255)
                    cggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTAc                                                                     <  2:6520994/93‑1 (MQ=255)
                     ggcAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACa                                                                    >  1:7125009/1‑93 (MQ=255)
                        aaTCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATc                                                                 >  2:143625/1‑93 (MQ=255)
                        aaTCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATc                                                                 <  1:6169280/93‑1 (MQ=255)
                         aTCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCg                                                                >  1:5333583/1‑93 (MQ=255)
                         aTCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCg                                                                >  1:1228170/1‑93 (MQ=255)
                          tCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCgg                                                               >  1:4954042/1‑93 (MQ=255)
                           cccGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCggt                                                              >  2:3210158/1‑92 (MQ=255)
                           cccGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGc                                                              <  1:5211258/93‑1 (MQ=255)
                            ccGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCt                                                             >  2:208573/1‑93 (MQ=255)
                             cGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCtt                                                            >  2:9580728/1‑93 (MQ=255)
                              gTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCttt                                                           <  2:4769191/93‑1 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGGCACCAACCGGTTCGGAGCGCGTACAATCGGCTTTc                                                          <  1:2201414/93‑1 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGATCCCAACCGGTTCGGAGCGCGTACAATCGGCtttg                                                          >  1:907327/1‑92 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTc                                                          >  2:8169199/1‑93 (MQ=255)
                               ttGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTc                                                          <  1:596408/93‑1 (MQ=255)
                                tGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCg                                                         <  2:7231540/93‑1 (MQ=255)
                                tGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCg                                                         >  2:1092135/1‑93 (MQ=255)
                                 ggTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCgg                                                        >  1:7993587/1‑93 (MQ=255)
                                 ggTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCgg                                                        >  2:5364549/1‑93 (MQ=255)
                                  gTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCggg                                                       <  2:7825466/93‑1 (MQ=255)
                                   tGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGa                                                      >  1:2680234/1‑93 (MQ=255)
                                   tGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGa                                                      >  1:2406009/1‑93 (MQ=255)
                                    gcgcgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGaa                                                     <  1:3804374/93‑1 (MQ=255)
                                       cgctcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGa                                                  <  1:3166557/93‑1 (MQ=255)
                                       cgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGa                                                  <  2:9753316/93‑1 (MQ=255)
                                       cgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGa                                                  <  1:6913559/93‑1 (MQ=255)
                                        gcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCCGATTTCGGGAAAGAt                                                 >  2:1924388/1‑93 (MQ=255)
                                        gcccTGCGTCATCGTTATCCTTATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAt                                                 <  1:8642428/90‑1 (MQ=255)
                                         cgcgGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATg                                                <  2:4867797/93‑1 (MQ=255)
                                         cgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATg                                                >  1:9389168/1‑93 (MQ=255)
                                          gcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGt                                               >  1:880074/1‑93 (MQ=255)
                                          gcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGt                                               >  2:8816742/1‑93 (MQ=255)
                                          gcTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGt                                               >  1:310591/1‑93 (MQ=255)
                                           cTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGtt                                              <  2:5251717/93‑1 (MQ=255)
                                           cTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATAGGCTTTCGGGAAGGATGtt                                              >  2:6335374/1‑93 (MQ=255)
                                            tGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTc                                             >  2:5079483/1‑93 (MQ=255)
                                             gCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCa                                            <  2:1582721/93‑1 (MQ=255)
                                              cGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAg                                           >  1:7353045/1‑93 (MQ=255)
                                               gTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGc                                          >  2:2609013/1‑93 (MQ=255)
                                               gTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGc                                          >  2:70306/1‑93 (MQ=255)
                                                 cATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCAc                                        <  1:5481426/93‑1 (MQ=255)
                                                   tCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGCGCGCGTACAATCGGCTTTCGGGTAGGATGTTCAGCCCGt                                      >  2:4650077/1‑93 (MQ=255)
                                                   tCGTTATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGt                                      >  2:3538941/1‑93 (MQ=255)
                                                      ttATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                   >  1:2926933/1‑93 (MQ=255)
                                                      ttATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTa                                   >  1:4370253/1‑93 (MQ=255)
                                                       tATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTAt                                  >  2:8522419/1‑93 (MQ=255)
                                                       tATCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTAt                                  >  2:4327721/1‑93 (MQ=255)
                                                        aTCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATc                                 >  2:9216719/1‑93 (MQ=255)
                                                         tCCTGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCt                                <  2:7457983/93‑1 (MQ=255)
                                                            tGATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCATTCGGGCAGGATCTTCAGCACGTTTATCTGcc                             >  2:2849178/1‑93 (MQ=255)
                                                             gATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCg                            >  1:1068699/1‑93 (MQ=255)
                                                             gATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCg                            <  1:8583318/93‑1 (MQ=255)
                                                             gATTTACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTAAGCACGTTTATCTGCCg                            >  2:2160552/1‑93 (MQ=255)
                                                              aTTTAACGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGt                           <  2:661164/93‑1 (MQ=255)
                                                               tttACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGGCGTa                          >  2:5510941/1‑93 (MQ=255)
                                                               gttACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTa                          <  2:6243028/92‑1 (MQ=255)
                                                                ttACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTAt                         <  1:6354601/93‑1 (MQ=255)
                                                                ttACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTAt                         >  1:6493216/1‑93 (MQ=255)
                                                                 tACCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATg                        >  1:9880595/1‑93 (MQ=255)
                                                                  aCCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                       >  2:9445125/1‑93 (MQ=255)
                                                                  aCCGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGa                       <  2:6805480/93‑1 (MQ=255)
                                                                   ccGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGAt                      <  1:4862494/93‑1 (MQ=255)
                                                                    cGATTACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATc                     <  1:5763011/93‑1 (MQ=255)
                                                                       ttACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGc                  >  1:1570943/1‑93 (MQ=255)
                                                                       ttACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGc                  <  1:3889609/93‑1 (MQ=255)
                                                                       ttACCGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGc                  <  1:680615/93‑1 (MQ=255)
                                                                          ccGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGccca               >  2:3041956/1‑92 (MQ=255)
                                                                           cGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGa              <  1:5863055/93‑1 (MQ=255)
                                                                           cGTAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGa              >  2:8275038/1‑93 (MQ=255)
                                                                             tAACAACCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGAGGTTCACAACGTTTATCTGCCGTATGATCTGCCCgttc            >  1:419603/1‑90 (MQ=255)
                                                                               acaacCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCa          >  1:2239246/1‑93 (MQ=255)
                                                                                caacCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCAc         <  2:8074252/93‑1 (MQ=255)
                                                                                 aacCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACt        >  2:2253088/1‑93 (MQ=255)
                                                                                  acCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTc       >  1:7679601/1‑93 (MQ=255)
                                                                                  acCATGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTc       >  2:6452562/1‑93 (MQ=255)
                                                                                     aTGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAc    >  1:951439/1‑93 (MQ=255)
                                                                                      tGACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAAcc   <  1:2079648/93‑1 (MQ=255)
                                                                                       gACACCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCg  >  2:3677987/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCG  >  NC_000913/3806334‑3806513

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