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.6% 803.8 / NA 249Q749R (CAG→CGG) pflBformate acetyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (130/118);  minor base A (1/0);  total (131/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATC  >  NC_000913/950442‑950621
                                                                                          |                                                                                         
ttCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgc                                                                                         <  1:5090883/93‑1 (MQ=255)
ttCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgc                                                                                         <  1:5885055/93‑1 (MQ=255)
   gTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgc                                                                                      <  1:3076806/93‑1 (MQ=255)
   gTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgc                                                                                      <  1:2525641/93‑1 (MQ=255)
    tGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgct                                                                                     >  2:4701929/1‑93 (MQ=255)
    tGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgct                                                                                     >  1:5775917/1‑93 (MQ=255)
     ggAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCccgctgctg                                                                                    >  1:1132079/1‑93 (MQ=255)
      gAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGCAATAACGTCCCGCTGCTGt                                                                                   >  1:4585535/1‑93 (MQ=255)
       aGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGtt                                                                                  >  2:2359364/1‑93 (MQ=255)
        gCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCCAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGttt                                                                                 >  1:4434445/1‑92 (MQ=255)
          cTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCtt                                                                               >  2:4388285/1‑93 (MQ=255)
          cTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCtt                                                                               >  2:3416603/1‑93 (MQ=255)
          cTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGATGTTCtt                                                                               >  1:837672/1‑93 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              >  1:2340440/1‑93 (MQ=255)
           tttATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCttt                                                                              >  1:3361875/1‑93 (MQ=255)
            ttATTGTACGCTTTTTACTGTACGATTTCAGTCAAGTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTa                                                                             <  2:3610079/93‑1 (MQ=255)
               ttGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTc                                                                          >  1:1682248/1‑93 (MQ=255)
                tGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCa                                                                         >  1:1437803/1‑93 (MQ=255)
                 gTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAg                                                                        <  2:5965547/93‑1 (MQ=255)
                  tACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGc                                                                       <  2:1854939/93‑1 (MQ=255)
                  tACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGc                                                                       <  2:1213080/93‑1 (MQ=255)
                    cGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGa                                                                     <  1:4852916/93‑1 (MQ=255)
                     gCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAg                                                                    >  1:5238703/1‑93 (MQ=255)
                     gCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAg                                                                    <  1:2711533/93‑1 (MQ=255)
                      cTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGt                                                                   >  1:2229807/1‑93 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  <  2:3052135/93‑1 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  <  2:1108423/93‑1 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  >  1:6133207/1‑93 (MQ=255)
                       tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGtt                                                                  <  1:3622531/93‑1 (MQ=255)
                        ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTg                                                                 <  1:623018/93‑1 (MQ=255)
                        ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGGTg                                                                 >  2:149506/1‑93 (MQ=255)
                            aCTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAAc                                                             >  1:2105407/1‑93 (MQ=255)
                              tGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                           >  2:4783648/1‑93 (MQ=255)
                              tGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                                                           <  2:600372/93‑1 (MQ=255)
                               gTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                                                          <  1:2627913/93‑1 (MQ=255)
                                tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTAc                                                         <  2:4291960/93‑1 (MQ=255)
                                 aCGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACt                                                        <  2:4478958/93‑1 (MQ=255)
                                  cGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTg                                                       <  1:5922928/93‑1 (MQ=255)
                                   gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                      <  2:5761701/93‑1 (MQ=255)
                                   gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGc                                                      <  1:5574461/93‑1 (MQ=255)
                                       tCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                                  <  2:1820409/93‑1 (MQ=255)
                                       tCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAg                                                  <  1:4384616/93‑1 (MQ=255)
                                        cAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                                                 >  2:5331248/1‑93 (MQ=255)
                                         aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                <  2:3895670/93‑1 (MQ=255)
                                         aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                >  1:1407763/1‑93 (MQ=255)
                                         aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGcc                                                >  1:4775270/1‑93 (MQ=255)
                                             aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGACACGTACTGCGTAGCCAGAt                                            <  2:6139322/93‑1 (MQ=255)
                                             aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGAt                                            <  1:184751/93‑1 (MQ=255)
                                              aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCCGATa                                           >  2:6177224/1‑93 (MQ=255)
                                              aaTCTAATTACATAGAATGAGTGAAGGTACGAGTAATAACGTCAAACTGCTGTTCTTTGGTCAGCGAGGCTATACGTACTGCGTAGCCAGAtt                                           >  2:1133782/1‑92 (MQ=255)
                                                tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                         <  2:2864042/93‑1 (MQ=255)
                                                tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATaca                                         >  1:5159098/1‑93 (MQ=255)
                                                  tAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                       <  2:1333082/93‑1 (MQ=255)
                                                  tAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg                                       >  2:250423/1‑93 (MQ=255)
                                                   aaTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg                                      <  1:1325244/93‑1 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     >  1:384414/1‑93 (MQ=255)
                                                    aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGa                                     >  2:1601047/1‑93 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    <  1:2944436/93‑1 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    <  2:3686535/93‑1 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    <  2:3165056/93‑1 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    <  2:5227529/93‑1 (MQ=255)
                                                     ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGAt                                    <  1:5362522/93‑1 (MQ=255)
                                                      tACATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGGGTAGCCAGATTCACGGATg                                   <  2:6153824/93‑1 (MQ=255)
                                                        cATAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGt                                 <  1:958657/93‑1 (MQ=255)
                                                         aTAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTc                                >  1:6087597/1‑93 (MQ=255)
                                                         aTAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTGAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTc                                <  2:4562093/93‑1 (MQ=255)
                                                          tAGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCa                               <  1:734339/93‑1 (MQ=255)
                                                           aGATTGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAg                              <  1:3637174/93‑1 (MQ=255)
                                                              ttGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTg                           >  1:1209589/1‑93 (MQ=255)
                                                               tGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTgg                          >  1:3957283/1‑92 (MQ=255)
                                                               tGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGc                          >  1:1880835/1‑93 (MQ=255)
                                                               tGAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGc                          <  1:5692112/93‑1 (MQ=255)
                                                                gAGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCg                         >  2:5166105/1‑93 (MQ=255)
                                                                 aGTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCgg                        >  2:771436/1‑93 (MQ=255)
                                                                  gTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                       <  2:2220273/93‑1 (MQ=255)
                                                                  gTGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGa                       <  1:4359428/93‑1 (MQ=255)
                                                                   tGAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTCAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGat                      <  2:3495079/93‑1 (MQ=255)
                                                                    gAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGata                     <  1:2062875/93‑1 (MQ=255)
                                                                    gAAGGTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGata                     <  2:418279/93‑1 (MQ=255)
                                                                       ggTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAttt                  >  2:857877/1‑93 (MQ=255)
                                                                        gTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtttt                 <  2:3991273/93‑1 (MQ=255)
                                                                        gTACGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATAtttt                 >  2:3803085/1‑93 (MQ=255)
                                                                          aCGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc               >  1:2390099/1‑93 (MQ=255)
                                                                          aCGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTc               <  2:5849376/93‑1 (MQ=255)
                                                                          aCGAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATAAACGGATGGTCAGCTGCGGAGATTTTTc               >  1:2615478/1‑93 (MQ=255)
                                                                            gAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCg             >  2:2627838/1‑93 (MQ=255)
                                                                            gAGTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCg             >  1:5436879/1‑93 (MQ=255)
                                                                              gTAATAACGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCggg           <  2:139541/93‑1 (MQ=255)
                                                                               taataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGt          >  2:4044379/1‑93 (MQ=255)
                                                                                aataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt         <  2:5580027/93‑1 (MQ=255)
                                                                                aataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGtt         <  1:5050738/93‑1 (MQ=255)
                                                                                aataaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGGTGAAACGTACTGCGTAGCCAGATACACGGATGGGCAGCTGCGGATATTTTTCCGGGtg         >  2:2437408/1‑92 (MQ=21)
                                                                                   aaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc      >  1:2252014/1‑93 (MQ=255)
                                                                                   aaCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTc      <  2:5213214/93‑1 (MQ=255)
                                                                                    aCGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTcc     <  1:2405553/93‑1 (MQ=255)
                                                                                     cGTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCa    <  1:4139928/93‑1 (MQ=255)
                                                                                      gTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt   <  1:3302065/93‑1 (MQ=255)
                                                                                      gTCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt   >  1:3808210/1‑93 (MQ=255)
                                                                                       tCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATc  >  2:5532260/1‑93 (MQ=255)
                                                                                       tCCCGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATc  <  1:5030669/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTCGTGGAGCCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATC  >  NC_000913/950442‑950621

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