Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,886,722 G→A G122S (GGC→AGC)  b3708 → tryptophanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,886,7220GA99.7% 2091.9 / ‑6.1 574G122S (GGC→AGC) b3708tryptophanase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (300/269);  total (302/272)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 590 total aligned reads displayed.

GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGT  >  NC_000913/3886632‑3886811
                                                                                          |                                                                                         
ggTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGc                                                                                         <  1:5619830/93‑1 (MQ=255)
 gTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGcc                                                                                        <  2:9744674/93‑1 (MQ=255)
  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                       >  2:2609738/1‑93 (MQ=255)
  ttATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCt                                                                                       >  2:5761065/1‑93 (MQ=255)
   tATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                      >  2:128666/1‑93 (MQ=255)
   tATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                      >  2:4437677/1‑93 (MQ=255)
   tATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTg                                                                                      <  2:5060190/93‑1 (MQ=255)
    aTCAATACACCATTCCGACTCACCACGGCCGTGGCGCAGAGCAAATCTATATTCCGGGACGGATTAAAAAACGCGAGCAGGAAAAACGCCTgg                                                                                     >  1:7572499/1‑93 (MQ=255)
     tCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGa                                                                                    >  1:3970034/1‑93 (MQ=255)
     tCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGa                                                                                    <  2:7890785/93‑1 (MQ=255)
        aTACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCg                                                                                 >  1:2111794/1‑93 (MQ=255)
         tccacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgc                                                                                <  2:8301914/91‑1 (MQ=255)
           cacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcag                                                                              >  1:304673/1‑93 (MQ=255)
           cacCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcag                                                                              >  1:3647691/1‑93 (MQ=255)
            cccATGCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGAGTAAAAAACGCTAGCAGGAAAAAAGCCTGGATCgcagc                                                                             <  2:8573358/92‑1 (MQ=255)
            acCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                             <  2:6770073/93‑1 (MQ=255)
            acCATCCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagc                                                                             <  1:1940692/93‑1 (MQ=255)
             ccATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagca                                                                            >  1:3218522/1‑93 (MQ=255)
             ccATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCgcagca                                                                            >  2:7212777/1‑93 (MQ=255)
                ttCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                         >  2:2946898/1‑93 (MQ=255)
                ttCCGACTAACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCaaaa                                                                         >  2:939105/1‑93 (MQ=255)
                 tCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAt                                                                        <  1:7389369/93‑1 (MQ=255)
                 tCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAt                                                                        >  2:1559084/1‑93 (MQ=255)
                  ccGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtg                                                                       <  1:836472/93‑1 (MQ=255)
                   cGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtgg                                                                      <  1:4955683/93‑1 (MQ=255)
                    gACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggg                                                                     >  2:2210009/1‑92 (MQ=255)
                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                    <  2:6474092/93‑1 (MQ=255)
                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                    <  1:5145871/93‑1 (MQ=255)
                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                    <  2:7556719/93‑1 (MQ=255)
                     aCTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAAtggtg                                                                    <  2:4635223/93‑1 (MQ=255)
                       tCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGc                                                                  >  1:8061564/1‑93 (MQ=255)
                        cACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCg                                                                 >  1:1721923/1‑93 (MQ=255)
                         aCCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGt                                                                <  2:3118950/93‑1 (MQ=255)
                          ccAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGtt                                                               >  1:9062601/1‑93 (MQ=255)
                          ccAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGtt                                                               >  1:7934112/1‑93 (MQ=255)
                             gggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTtctc                                                            <  1:2453783/93‑1 (MQ=255)
                              ggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                                           <  1:7245328/93‑1 (MQ=255)
                              ggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                                           >  1:2764066/1‑93 (MQ=255)
                              ggCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTctct                                                           >  1:8205138/1‑93 (MQ=255)
                               gCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTa                                                          <  1:9801370/93‑1 (MQ=255)
                                ccGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTaa                                                         >  2:7416915/1‑93 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTCAc                                                        >  1:1020933/1‑93 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                        <  1:6514170/93‑1 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                        >  2:2448211/1‑93 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                        >  2:4904064/1‑93 (MQ=255)
                                 cGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAAc                                                        <  2:7619583/93‑1 (MQ=255)
                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAACTGGTGGCGTTCTCTAACTAtt                                                    <  1:1834819/93‑1 (MQ=255)
                                     gcgcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTAtt                                                    <  2:8815295/93‑1 (MQ=255)
                                       gcAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTc                                                  <  1:1078292/93‑1 (MQ=255)
                                         agagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                                                >  2:4189346/1‑93 (MQ=255)
                                         agagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCtt                                                <  2:9036987/93‑1 (MQ=255)
                                          gagCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCttt                                               >  1:9542434/1‑93 (MQ=255)
                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                             >  1:5808356/1‑93 (MQ=255)
                                            gCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGa                                             >  1:626474/1‑93 (MQ=255)
                                              aaaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                                           >  1:4861429/1‑93 (MQ=255)
                                              aaaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                                           >  2:3697735/1‑93 (MQ=255)
                                              aaaTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGAATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATa                                           <  2:1037593/93‑1 (MQ=255)
                                                aTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATAcc                                         >  2:5887875/1‑93 (MQ=255)
                                                aTCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGGTCTCTAACTATTTCTTTGATAcc                                         >  2:5610752/1‑93 (MQ=255)
                                                 tgtataTTACGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                        <  1:8852670/91‑1 (MQ=255)
                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                        <  1:14530/93‑1 (MQ=255)
                                                 tCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCa                                        >  1:3436646/1‑93 (MQ=255)
                                                  cTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                                       >  1:5710619/1‑93 (MQ=255)
                                                  cTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                                       <  1:5545676/93‑1 (MQ=255)
                                                  cTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCAc                                       <  1:9977792/93‑1 (MQ=255)
                                                    atatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGc                                     <  1:7127604/93‑1 (MQ=255)
                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    <  1:9307461/93‑1 (MQ=255)
                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    <  2:8680402/93‑1 (MQ=255)
                                                     tatTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCa                                    <  2:9646045/93‑1 (MQ=255)
                                                       ttCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAgg                                  <  1:8723491/93‑1 (MQ=255)
                                                        tCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg                                 >  1:8526582/1‑93 (MQ=255)
                                                        tCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAggg                                 >  1:5898039/1‑93 (MQ=255)
                                                         ccGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGc                                >  1:3757116/1‑93 (MQ=255)
                                                          cGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGcc                               >  2:4313452/1‑93 (MQ=255)
                                                          cGGTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGcc                               >  1:9386183/1‑93 (MQ=255)
                                                           ggTACTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCa                              >  1:4438775/1‑93 (MQ=255)
                                                            gTACTGATTAAAAACCGCGAGCATGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCAt                             <  2:5531601/93‑1 (MQ=255)
                                                             gaCTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATa                            <  2:8234520/92‑1 (MQ=255)
                                                              aCTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAg                           >  1:6738812/1‑93 (MQ=255)
                                                               cTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc                          <  2:3497677/93‑1 (MQ=255)
                                                               cTGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGc                          <  1:2637483/93‑1 (MQ=255)
                                                                tGATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGcc                         >  2:3886761/1‑93 (MQ=255)
                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa                        >  1:1286110/1‑93 (MQ=255)
                                                                 gATTAAAAAACGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCa                        <  1:4101829/93‑1 (MQ=255)
                                                                       aaaaCGCGAGCAGGAAAAAAGGCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCaa                  <  1:768898/93‑1 (MQ=255)
                                                                       aaaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCaa                  <  2:6638220/93‑1 (MQ=255)
                                                                       aaaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCaa                  <  1:2557935/93‑1 (MQ=255)
                                                                        aaaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAAc                 >  1:6060176/1‑93 (MQ=255)
                                                                        aaaCGCAAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAAc                 <  2:527170/93‑1 (MQ=255)
                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACg                >  2:6222684/1‑93 (MQ=255)
                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACg                >  1:6900476/1‑93 (MQ=255)
                                                                         aaCGCGAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACg                <  2:7033101/93‑1 (MQ=255)
                                                                           cgcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGc              >  2:9478331/1‑93 (MQ=255)
                                                                            gcgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCt             <  1:5401090/93‑1 (MQ=255)
                                                                             cgAGCAGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTg            <  2:8373690/93‑1 (MQ=255)
                                                                                  aGGAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCg       >  2:8104026/1‑93 (MQ=255)
                                                                                   ggAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGt      >  1:6582858/1‑93 (MQ=255)
                                                                                   ggAAAAAAGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGt      >  2:1920783/1‑93 (MQ=255)
                                                                                       aaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTGTTTGATACCAGGCAGGGCGCTAGCCAGATAAACGGCTGTAgcgtgcgt  >  1:4282425/1‑93 (MQ=255)
                                                                                       aaaaGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGt  >  2:3373618/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGTTATCAATACACCATTCCGACTCACCAGGGCCGTGGCGCAGAGCAAATCTATATTCCGGTACTGATTAAAAAACGCGAGCAGGAAAAAGGCCTGGATCGCAGCAAAATGGTGGCGTTCTCTAACTATTTCTTTGATACCACGCAGGGCCATAGCCAGATCAACGGCTGTACCGTGCGT  >  NC_000913/3886632‑3886811

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: