Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,305,029 G→A G263G (GGC→GGT mtr ← tryptophan:H(+) symporter Mtr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,305,0290GA98.6% 1059.9 / NA 277G263G (GGC→GGTmtrtryptophan:H(+) symporter Mtr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (157/116);  minor base C (0/4);  total (157/120)
Fisher's exact test for biased strand distribution p-value = 3.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

CTCGCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCGCCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGCGTACCGT  >  NC_000913/3304931‑3305126
                                                                                                  |                                                                                                 
cTCGCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTaccacc                                                                                                 >  1:5216366/1‑101 (MQ=255)
cTCGCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTaccacc                                                                                                 <  1:5215967/101‑1 (MQ=255)
   gCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCtt                                                                                              <  2:5215970/101‑1 (MQ=255)
   gCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCtt                                                                                              >  1:5216369/1‑101 (MQ=255)
   gCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCtt                                                                                              <  2:5215990/101‑1 (MQ=255)
   gCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCtt                                                                                              <  2:5215878/101‑1 (MQ=255)
   gCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCCTAACGCCTGTACCAGCACATCAATATTACCACCCtt                                                                                              <  1:5215979/101‑1 (MQ=255)
    cTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTtc                                                                                             >  2:5216371/1‑101 (MQ=255)
     tACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTtct                                                                                            <  2:5215959/101‑1 (MQ=255)
     tACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTtct                                                                                            >  1:5216375/1‑101 (MQ=255)
      aCCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTtctc                                                                                           <  1:5215931/101‑1 (MQ=255)
       ccGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTctct                                                                                          <  2:5216102/101‑1 (MQ=255)
        cGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTg                                                                                         <  1:5216042/101‑1 (MQ=255)
         gCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGc                                                                                        <  1:5215854/101‑1 (MQ=255)
         gCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGc                                                                                        >  1:5216377/1‑101 (MQ=255)
         gCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGc                                                                                        <  2:5215898/101‑1 (MQ=255)
            aaGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAAt                                                                                     >  1:5216382/1‑101 (MQ=255)
            aaGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAAt                                                                                     >  2:5216381/1‑101 (MQ=255)
             aGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATa                                                                                    >  1:5216386/1‑101 (MQ=255)
             aGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATa                                                                                    <  2:5216195/101‑1 (MQ=255)
             aGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATa                                                                                    <  1:5215669/101‑1 (MQ=255)
              gTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATAc                                                                                   >  2:5216388/1‑101 (MQ=255)
                ttGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCg                                                                                 >  1:5216393/1‑101 (MQ=255)
                 tGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGa                                                                                >  2:5216394/1‑101 (MQ=255)
                 tGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGa                                                                                <  2:5215935/101‑1 (MQ=255)
                     aaCACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATaaa                                                                            >  2:5216399/1‑101 (MQ=255)
                      acacGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAAc                                                                           >  1:5216400/1‑101 (MQ=255)
                       cacGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACt                                                                          >  1:5216402/1‑101 (MQ=255)
                        acGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTc                                                                         <  2:5215950/101‑1 (MQ=255)
                         cGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTcc                                                                        <  2:5216022/101‑1 (MQ=255)
                          gACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCg                                                                       <  2:5215980/101‑1 (MQ=255)
                          gACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCg                                                                       >  2:5216406/1‑101 (MQ=255)
                            ccAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGa                                                                     >  1:5216409/1‑101 (MQ=255)
                             cagcagcagATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGAc                                                                    >  1:5216412/1‑101 (MQ=255)
                              agcagcagATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGa                                                                     <  2:5215813/99‑1 (MQ=255)
                                cagcagATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGAcgcg                                                                 >  2:5216415/1‑101 (MQ=255)
                                 agcagATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCgg                                                                >  1:5216417/1‑101 (MQ=255)
                                   cagATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGa                                                              >  1:5216422/1‑101 (MQ=255)
                                    agATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGAt                                                             >  2:5216423/1‑101 (MQ=255)
                                    agATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGAt                                                             <  2:5216028/101‑1 (MQ=255)
                                      aTCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGAc                                                                    >  2:5216426/1‑92 (MQ=255)
                                       tCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGtt                                                          <  1:5216139/101‑1 (MQ=255)
                                       tCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGtt                                                          <  1:5215777/101‑1 (MQ=255)
                                       tCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGtt                                                          <  2:5216004/101‑1 (MQ=255)
                                        ccAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTa                                                         >  1:5216429/1‑101 (MQ=255)
                                         cAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTAc                                                        >  1:5216431/1‑101 (MQ=255)
                                          aGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTAcc                                                       <  2:5216089/101‑1 (MQ=255)
                                           gACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTAccc                                                      >  1:5216436/1‑101 (MQ=255)
                                            aCTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCa                                                     >  2:5216440/1‑101 (MQ=255)
                                             cTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAt                                                    <  2:5216059/101‑1 (MQ=255)
                                              tACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAtc                                                   <  1:5215903/101‑1 (MQ=255)
                                              tACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAtc                                                   <  2:5215883/101‑1 (MQ=255)
                                                 ggCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAtcgtc                                                >  2:5216443/1‑101 (MQ=255)
                                                 ggCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAtcgtc                                                <  1:5215876/101‑1 (MQ=255)
                                                  gCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAtcgtcg                                               <  2:5216077/101‑1 (MQ=255)
                                                    tGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGcc                                             >  2:5216446/1‑101 (MQ=255)
                                                      ttCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAg                                           <  1:5215713/101‑1 (MQ=255)
                                                      ttCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAg                                           >  1:5216450/1‑101 (MQ=255)
                                                       tCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGc                                          <  2:5215944/101‑1 (MQ=255)
                                                        cAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCa                                         >  2:5216452/1‑101 (MQ=255)
                                                          gTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAAc                                       >  1:5216456/1‑101 (MQ=255)
                                                           tACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAAcc                                      >  2:5216457/1‑101 (MQ=255)
                                                             cgccgcTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAg                                    >  2:5216461/1‑101 (MQ=255)
                                                              gccgcTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGa                                   >  1:5216462/1‑101 (MQ=255)
                                                              gccgcTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGa                                   <  1:5216038/101‑1 (MQ=255)
                                                                 gcTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCGGATgg                                >  1:5216464/1‑101 (MQ=255)
                                                                  cTTAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGt                               <  1:5216065/101‑1 (MQ=255)
                                                                   ttAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGta                              >  2:5216468/1‑101 (MQ=255)
                                                                    tAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGtat                             <  1:5216206/101‑1 (MQ=255)
                                                                    tAACGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGtat                             >  2:5216472/1‑101 (MQ=255)
                                                                     aaCGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGtata                            <  1:5216097/101‑1 (MQ=255)
                                                                     aaCGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGtata                            >  1:5216477/1‑101 (MQ=255)
                                                                     aaCGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGtata                            >  1:5216473/1‑101 (MQ=255)
                                                                      aCGCCTGTACCAGCACATCAATATTACCCCCCTTCTCTGCAATACCGAt                                                                               <  2:5216128/49‑1 (MQ=255)
                                                                       cGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAca                          >  1:5216482/1‑101 (MQ=255)
                                                                       cGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAca                          <  1:5216241/101‑1 (MQ=255)
                                                                       cGCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAca                          >  2:5216480/1‑101 (MQ=255)
                                                                        gCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcag                         <  2:5215748/101‑1 (MQ=255)
                                                                        gCCTGTACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcag                         <  1:5216318/101‑1 (MQ=255)
                                                                             tACCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcagcgcca                    >  1:5216487/1‑101 (MQ=255)
                                                                              aCCAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcagcgccag                   >  2:5216489/1‑101 (MQ=255)
                                                                               ccAGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcagcgccagc                  <  1:5216255/101‑1 (MQ=255)
                                                                                cAGCACATCAATATTATCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCAtcg                                                  >  1:5216493/1‑68 (MQ=255)
                                                                                 aGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcagcgccagcgc                >  1:5216495/1‑101 (MQ=255)
                                                                                 aGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcagcgccagcgc                >  2:5216497/1‑101 (MQ=255)
                                                                                 aGCACATCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATAcagcgccagcgc                <  1:5216001/101‑1 (MQ=255)
                                                                                   cacaTCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAgcgccagcgcca              >  2:5216501/1‑101 (MQ=255)
                                                                                   cacaTCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAgcgccagcgcca              >  1:5216500/1‑101 (MQ=255)
                                                                                    acaTCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCAt             <  2:5216066/101‑1 (MQ=255)
                                                                                    acaTCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCAt             >  2:5216506/1‑101 (MQ=255)
                                                                                       tCAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAg          >  1:5216508/1‑101 (MQ=255)
                                                                                        cAATATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGc         <  2:5216287/101‑1 (MQ=255)
                                                                                         aaTATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGcg        >  1:5216512/1‑101 (MQ=255)
                                                                                         aaTATTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAAACAGATGGTATACAGCGCCAGCGCCATCAGcg        <  1:5216021/101‑1 (MQ=255)
                                                                                          atatTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGcgt       >  2:5216514/1‑101 (MQ=255)
                                                                                          atatTACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGcgt       >  1:5216516/1‑101 (MQ=255)
                                                                                            atTACCACCCTTCTATGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGcgtac     >  2:5216519/1‑101 (MQ=255)
                                                                                              tACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGc                                              >  1:5216527/1‑58 (MQ=255)
                                                                                              tACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGcgtaccg   >  2:5216525/1‑101 (MQ=255)
                                                                                              tACCACCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGcgtaccg   >  1:5216524/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CTCGCTACCGCAAAGTTTGAGAACACGACCAGCAGCAGATCCAGACTACGGCTGTTCAGTACGCCGCTTAACGCCTGTACCAGCACATCAATATTACCGCCCTTCTCTGCAATACCGATAAACTCCGGACGCGGGATGTTACCCATCGTCGCCAGCAACCAGATGGTATACAGCGCCAGCGCCATCAGCGTACCGT  >  NC_000913/3304931‑3305126

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