Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 15.7% intergenic (+31/‑186) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2610GT15.7% 540.2 / 72.2 186intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (67/89);  new base T (16/14);  total (83/103)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4060162‑4060360
                                                                                                   |                                                                                                    
  tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCg                                                                                                    <  2:377970/100‑1 (MQ=255)
     tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                 >  1:3602395/1‑100 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                >  1:2202073/1‑100 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                >  1:1207451/1‑100 (MQ=255)
       gTAAACGTCGCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                               <  2:2635980/100‑1 (MQ=255)
        tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                              <  2:1694593/100‑1 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAggn                                                                                           <  2:1140073/100‑2 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          <  2:553040/100‑1 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          <  1:1469290/100‑1 (MQ=255)
           aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACAGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          >  1:1994503/1‑100 (MQ=255)
            cttCACCTGACGGAAAACGATCGTCGCCGCGTCAACCGTGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                          <  1:2302147/97‑1 (MQ=38)
            cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                         <  2:827315/100‑1 (MQ=255)
               cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                      <  2:3191652/100‑1 (MQ=255)
               cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                      <  2:3090349/100‑1 (MQ=255)
                 tcTGACGGAATACGATCGTCTCGTCTCCTACCTTAGAACTAAATAGGATTAATTTCGTTTTTTTGTCATAACATACCTGCCATCGCTGGCATGtgtttga                                                                                    <  1:2628024/99‑3 (MQ=255)
                 ccggACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCAGCGATGGCAGGtttttttt                                                                                   <  1:1108300/97‑1 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    <  1:1912054/100‑1 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    >  2:3183332/1‑100 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    <  2:1789672/100‑1 (MQ=255)
                  cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                   <  1:232445/100‑3 (MQ=255)
                  cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                   <  1:1624306/100‑1 (MQ=255)
                   tGACGGAAAACGATCGTCGCCGCTTCAACTGTTTATCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                  <  2:59483/100‑4 (MQ=21)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttn                                                                                  <  1:276393/100‑2 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCTAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                 <  1:3024638/100‑1 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatga                                                                                 <  1:3390481/100‑5 (MQ=38)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                 <  1:3492672/100‑1 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                 <  1:3462869/100‑1 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCCTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                <  2:2629052/100‑1 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                <  1:990286/100‑1 (MQ=255)
                      cgatctaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                               <  2:2430683/94‑1 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                               <  1:61324/100‑1 (MQ=255)
                        gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                             <  1:2669222/100‑1 (MQ=255)
                            aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGccagatcg                                                                         >  1:2430669/1‑94 (MQ=255)
                            aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                         >  2:3324197/1‑100 (MQ=255)
                             cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                        <  2:2952664/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgc                                                                     >  1:78541/1‑100 (MQ=255)
                                 cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                    <  1:1058609/100‑1 (MQ=255)
                                  gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTGTTTAGTCATAAACAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCGCTCGTAAacgc                                                                   >  2:1407513/1‑100 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcac                                                                 >  2:2867496/1‑100 (MQ=255)
                                     gccgcGCCAAGCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAAAAAAAACCTGCCAGCGATGGCAGTTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                <  2:2968804/100‑1 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                >  2:2804348/1‑100 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAAGTTCGTTTGTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCATGTCGTGAGCGcaca                                                                >  1:884796/1‑100 (MQ=255)
                                       cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                >  1:2494183/1‑98 (MQ=255)
                                       cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                              <  2:2320294/100‑1 (MQ=255)
                                       cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                              >  2:505118/1‑100 (MQ=255)
                                          gcCAACCGCGCACCGAAAGACGAGTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                >  1:3586855/1‑95 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                         <  2:1216700/100‑1 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                         >  2:1916034/1‑100 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                        <  2:1774655/100‑1 (MQ=255)
                                               ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      <  2:436566/100‑1 (MQ=255)
                                               ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      <  1:1452397/100‑1 (MQ=255)
                                                cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCTTGTTTTGAGCGGAGACTAATGCGa                                                       >  1:2673119/1‑98 (MQ=255)
                                                cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                     <  1:2066069/100‑1 (MQ=255)
                                                 gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    <  1:2970914/100‑1 (MQ=255)
                                                 gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    <  2:3433877/100‑1 (MQ=255)
                                                 gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    >  1:3109929/1‑100 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  1:1748973/1‑100 (MQ=255)
                                                   gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                  <  1:658838/100‑1 (MQ=255)
                                                    cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                 >  1:1257535/1‑100 (MQ=255)
                                                      ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                               >  2:1347312/1‑100 (MQ=255)
                                                      ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                               <  1:743147/100‑1 (MQ=255)
                                                         aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                            >  2:3503766/1‑100 (MQ=255)
                                                           aGACGATTAATTTCGTTTTTTAGTCATAAGAAACCTGCAAGCGATGGCAGGTTTTTTTTTGCCTGTTGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                          <  1:2241656/100‑1 (MQ=255)
                                                            gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTAGTGAGCGCAGACCAATGCGAGTAtctacat                                             >  1:132157/1‑92 (MQ=255)
                                                             aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                        <  2:3728979/100‑1 (MQ=255)
                                                                aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                     <  1:413669/100‑1 (MQ=255)
                                                                   aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                  <  2:1994513/100‑1 (MQ=255)
                                                                    aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                 <  2:2706381/100‑1 (MQ=255)
                                                                     tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                <  2:942279/100‑1 (MQ=255)
                                                                     tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                <  2:3103607/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                            <  2:1488005/100‑1 (MQ=255)
                                                                          ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           <  2:1876176/100‑1 (MQ=255)
                                                                          ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           <  2:3149691/100‑1 (MQ=255)
                                                                           tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtgagatcgg                          >  1:1359894/1‑93 (MQ=255)
                                                                            ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         <  2:3250113/100‑1 (MQ=255)
                                                                                 gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    >  2:723761/6‑100 (MQ=38)
                                                                                 gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    >  2:3532286/6‑100 (MQ=38)
                                                                                 gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     >  2:2087304/1‑100 (MQ=255)
                                                                                 gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    <  2:2455102/100‑1 (MQ=255)
                                                                                 gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    <  2:1870224/100‑1 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTTTTTTTGTAAGCGAAAAACAAATCTATTATATCACTTATAGGTTTTTGGGTACGTCAATaga                   >  1:1101294/1‑99 (MQ=255)
                                                                                   aaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                  >  2:1584656/4‑100 (MQ=38)
                                                                                    aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  2:1606795/3‑100 (MQ=255)
                                                                                    aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  2:2940745/3‑100 (MQ=255)
                                                                                     tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                <  1:3061821/100‑1 (MQ=255)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                   >  2:3212628/2‑97 (MQ=255)
                                                                                      aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc               >  2:3416846/1‑100 (MQ=255)
                                                                                      aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCGCCCTTATAAGTCTGTGGTTACGTCAATagagag                <  1:1140744/99‑1 (MQ=255)
                                                                                       aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt              >  2:2052796/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCAGCGCTGGCAGGTTTTTTTTTGACTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt             >  2:2692103/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt             >  2:1895573/1‑100 (MQ=255)
                                                                                         aaaCCTGCCATCGCTGGCAGGTTTTTTTTTTTTTTTTTTCAGCTCACACCAATGCGAGTATCTCCCTTATAAATGTGTCGTTACGTTAAtttaatcctat            >  1:371467/1‑90 (MQ=12)
                                                                                         aaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt            >  2:3610737/1‑100 (MQ=255)
                                                                                          aaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg           <  1:2747910/100‑1 (MQ=255)
                                                                                              tGCCAGCGATGGCAGGTTTTTTTTTTTTTTTTTTTTTTGCACACAAATGCGAGTATCTCTATTTTAAGTTTGTGGTTACGTCAATATAGAGCTTTGtta        >  1:1858272/1‑97 (MQ=255)
                                                                                              tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac       <  1:3542778/100‑1 (MQ=255)
                                                                                               gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca      <  1:3690376/100‑1 (MQ=255)
                                                                                               gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca      >  2:2542998/1‑100 (MQ=255)
                                                                                                  aGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt    <  2:2202084/99‑1 (MQ=255)
                                                                                                  aGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa  <  2:1658729/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4060162‑4060360

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: