Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 950,532 T→C Q749R (CAG→CGG)  pflB ← formate acetyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913950,5320TC99.1% 1093.5 / ‑1.7 343Q749R (CAG→CGG) pflBformate acetyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (160/180);  total (160/183)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 361 total aligned reads displayed.

TTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATC  >  NC_000913/950442‑950621
                                                                                          |                                                                                         
 tCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgct                                                                                        <  2:4805094/93‑1 (MQ=255)
  cGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctg                                                                                       <  2:2899250/93‑1 (MQ=255)
    tGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgct                                                                                     >  2:3614576/1‑93 (MQ=255)
    tGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgct                                                                                     >  1:6657071/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACTATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  >  2:1396461/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  >  1:3998798/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  <  2:10014347/93‑1 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  >  2:7661693/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  >  2:2131658/1‑93 (MQ=255)
        gCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTc                                                                                 >  1:2536936/1‑93 (MQ=255)
         ccTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCt                                                                                >  1:7453813/1‑93 (MQ=255)
         ccTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCt                                                                                <  2:5423117/93‑1 (MQ=255)
          cTTTATTGTACTCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCtt                                                                               <  2:272319/93‑1 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              <  1:2834347/93‑1 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              >  2:5842688/1‑93 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              >  1:9405142/1‑93 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              <  1:9352070/93‑1 (MQ=255)
            ttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTa                                                                             >  1:8367577/1‑93 (MQ=255)
            ttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTa                                                                             <  1:5781343/93‑1 (MQ=255)
              aTTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGt                                                                           <  1:8849251/93‑1 (MQ=255)
              aTTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGt                                                                           <  2:6702417/93‑1 (MQ=255)
               ttGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTc                                                                          >  2:487835/1‑93 (MQ=255)
                tGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCa                                                                         >  1:5340064/1‑93 (MQ=255)
                tGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCa                                                                         >  1:5834132/1‑93 (MQ=255)
                 gTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAg                                                                        >  2:9402500/1‑93 (MQ=255)
                  tACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGc                                                                       <  1:2236424/93‑1 (MQ=255)
                  tACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGc                                                                       >  1:704563/1‑93 (MQ=255)
                      cTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGt                                                                   <  2:442501/93‑1 (MQ=255)
                         tttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGa                                                                <  2:3889076/93‑1 (MQ=255)
                           tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                              <  2:878150/93‑1 (MQ=255)
                           tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                              <  2:1191684/93‑1 (MQ=255)
                           tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                              <  1:7990674/93‑1 (MQ=255)
                           tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                                                              <  2:7484904/93‑1 (MQ=255)
                             cTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                            >  1:6737301/1‑93 (MQ=255)
                              tGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                           <  2:6309398/93‑1 (MQ=255)
                               gTACGATTTCAGTCGACTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                          <  1:8954420/93‑1 (MQ=255)
                               gTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATACCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                          <  1:4249254/93‑1 (MQ=255)
                               gTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                          <  1:7302840/93‑1 (MQ=255)
                                tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTAc                                                         <  1:5982978/93‑1 (MQ=255)
                                 aCGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACt                                                        <  1:7782636/93‑1 (MQ=255)
                                  agATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTg                                                       <  1:7664887/92‑1 (MQ=255)
                                   gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                      <  1:5174681/93‑1 (MQ=255)
                                    aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                                     >  2:3731192/1‑93 (MQ=255)
                                    aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                                                     >  2:6235904/1‑93 (MQ=255)
                                     tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                                    >  2:1580127/1‑93 (MQ=255)
                                     tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                                    >  1:1963571/1‑93 (MQ=255)
                                     tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGt                                                    >  1:842915/1‑93 (MQ=255)
                                        cAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                 >  1:7136885/1‑93 (MQ=255)
                                         aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                >  1:7637795/1‑93 (MQ=255)
                                          gTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAAGGTACTGCGTAGCCa                                               >  2:7843473/1‑93 (MQ=255)
                                           tCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAg                                              >  1:3161430/1‑93 (MQ=255)
                                            cAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                             <  2:8073288/93‑1 (MQ=255)
                                            cAAATCTAATTACATAGATTGAGTGAAGGTACGAGCAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                             >  1:2955450/1‑93 (MQ=255)
                                             aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGAt                                            >  2:2585146/1‑93 (MQ=255)
                                              aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                           <  2:7311409/93‑1 (MQ=255)
                                               aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATac                                          >  2:5156867/1‑93 (MQ=255)
                                               aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATac                                          <  2:3675276/93‑1 (MQ=255)
                                               aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATac                                          >  1:7574692/1‑93 (MQ=255)
                                                tatAATTAGATAGATTGAGTGAAGCTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                         <  2:3255558/91‑1 (MQ=255)
                                                 cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATacac                                        <  2:5105805/93‑1 (MQ=255)
                                                 cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATacac                                        >  2:8258561/1‑93 (MQ=255)
                                                 cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATacac                                        <  2:8597587/93‑1 (MQ=255)
                                                  tAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                       >  2:5028886/1‑93 (MQ=255)
                                                   aaTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                      <  2:1784483/93‑1 (MQ=255)
                                                   aaTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                      >  1:4359063/1‑93 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     <  1:1021535/93‑1 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     <  1:4566047/93‑1 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     <  1:3563884/93‑1 (MQ=255)
                                                      tACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATg                                   >  1:632267/1‑93 (MQ=255)
                                                      tACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATg                                   <  2:8436895/93‑1 (MQ=255)
                                                         aTAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTc                                <  2:4048076/93‑1 (MQ=255)
                                                         aTAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTc                                <  1:2466850/93‑1 (MQ=255)
                                                          tAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                               <  2:1027292/93‑1 (MQ=255)
                                                          tAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                               <  2:4376700/93‑1 (MQ=255)
                                                          tAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                               >  1:5072302/1‑93 (MQ=255)
                                                           aGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                              <  1:6107145/93‑1 (MQ=255)
                                                            gATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCATCGAGGTGAAACGTACTGCGTAGCTATATACTCGGATGGTCAGc                             >  2:8986275/1‑93 (MQ=255)
                                                            gATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGc                             <  1:1909188/93‑1 (MQ=255)
                                                              ttGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTg                           >  2:557401/1‑93 (MQ=255)
                                                               tGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGc                          <  1:5495957/93‑1 (MQ=255)
                                                                   tGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACGGCGTAGCCAGATACACGGATGGTCAGCTGCGGat                      <  1:811192/93‑1 (MQ=255)
                                                                      aGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtt                   >  2:7919765/1‑93 (MQ=255)
                                                                        gTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtttt                 >  2:91099/1‑93 (MQ=255)
                                                                        gTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtttt                 >  1:1048183/1‑93 (MQ=255)
                                                                        gTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtttt                 >  1:3653431/1‑93 (MQ=255)
                                                                         tACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                <  2:4711600/93‑1 (MQ=255)
                                                                         tACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttttt                >  1:1556155/1‑93 (MQ=255)
                                                                           cGAGTAATAACGTCCCGCTGCTTTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACCCGGATGGTCAGCTGCGGATATTTTTcc              >  1:1773801/1‑93 (MQ=25)
                                                                           cGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTcc              <  2:6814339/93‑1 (MQ=255)
                                                                             aGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg            <  2:862625/93‑1 (MQ=255)
                                                                              gTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCggg           <  1:8658454/93‑1 (MQ=255)
                                                                              gTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCggg           >  1:1845146/1‑93 (MQ=255)
                                                                               taataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGt          >  1:1290378/1‑93 (MQ=255)
                                                                                aataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt         >  2:4609886/1‑93 (MQ=255)
                                                                                 ataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt        >  1:1385752/1‑93 (MQ=255)
                                                                                   agcGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc      <  1:3862376/91‑1 (MQ=255)
                                                                                      gTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt   >  2:4527642/1‑93 (MQ=255)
                                                                                      gTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt   >  2:2234632/1‑93 (MQ=255)
                                                                                       tCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATc  <  2:1943761/93‑1 (MQ=255)
                                                                                       tCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATc  >  1:9780359/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATC  >  NC_000913/950442‑950621

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: