Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,287,820 G→A R269C (CGC→TGC)  purU ← formyltetrahydrofolate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,8200GA98.3% 386.8 / ‑4.6 115R269C (CGC→TGC) purUformyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (62/51);  total (63/52)
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 115 total aligned reads displayed.

CTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATG  >  NC_000913/1287723‑1287910
                                                                                                 |                                                                                          
cTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACt                                                                                          <  1:683590/100‑1 (MQ=255)
  ggTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACtgt                                                                                        >  1:459083/1‑100 (MQ=255)
   gTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACtgtg                                                                                       <  2:2881083/100‑1 (MQ=255)
       cGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCa                                                                                    >  2:2062735/1‑99 (MQ=255)
       cGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAg                                                                                   >  2:1137053/1‑100 (MQ=255)
       cGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAg                                                                                   <  1:1898383/100‑1 (MQ=255)
        gTAATTTTCATCCGTGACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGt                                                                                  >  2:2156700/1‑100 (MQ=255)
        gTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTTTTTTTTTTTAAAGAATAAACGTGCGATTACCGTAAACAAAGACGCACTGTGCCAGt                                                                                  >  2:2592955/1‑100 (MQ=255)
        gTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGt                                                                                  >  2:2782834/1‑100 (MQ=255)
          aaTTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTAc                                                                                >  1:2364810/1‑100 (MQ=255)
              ttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTg                                                                            <  1:937427/100‑1 (MQ=255)
              ttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTg                                                                            <  1:171999/100‑1 (MQ=255)
                cATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGta                                                                          <  1:790987/100‑1 (MQ=255)
                 aTCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGtat                                                                         >  2:478568/1‑100 (MQ=255)
                    cGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGt                                                                      >  2:2402915/1‑100 (MQ=255)
                     gTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTg                                                                     >  1:2019309/1‑100 (MQ=255)
                      tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGc                                                                    <  2:954372/100‑1 (MQ=255)
                       aaCGGATTAAAGGTAACCAGTTATTTTTGTTTTTTTTTATAGAATAATAGTTCGATTACCGTAAACAATGAAGTACTGTGACAGTACTTTGTATATTg                                                                     >  2:2446679/1‑98 (MQ=255)
                        aCGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCAc                                                                  >  1:1053304/1‑100 (MQ=255)
                        aCGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCAc                                                                  <  2:1645436/100‑1 (MQ=255)
                        aCGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCATAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCAc                                                                  >  2:1122424/1‑100 (MQ=255)
                         cGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACg                                                                 <  2:211390/100‑1 (MQ=255)
                           gATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGAc                                                               <  2:605021/100‑1 (MQ=255)
                           gATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGAc                                                               >  1:1201373/1‑100 (MQ=255)
                           gATTAAAAGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGAc                                                               >  1:1522770/1‑100 (MQ=255)
                            aTTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACt                                                              <  1:2356/100‑1 (MQ=255)
                               aaaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTATTGAAAACAAAGACGCACTGTGCCATTACTTTGTATAGTGCACGACTTaa                                                           >  2:78234/1‑100 (MQ=255)
                               aaaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTaa                                                           >  1:2253642/1‑100 (MQ=255)
                               aaaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTaa                                                           <  1:223152/100‑1 (MQ=255)
                                aaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAg                                                          <  1:1286654/100‑1 (MQ=255)
                                aaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAg                                                          <  1:2792276/100‑1 (MQ=255)
                                aaGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAg                                                          <  1:116188/100‑1 (MQ=255)
                                  ggTAACCAGTTATTTTTTTTTTTTTTTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGAc                                                        >  2:2223371/1‑100 (MQ=255)
                                  ggTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGAc                                                        <  2:2034581/100‑1 (MQ=255)
                                  ggTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGAc                                                        >  2:384805/1‑100 (MQ=255)
                                  ggTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGCTCAATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGAc                                                        <  2:2032939/100‑1 (MQ=255)
                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGt                                                      >  1:2164793/1‑100 (MQ=255)
                                     aaCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGtt                                                     >  1:1660100/1‑100 (MQ=255)
                                         aGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtc                                                 <  2:2221105/100‑1 (MQ=255)
                                         aGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtc                                                 >  2:800044/1‑100 (MQ=255)
                                         aGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtc                                                 >  2:2667311/1‑100 (MQ=255)
                                           ttATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtctc                                               <  2:925663/100‑1 (MQ=255)
                                           ttATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtctc                                               <  2:509998/100‑1 (MQ=255)
                                            tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCg                                              <  1:2046598/100‑1 (MQ=255)
                                             aTTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGa                                             >  2:44936/1‑100 (MQ=255)
                                               ttttGCTGGCGATTACCGAATAATCGTTCGATTACCGTAAAAAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACg                                           <  1:1259568/100‑1 (MQ=255)
                                               ttttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACg                                           <  1:2782855/100‑1 (MQ=255)
                                               ttttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACg                                           <  2:578946/100‑1 (MQ=255)
                                               ttttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACg                                           <  2:1838805/100‑1 (MQ=255)
                                                 ttGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTc                                         <  1:2624625/100‑1 (MQ=255)
                                                  tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCa                                        <  1:1307990/100‑1 (MQ=255)
                                                   gCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCAc                                       >  1:829858/1‑100 (MQ=255)
                                                   gCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCAc                                       <  1:2062748/100‑1 (MQ=255)
                                                    cTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACg                                      <  2:2164780/100‑1 (MQ=255)
                                                    cTGGCGATTAAAGAATAATCGTTCGATTACCGTAAAAAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACg                                      >  1:2032185/1‑100 (MQ=255)
                                                    cTGGCGATTAAAGAATAATCGTTCGATTACCGCAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACg                                      >  1:819416/1‑100 (MQ=255)
                                                     tGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTACGACGTTTTTCTCGACGTCACGa                                     >  2:1945731/1‑100 (MQ=255)
                                                     tGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGa                                     >  2:2530155/1‑100 (MQ=255)
                                                           ttAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                               >  2:2026124/1‑100 (MQ=255)
                                                           ttAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                               <  1:1003308/100‑1 (MQ=255)
                                                             aaaGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                             <  2:784324/100‑1 (MQ=255)
                                                             aaaGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                             <  1:2144288/100‑1 (MQ=255)
                                                             aaaGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                             >  1:679508/1‑100 (MQ=255)
                                                             aaaGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgc                              >  1:2199251/1‑99 (MQ=255)
                                                             aaaGAATAATAGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACCTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                             <  2:1892871/100‑1 (MQ=255)
                                                              aaGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcgc                            >  1:608660/1‑100 (MQ=255)
                                                                gAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcat                          >  2:2654550/1‑100 (MQ=255)
                                                                gAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcat                          >  1:1571533/1‑100 (MQ=255)
                                                                 aataatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatc                         <  1:1137057/100‑1 (MQ=255)
                                                                  ataatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatc                         <  1:2554854/99‑1 (MQ=255)
                                                                  ataatCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTGTCTCGACGTCACGACCTGCGCGcatc                         <  1:1759292/99‑1 (MQ=255)
                                                                     atCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCatat                     >  1:1914741/1‑100 (MQ=255)
                                                                     atCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCatat                     >  2:145540/1‑100 (MQ=255)
                                                                       cGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCt                   >  2:37977/1‑100 (MQ=255)
                                                                        gTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCtt                  >  1:453671/1‑100 (MQ=255)
                                                                         ttCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTc                 <  1:410885/100‑1 (MQ=255)
                                                                         ttCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTc                 <  1:314442/100‑1 (MQ=255)
                                                                           cGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg               >  1:1102445/1‑100 (MQ=255)
                                                                           cGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg               <  2:438863/100‑1 (MQ=255)
                                                                           cGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg               >  2:2325625/1‑100 (MQ=255)
                                                                             aTTACCTTAAACAAAGACGCACTGTGCCAGTACTTTTTATAGTGCACGACTGAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTcc                <  2:1259415/97‑2 (MQ=255)
                                                                              ttACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtg            <  2:2353147/100‑1 (MQ=255)
                                                                               tACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgt           <  1:840943/100‑1 (MQ=255)
                                                                                aCCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg          >  2:559015/1‑100 (MQ=255)
                                                                                 ccGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg          >  1:2633682/1‑99 (MQ=255)
                                                                                 ccGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTgtgt         >  1:282811/1‑100 (MQ=255)
                                                                                 ccGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTgtgt         >  2:2115997/1‑100 (MQ=255)
                                                                                 ccGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTgtgt         >  1:593924/1‑100 (MQ=255)
                                                                                   gTAAACAGAGACGCGCTGTGCCGGGACTTTGTGTAGTGCACGACGTAAGACTTTTGTCTCTACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAg       >  2:1687333/1‑100 (MQ=255)
                                                                                   gTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAg       >  2:1710472/1‑100 (MQ=255)
                                                                                     aaaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGt     <  1:1122428/100‑1 (MQ=255)
                                                                                     aaaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGt     >  1:1850465/1‑100 (MQ=255)
                                                                                     aaaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGt     >  2:827342/1‑100 (MQ=255)
                                                                                      aaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGATGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTa    <  1:1041423/100‑1 (MQ=255)
                                                                                      aaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTa    >  2:2857226/1‑100 (MQ=255)
                                                                                      aaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTa    >  2:2098578/1‑100 (MQ=255)
                                                                                      aaCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTa    <  1:1838969/100‑1 (MQ=255)
                                                                                       aCAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTAt   <  2:1331343/100‑1 (MQ=255)
                                                                                        cAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATg  >  2:1342929/1‑100 (MQ=255)
                                                                                        cAAAGACGCACTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATg  <  1:2846281/100‑1 (MQ=255)
                                                                                                 |                                                                                          
CTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATG  >  NC_000913/1287723‑1287910

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