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.2% 798.3 / ‑2.6 248Q749R (CAG→CGG) pflBformate acetyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base C (127/118);  minor base A (2/0);  total (130/118)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 263 total aligned reads displayed.

TTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATC  >  NC_000913/950442‑950621
                                                                                          |                                                                                         
ttCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgc                                                                                         >  1:4776348/1‑93 (MQ=255)
 tCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAAAGTCccgct                                                                                        >  1:6489942/1‑93 (MQ=255)
   gTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgc                                                                                      <  1:1028002/93‑1 (MQ=255)
   gTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgc                                                                                      <  2:1055924/93‑1 (MQ=255)
     ggAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgctg                                                                                    >  2:3049886/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCGGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  <  2:4747437/93‑1 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  >  1:5043033/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  <  1:1854280/93‑1 (MQ=255)
        gCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTc                                                                                 >  2:3880001/1‑93 (MQ=255)
        gCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTc                                                                                 >  1:3152200/1‑93 (MQ=255)
        gCCTTTATTGGACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTc                                                                                 >  1:4075486/1‑93 (MQ=255)
         ccTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCt                                                                                >  2:5333188/1‑93 (MQ=255)
         ccTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCt                                                                                <  1:660767/93‑1 (MQ=255)
          cTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCtt                                                                               <  2:404063/93‑1 (MQ=255)
          cTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCtt                                                                               >  1:1177819/1‑93 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              >  2:4347854/1‑93 (MQ=255)
             tATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAg                                                                            <  1:1219365/93‑1 (MQ=255)
             tATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAg                                                                            <  2:6463609/93‑1 (MQ=255)
             tATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAg                                                                            >  2:5269149/1‑93 (MQ=255)
              aTTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGt                                                                           >  2:1614382/1‑93 (MQ=255)
                  tACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGc                                                                       <  1:2212552/93‑1 (MQ=255)
                  tACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGc                                                                       <  2:2541537/93‑1 (MQ=255)
                   aCGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCg                                                                      <  1:166498/93‑1 (MQ=255)
                    cGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGa                                                                     >  2:4376771/1‑93 (MQ=255)
                     gCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAg                                                                    >  1:2775235/1‑93 (MQ=255)
                     gCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAg                                                                    >  2:150033/1‑93 (MQ=255)
                     gCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAg                                                                    >  2:5096907/1‑93 (MQ=255)
                     gCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCAGCTGCTGTTCTTTAGTCAGCGAg                                                                    >  2:6531763/1‑93 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  <  1:4104282/93‑1 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  >  2:792587/1‑93 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  <  2:4597594/93‑1 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  >  2:4321992/1‑93 (MQ=255)
                         tttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGa                                                                >  2:1381518/1‑93 (MQ=255)
                            aCTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAAc                                                             >  2:6024592/1‑93 (MQ=255)
                            aCTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAAc                                                             >  2:1698961/1‑93 (MQ=255)
                            aCTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGGTGCAAc                                                             >  1:5968892/1‑93 (MQ=255)
                             cTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                            >  2:3682796/1‑93 (MQ=255)
                             cTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                                                            <  2:752992/93‑1 (MQ=255)
                               gTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                          >  2:4500420/1‑93 (MQ=255)
                                tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTAc                                                         <  1:5325220/93‑1 (MQ=255)
                                 aCGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACt                                                        <  1:3215489/93‑1 (MQ=255)
                                   gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                      >  2:1534205/1‑93 (MQ=255)
                                   gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                      <  2:5576989/93‑1 (MQ=255)
                                       tCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                                  >  2:3498170/1‑93 (MQ=255)
                                        cAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                 >  1:2268521/1‑93 (MQ=255)
                                        cAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                 >  2:6642690/1‑93 (MQ=255)
                                         aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                >  2:111448/1‑93 (MQ=255)
                                         aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCACGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                >  2:6195965/1‑93 (MQ=255)
                                          gTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCa                                               >  1:6127551/1‑93 (MQ=255)
                                          gTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCa                                               >  1:3829659/1‑93 (MQ=255)
                                           tCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAg                                              >  2:4058874/1‑93 (MQ=255)
                                           tCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAg                                              <  1:4500097/93‑1 (MQ=255)
                                            cAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGa                                             <  1:5597152/93‑1 (MQ=255)
                                              aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                           >  2:4282813/1‑93 (MQ=255)
                                              aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa                                           <  1:1587225/93‑1 (MQ=255)
                                               aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATac                                          <  1:3545971/93‑1 (MQ=255)
                                                tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                         <  2:5813498/93‑1 (MQ=255)
                                                tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                         >  2:2381752/1‑93 (MQ=255)
                                                 cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATacac                                        <  2:5212144/93‑1 (MQ=255)
                                                  tAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                       >  2:291251/1‑93 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     <  2:2750807/93‑1 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     <  2:5923300/93‑1 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTATTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    <  2:6556074/93‑1 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    >  1:3070104/1‑93 (MQ=255)
                                                       aCATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg                                  >  1:5456128/1‑93 (MQ=255)
                                                        cATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                 >  1:1601936/1‑93 (MQ=255)
                                                        cATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                 >  2:5126288/1‑93 (MQ=255)
                                                          tAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                               >  1:4948378/1‑93 (MQ=255)
                                                            gATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGc                             >  1:3012376/1‑93 (MQ=255)
                                                            gATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGc                             >  2:3667591/1‑93 (MQ=255)
                                                            gATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGc                             >  2:4220437/1‑93 (MQ=255)
                                                             aTTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                            >  1:4188139/1‑93 (MQ=255)
                                                             aTTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                            <  1:5813958/93‑1 (MQ=255)
                                                             aTTGAGTGAAGGTACGAGTAATAACGTCCAGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCt                            >  1:4402563/1‑93 (MQ=255)
                                                               tGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGc                          <  1:1956723/93‑1 (MQ=255)
                                                                gAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                         <  1:1390298/93‑1 (MQ=255)
                                                                 aGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCgg                        <  1:2360477/93‑1 (MQ=255)
                                                                  gTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                       >  2:3622116/1‑93 (MQ=255)
                                                                    gAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGata                     >  1:2450187/1‑93 (MQ=255)
                                                                    gAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGata                     <  1:4370117/93‑1 (MQ=255)
                                                                    gAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGata                     >  2:2683437/1‑93 (MQ=255)
                                                                    gAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGata                     >  1:6698302/1‑93 (MQ=255)
                                                                      aGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtt                   <  2:961625/93‑1 (MQ=255)
                                                                      aGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtt                   >  1:6202677/1‑93 (MQ=255)
                                                                          aCGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc               >  2:2257842/1‑93 (MQ=255)
                                                                           cGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATCTTTTTcc              >  1:4442539/1‑93 (MQ=255)
                                                                             aGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCgg            >  1:3390332/1‑93 (MQ=255)
                                                                              gTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCggg           >  2:334864/1‑93 (MQ=255)
                                                                                aataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt         <  2:5311026/93‑1 (MQ=255)
                                                                                aataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt         <  1:5387522/93‑1 (MQ=255)
                                                                                 ataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGttt        <  1:4887981/93‑1 (MQ=255)
                                                                                  taaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtttt       <  2:3162429/93‑1 (MQ=255)
                                                                                  taaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtttt       >  1:2908285/1‑93 (MQ=255)
                                                                                  taaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtttt       <  2:1242700/93‑1 (MQ=255)
                                                                                   aaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc      >  1:6718212/1‑93 (MQ=255)
                                                                                   aaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc      <  2:1845468/93‑1 (MQ=255)
                                                                                   aaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc      >  2:5397278/1‑93 (MQ=255)
                                                                                    aCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc     <  2:4770705/93‑1 (MQ=255)
                                                                                      gTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt   >  1:2106433/1‑93 (MQ=255)
                                                                                       tCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATc  >  1:1060322/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATC  >  NC_000913/950442‑950621

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