Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,180,783 A→C 5.9% D230A (GAC→GCC)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,180,7830AC 5.9% 233.5 / 2.6 69D230A (GAC→GCC) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base A (25/39);  new base C (2/3);  total (27/42)
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

TCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGC  >  NC_000913/2180637‑2180931
                                                                                                                                                  |                                                                                                                                                    
tcGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACg                                                                                                                                                    >  2:355487/1‑149 (MQ=255)
     ccATTGACGGAGAAAGCTGGTCAGCGCTTGTTGATTCCCTGCCTTCAATTGCGTGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGCCGCCAcc                                                                                                                                               <  2:122653/149‑1 (MQ=255)
      cATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCg                                                                                                                                              >  1:17321/1‑149 (MQ=255)
       aTTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCTACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAGTGCCGACGGCACCGa                                                                                                                                             >  2:169635/1‑149 (MQ=255)
         tGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGCCGGCACCGAAt                                                                                                                                           <  4:389778/149‑1 (MQ=255)
             ggAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGc                                                                                                                                       <  2:193668/149‑1 (MQ=255)
              gagaAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGcc                                                                                                                                      <  2:273041/149‑1 (MQ=255)
                   aGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGa                                                                                                                                 >  4:227292/1‑149 (MQ=255)
                    tcTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGaa                                                                                                                                <  3:50429/148‑1 (MQ=255)
                    gCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGaa                                                                                                                                <  3:301612/149‑1 (MQ=255)
                     cTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAAc                                                                                                                               >  4:47914/1‑149 (MQ=255)
                     cTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAAc                                                                                                                               <  3:137704/149‑1 (MQ=255)
                        gTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGCCGGCACCGAATCCGCCAGCGAACAGc                                                                                                                            >  1:325566/1‑149 (MQ=255)
                           aGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCttt                                                                                                                         <  1:74372/149‑1 (MQ=255)
                            gCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTa                                                                                                                        <  1:147701/149‑1 (MQ=255)
                               aTTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCa                                                                                                                     <  4:416136/149‑1 (MQ=255)
                               aTTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCa                                                                                                                     <  2:107327/149‑1 (MQ=255)
                               aTTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCa                                                                                                                     >  4:379445/1‑149 (MQ=255)
                               aTTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCa                                                                                                                     <  4:379837/149‑1 (MQ=255)
                                       ttCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTg                                                                                                             <  4:313464/149‑1 (MQ=255)
                                       ttCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTg                                                                                                             <  2:402218/149‑1 (MQ=255)
                                       ttCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTg                                                                                                             <  2:342701/149‑1 (MQ=255)
                                           cTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTTCGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAAGAGCTTTACCAGGTCGTTGgcgc                                                                                                         >  3:347078/1‑149 (MQ=255)
                                           cTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAAc                                                                                                                               <  4:410080/127‑1 (MQ=255)
                                           cTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAAc                                                                                                                               >  3:410080/1‑127 (MQ=255)
                                           cTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGgcgc                                                                                                         >  4:354440/1‑149 (MQ=255)
                                           cTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGgcgc                                                                                                         >  4:164086/1‑149 (MQ=255)
                                           cTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGgcgc                                                                                                         >  1:135941/1‑149 (MQ=255)
                                            tGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCa                                                                                                        >  2:289683/1‑149 (MQ=255)
                                              ccTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAgg                                                                                                      >  2:462246/1‑149 (MQ=255)
                                               cTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGCCGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGt                                                                                                     <  2:325566/149‑1 (MQ=255)
                                                 tCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCATCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTAc                                                                                                   <  4:381200/149‑1 (MQ=255)
                                                          cgACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTAt                                                                                          <  2:140496/149‑1 (MQ=255)
                                                           gACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTAtt                                                                                         <  1:250437/149‑1 (MQ=255)
                                                           gACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTAtt                                                                                         <  4:465759/149‑1 (MQ=255)
                                                             cATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGa                                                                                       <  3:372664/149‑1 (MQ=255)
                                                             cATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGa                                                                                       <  3:343145/149‑1 (MQ=255)
                                                              aTGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAg                                                                                      <  2:353410/149‑1 (MQ=255)
                                                                 cacaGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGtcc                                                                                   <  2:180417/149‑1 (MQ=255)
                                                                    aGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGtccgtt                                                                                >  4:321881/1‑149 (MQ=255)
                                                                    aGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGtccgtt                                                                                >  4:183079/1‑149 (MQ=255)
                                                                      cAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGtccgttcc                                                                              <  3:47914/149‑1 (MQ=255)
                                                                       aaaaaCGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGtccgttccg                                                                             >  2:372817/1‑149 (MQ=255)
                                                                          aaCGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTTAGTCCGTTCCGTgc                                                                          <  3:43188/149‑1 (MQ=255)
                                                                          aaCGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTgc                                                                          <  4:347078/149‑1 (MQ=255)
                                                                          aaCGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTgc                                                                          <  3:379445/149‑1 (MQ=255)
                                                                               aTCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCAt                                                                     >  1:825/1‑149 (MQ=255)
                                                                               aTCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCAt                                                                     >  4:233215/1‑149 (MQ=255)
                                                                                      cGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTg                                                              <  4:472760/149‑1 (MQ=255)
                                                                                      cGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTg                                                              <  4:559391/149‑1 (MQ=255)
                                                                                        ttCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGt                                                            <  2:386873/149‑1 (MQ=255)
                                                                                        ttCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGt                                                            <  2:543437/149‑1 (MQ=255)
                                                                                             cgcCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAAc                                                       <  3:139135/149‑1 (MQ=255)
                                                                                              gcCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACt                                                      <  1:227424/149‑1 (MQ=255)
                                                                                                 tCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGc                                                   <  4:475505/149‑1 (MQ=255)
                                                                                                  cACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGca                                                  >  2:300589/1‑149 (MQ=255)
                                                                                                        cggcgcggcTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACCGGCACAAAcc                                            <  4:417396/149‑1 (MQ=255)
                                                                                                         ggcgcggcTGGAGATTGCGCTGAAAATTGTTCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATGGCGCTGGTTGAACTGGCACAAACCg                                           >  2:297933/1‑149 (MQ=255)
                                                                                                            gcggcTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGAcc                                        >  2:130048/1‑149 (MQ=255)
                                                                                                                 tGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAAt                                   <  3:2004/149‑1 (MQ=255)
                                                                                                                  ggAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATc                                  >  2:176543/1‑149 (MQ=255)
                                                                                                                    agaTTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGc                                <  3:227292/149‑1 (MQ=255)
                                                                                                                     gaTTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCt                               >  1:300018/1‑149 (MQ=255)
                                                                                                                              gAAAAATATGGGCAAAGCCGCCGGCACCGAATCCGCCAGGCAACAGCTTTTCCATGGCGTTGGCGGCGGTGCCCGCACTAGTGCGGCCGTTCCGTGCGCCATTGCGCTGGGTGGAATGGCACAAACCGACCCGAATCCCTGCGcccgcc                      >  4:342099/1‑145 (MQ=255)
                                                                                                                                aaaTTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCtg                    >  2:103672/1‑149 (MQ=255)
                                                                                                                                   ttGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTgc                 <  3:233215/149‑1 (MQ=255)
                                                                                                                                   ttGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTgc                 >  3:398216/1‑149 (MQ=255)
                                                                                                                                         gcAATGCCGACGGCACCGAATCCGCCAGCGCACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAAc           <  1:378056/149‑1 (MQ=255)
                                                                                                                                              gCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTgg      <  3:39478/149‑1 (MQ=255)
                                                                                                                                                  aCGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTggcggc  <  3:354440/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGTGCCATTGACGGAGAAAGCTGGTCAGCGATTGTTGATTCCCTGCCTTCAATTGCGCGACATGCACAGCAAAAACGCATCACTACGTTCAGCGCCTCACTGGCGGCGCGGCTGGAGATTGCGCTGAAAATTGTGCGCAATGCCGACGGCACCGAATCCGCCAGCGAACAGCTTTACCAGGTCGTTGGCGCAGGTACCAGCACTATTGAGTCCGTTCCGTGCGCCATTGCGCTGGTTGAACTGGCACAAACCGACCCGAATCGCTGCGCCGTCCTGTGCGCTAACCTTGGCGGC  >  NC_000913/2180637‑2180931

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑