Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,828,954 T→C 16.3% S4P (TCC→CCC)  xanP → xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,828,9540TC16.3% 412.6 / 66.7 183S4P (TCC→CCC) xanPxanthine permease
Reads supporting (aligned to +/- strand):  ref base T (58/94);  new base C (12/18);  total (70/113)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 188 total aligned reads displayed.

TTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGT  >  NC_000913/3828855‑3829053
                                                                                                   |                                                                                                   
tttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGttt                                                                                                     >  1:1550695/1‑100 (MQ=255)
 ttAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTc                                                                                                    <  2:19560/100‑1 (MQ=255)
     atatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCAcc                                                                                                >  1:3064894/1‑100 (MQ=255)
      tatTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCAccc                                                                                               >  2:1153029/1‑100 (MQ=255)
       atTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCt                                                                                              <  1:4009770/100‑1 (MQ=255)
         ttAGATGAGGATGTGGCAAACGTTTGATTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCg                                                                                            <  2:826814/100‑1 (MQ=255)
          tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCTACCCTCGa                                                                                           >  2:377262/1‑100 (MQ=255)
          tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGa                                                                                           >  1:2356704/1‑100 (MQ=255)
           aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAg                                                                                          <  1:2463257/100‑1 (MQ=255)
           aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAg                                                                                          >  2:3336059/1‑100 (MQ=255)
           aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAg                                                                                          <  1:1153034/100‑1 (MQ=255)
             aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTc                                                                                        <  1:1643980/100‑1 (MQ=255)
             aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTc                                                                                        >  2:3493787/1‑100 (MQ=255)
             aTGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTc                                                                                        <  2:2471546/100‑1 (MQ=255)
               gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACTGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAg                                                                                      <  2:2229688/100‑1 (MQ=255)
               gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGTATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAg                                                                                      <  2:2977766/100‑1 (MQ=255)
               gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCa                                                                                       >  2:1235277/1‑99 (MQ=255)
                   aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGaaaa                                                                                  >  1:2674650/1‑100 (MQ=255)
                   aTGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGaaa                                                                                   <  1:2121886/99‑1 (MQ=255)
                    tgtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAAt                                                                                 <  1:2663137/100‑1 (MQ=255)
                     gtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATg                                                                                >  1:3084185/1‑100 (MQ=255)
                      tgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgc                                                                               >  2:3049495/1‑100 (MQ=255)
                      tagCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATg                                                                                <  2:4007085/97‑1 (MQ=255)
                       ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgcg                                                                              <  2:1837057/100‑1 (MQ=255)
                       ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgcg                                                                              >  1:3627736/1‑100 (MQ=255)
                       ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgcg                                                                              <  2:1838949/100‑1 (MQ=255)
                       ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                              <  1:57048/100‑1 (MQ=255)
                       ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                              >  1:1475004/1‑100 (MQ=255)
                        gCAAACGTTTGCTTTTCTCTTTCGGTCAGGTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgcgc                                                                             <  2:2388302/100‑1 (MQ=255)
                        gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATgcgc                                                                             >  1:310916/1‑100 (MQ=255)
                            aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                         <  2:3270032/100‑1 (MQ=255)
                            aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                         >  1:1299633/1‑100 (MQ=255)
                            aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                         <  1:681876/100‑1 (MQ=255)
                            aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                         <  1:3703419/100‑1 (MQ=255)
                            aCGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAAcc                                                                         <  1:1856385/100‑1 (MQ=255)
                               tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGt                                                                      <  2:2605249/100‑1 (MQ=255)
                                 tGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                    <  1:1227737/100‑1 (MQ=255)
                                 tGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTg                                                                    <  2:2096174/100‑1 (MQ=255)
                                  gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCTATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                   <  2:659875/100‑1 (MQ=255)
                                      ttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAg                                                               <  1:2643994/100‑1 (MQ=255)
                                      ttCTCTTGCGGTCGGCTAAAATGCCCGCTTTGCTGCCACGGGAGTGTTTTCGATGTCGGTTTCCGCCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAg                                                               <  2:355926/100‑1 (MQ=255)
                                         tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                            >  1:243830/1‑100 (MQ=255)
                                         tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                            <  1:4123739/100‑1 (MQ=255)
                                          ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTc                                                           >  1:2719270/1‑100 (MQ=255)
                                           tttCGGTTAGCTAAAATGCCCGCATTGCTACTACTGGCTTTGTTTCGATTTTTGTTCCCACCCTCGAGTCAGATAATGCGCAACCGGTTGCGCAGACTc                                                           <  2:2056311/99‑1 (MQ=255)
                                           tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGTGCAACCGGTTGCGCAGACTCa                                                          >  1:3225566/1‑100 (MQ=255)
                                           tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttct                                                          <  1:2189690/100‑7 (MQ=255)
                                            ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttctg                                                         <  1:3203105/100‑8 (MQ=255)
                                            ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttctg                                                         <  1:2095806/100‑8 (MQ=255)
                                            ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaa                                                         >  1:26130/1‑100 (MQ=255)
                                               ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                      <  2:2714994/100‑1 (MQ=255)
                                                gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACa                                                     >  1:3566076/1‑100 (MQ=255)
                                                gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACa                                                     >  1:3944901/1‑100 (MQ=255)
                                                  cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                   <  2:383228/100‑1 (MQ=255)
                                                   aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                                  <  1:834575/100‑1 (MQ=255)
                                                    gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGa                                                 >  1:4138503/1‑100 (MQ=255)
                                                    gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGa                                                 <  2:1628003/100‑1 (MQ=255)
                                                     cTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                                  <  2:917165/98‑1 (MQ=255)
                                                      tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                               <  1:2321238/100‑1 (MQ=255)
                                                         aaTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGa                                            >  1:3432608/1‑100 (MQ=255)
                                                            gCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                         >  1:1026679/1‑100 (MQ=255)
                                                            gCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                         <  2:4123100/100‑1 (MQ=255)
                                                            gCCCGCTTTGCTACCACCGGCTTGTTTTCCATGTCTGTTTCCACCCTCGAGGCAGAAAACGCTCAACCGGTTGCCCCGGCTCAAAAc                                                      >  2:2856293/1‑87 (MQ=255)
                                                             cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                        >  1:2511556/1‑100 (MQ=255)
                                                               cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAcc                                      <  1:3286547/100‑1 (MQ=255)
                                                                 cTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGt                                    <  1:1991055/100‑1 (MQ=255)
                                                                    tGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                                 >  2:1541624/1‑100 (MQ=255)
                                                                     gcaactaCGGGTTTGTTTTCCATTTCTGTTTCCTCACACGAGGCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCCAACTGATTTACCGTCtt                                 <  2:2443330/93‑1 (MQ=255)
                                                                     gCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTg                                <  1:3526973/100‑1 (MQ=255)
                                                                       tACCACGGGTTTGTGTTCGATGTCTGTTTCCACCCTCGAGGCAGAAAATGCTCTACCGGTTTCGCAAACTCAAAACAGCGCACTGATTTACCGCCTTGa                               <  2:1477368/99‑1 (MQ=255)
                                                                       tACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGaa                              <  1:1016848/100‑1 (MQ=255)
                                                                        aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAg                             <  2:3192733/100‑1 (MQ=255)
                                                                        aCCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAg                             <  1:1371689/100‑1 (MQ=255)
                                                                         ccACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGCCTCAAAACAGCGAACTGATTTACCGTCTTGAAGa                            >  1:1474805/1‑100 (MQ=255)
                                                                          cACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGAt                           <  2:4073331/100‑1 (MQ=255)
                                                                           aCGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                          >  2:47143/1‑100 (MQ=255)
                                                                           aCGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                          >  1:3904062/1‑100 (MQ=255)
                                                                           aCGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                          <  1:1702405/100‑1 (MQ=255)
                                                                            cGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGGAAATGCGCAACCGGTTGCGCAGACTCAAAGCAGCGGACGGATTTACCGTGTTGAAGa                            <  2:2087807/97‑1 (MQ=255)
                                                                            cGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCg                         <  1:959651/100‑1 (MQ=255)
                                                                             gggATTGTTTTCGATGTCTGTTTCCACACTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGt                        <  2:54296/100‑1 (MQ=255)
                                                                                  tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                   >  1:3845784/1‑100 (MQ=255)
                                                                                  tGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                   >  1:110594/1‑100 (MQ=255)
                                                                                  tGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgcc                   <  2:95891/100‑1 (MQ=255)
                                                                                   gTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg                  >  1:424145/1‑100 (MQ=255)
                                                                                   gTTTTCGATGTCGGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg                  >  1:2747211/1‑100 (MQ=255)
                                                                                    ttttCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACGGATTTACCGTCTTGAAGATCGTCcgccgc                 >  2:1025264/1‑100 (MQ=255)
                                                                                     tttCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCt                >  1:3775429/1‑100 (MQ=255)
                                                                                         gATGTCTGTTTCCACCCTCGCGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt            <  1:611256/100‑1 (MQ=255)
                                                                                          aTGTCTGTTTCCACCCTCGAGTCATAACATGCACAACCGGTTGCGCAGACTCAAAACAGCGAACTGATGTACCGTCTTGAATATCGTCCGCCGCTTCCt            >  2:2006161/1‑99 (MQ=255)
                                                                                          aTGTCTGTTTACACCCTCGAGTCATAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCt            <  2:1079287/99‑1 (MQ=255)
                                                                                             tCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa        >  1:1141391/1‑100 (MQ=255)
                                                                                             tCTGTTTCCACCCCCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa        >  1:2834775/1‑100 (MQ=255)
                                                                                              cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc       <  1:2637892/100‑1 (MQ=255)
                                                                                              cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc       <  2:2793123/100‑1 (MQ=255)
                                                                                              cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc       <  1:3838755/100‑1 (MQ=255)
                                                                                              cTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc       <  2:1263012/100‑1 (MQ=255)
                                                                                                gTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAccc     <  1:2172093/100‑1 (MQ=255)
                                                                                                  ttCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTg   <  2:3097540/100‑1 (MQ=255)
                                                                                                   tCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGt  <  1:3299481/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTTAAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCCTGT  >  NC_000913/3828855‑3829053

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