Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 245,050 C→T P242S (CCA→TCA)  yafJ → putative amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913245,0500CT92.7% 1121.4 / 39.5 357P242S (CCA→TCA) yafJputative amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (17/9);  new base T (158/173);  total (175/182)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 366 total aligned reads displayed.

GACTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACC  >  NC_000913/244960‑245140
                                                                                          |                                                                                          
gACTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCa                                                                                          <  1:4250392/93‑1 (MQ=255)
 gattCAGCTCGCAGACCACAGCGAATCATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAg                                                                                         <  2:8686345/91‑1 (MQ=255)
 aCTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAg                                                                                         <  1:5622823/93‑1 (MQ=255)
 aCTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAg                                                                                         <  2:5543676/93‑1 (MQ=255)
 aCTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAg                                                                                         <  2:6263919/93‑1 (MQ=255)
   ttCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGc                                                                                       >  2:8129630/1‑93 (MQ=255)
   ttCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGc                                                                                       >  2:9419826/1‑93 (MQ=255)
    tCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCg                                                                                      >  1:2871883/1‑93 (MQ=255)
      aGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGaa                                                                                    <  1:2325758/93‑1 (MQ=255)
      aGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGaa                                                                                    >  1:9887886/1‑93 (MQ=255)
      aGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGaa                                                                                    <  2:2937497/93‑1 (MQ=255)
      aGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGaa                                                                                    >  2:8868648/1‑93 (MQ=255)
       gCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGGGAAt                                                                                   >  1:7031102/1‑93 (MQ=255)
        cTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGGGAATg                                                                                  >  2:749212/1‑93 (MQ=255)
        cTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATg                                                                                  <  1:2178407/93‑1 (MQ=255)
        cTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATg                                                                                  <  2:1029229/93‑1 (MQ=255)
        cTCCCAGACCACACCGAATGATGTGGTCACGGTGATGGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATg                                                                                  <  2:2219168/93‑1 (MQ=255)
         tCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATgg                                                                                 <  2:9009339/93‑1 (MQ=255)
         tCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATgg                                                                                 >  1:7291394/1‑93 (MQ=255)
         tCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATgg                                                                                 >  2:1497122/1‑93 (MQ=255)
           gCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGGATGgcg                                                                               >  2:9853043/1‑93 (MQ=255)
               cacacaCCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTa                                                                           <  2:7222319/91‑1 (MQ=255)
               aCCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTa                                                                           >  1:6019424/1‑93 (MQ=255)
               aCCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTa                                                                           <  2:9558067/93‑1 (MQ=255)
                 cacaCCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTAtt                                                                         >  2:7566687/1‑93 (MQ=255)
                  acacCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTAttt                                                                        <  1:8760317/93‑1 (MQ=255)
                     ccGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTTTTTTGc                                                                     >  2:8396333/1‑93 (MQ=255)
                     ccGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGc                                                                     <  1:2002002/93‑1 (MQ=255)
                       gAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCt                                                                   <  1:6378882/93‑1 (MQ=255)
                        aaTGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTc                                                                  <  1:4781148/93‑1 (MQ=255)
                          tgatgTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCgg                                                                <  1:613615/93‑1 (MQ=255)
                          tgatgTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCgg                                                                <  2:4496889/93‑1 (MQ=255)
                           gatgTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCggg                                                               <  2:7091543/93‑1 (MQ=255)
                           gatgTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCggg                                                               <  1:7144488/93‑1 (MQ=255)
                           gatgTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCggg                                                               >  1:3544158/1‑93 (MQ=255)
                            atgTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCgggg                                                              <  1:2704216/93‑1 (MQ=255)
                              gtgGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAg                                                            <  2:2628167/93‑1 (MQ=255)
                              gtgGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAg                                                            <  2:2389942/93‑1 (MQ=255)
                               tgGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTTCCTCGGGGAGc                                                           >  2:9219647/1‑93 (MQ=255)
                               tgGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGc                                                           >  1:1653072/1‑93 (MQ=255)
                               tgGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGc                                                           >  2:657003/1‑93 (MQ=255)
                                 gTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCgt                                                         >  1:3074322/1‑93 (MQ=255)
                                   cACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCgtgt                                                       <  2:4985818/93‑1 (MQ=255)
                                       gTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGtt                                                   >  1:5493015/1‑93 (MQ=255)
                                       gTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGtt                                                   >  2:1246920/1‑93 (MQ=255)
                                        tGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTGTTTTGCCTCGGGGAGCGTGTAGttt                                                  <  2:7374704/93‑1 (MQ=255)
                                        tGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGttt                                                  <  1:4632986/93‑1 (MQ=255)
                                        tGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGttt                                                  <  2:2723573/93‑1 (MQ=255)
                                        tGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCCGGGAGCGTGTAGttt                                                  >  1:5926021/1‑93 (MQ=255)
                                         gATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTg                                                 <  1:4477787/93‑1 (MQ=255)
                                         gATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGGTTg                                                 >  2:6206690/1‑93 (MQ=255)
                                         gATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTAGTTTGCCTCGGGGAGCGTGTAGTTTg                                                 <  1:9632263/93‑1 (MQ=255)
                                          aTTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGa                                                <  1:9135695/93‑1 (MQ=255)
                                          aTTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGa                                                >  2:649223/1‑93 (MQ=255)
                                          aTTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGa                                                >  1:844431/1‑93 (MQ=255)
                                           ttGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGAt                                               >  2:2270508/1‑93 (MQ=255)
                                             gCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGc                                             <  1:3773745/93‑1 (MQ=255)
                                             gCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGc                                             >  1:7907772/1‑93 (MQ=255)
                                              cAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGcc                                            <  2:3308724/93‑1 (MQ=255)
                                               aaCACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCa                                           <  1:9374173/93‑1 (MQ=255)
                                               aaCACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCa                                           <  2:1528975/93‑1 (MQ=255)
                                                  acaGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGtt                                        <  2:520996/93‑1 (MQ=255)
                                                   caGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTtg                                       <  1:488263/93‑1 (MQ=255)
                                                    aGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTtgt                                      <  1:3611242/93‑1 (MQ=255)
                                                      ccGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTgg                                    <  1:7653112/93‑1 (MQ=255)
                                                         cTGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTg                                 <  2:6089213/93‑1 (MQ=255)
                                                          tGACGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGc                                >  1:4769479/1‑93 (MQ=255)
                                                            aCGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGcac                              <  1:8639067/93‑1 (MQ=255)
                                                            aCGGGCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGcac                              >  1:4068416/1‑93 (MQ=255)
                                                               ggCAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAAc                           >  1:1495741/1‑93 (MQ=255)
                                                                 cAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACtt                         <  2:1850661/93‑1 (MQ=255)
                                                                 cAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACtt                         <  1:5346238/93‑1 (MQ=255)
                                                                 cAATGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACtt                         <  1:8359673/93‑1 (MQ=255)
                                                                  aaTGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTc                        <  2:4777694/93‑1 (MQ=255)
                                                                   aTGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCg                       >  1:1867973/1‑93 (MQ=255)
                                                                    tGAAACCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGt                      >  1:8243659/1‑93 (MQ=255)
                                                                      aaaCCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTgg                    <  2:5466610/93‑1 (MQ=255)
                                                                      aaaCCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTgg                    >  2:4396007/1‑93 (MQ=255)
                                                                       aaCCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGc                   >  2:3594107/1‑93 (MQ=255)
                                                                        aCCTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCt                  >  2:4674811/1‑93 (MQ=255)
                                                                          cTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTa                >  1:7756358/1‑93 (MQ=255)
                                                                          cTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTa                >  2:5821855/1‑93 (MQ=255)
                                                                          cTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTa                >  1:7413955/1‑93 (MQ=255)
                                                                          cTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTa                <  1:7057830/93‑1 (MQ=255)
                                                                          cTGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTa                >  1:280687/1‑93 (MQ=255)
                                                                           tGGCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTaa               <  1:1118091/93‑1 (MQ=255)
                                                                            ggCAAAAGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGCGGCTGCACAACTTCGTGGCTTAAc              <  1:4359658/93‑1 (MQ=255)
                                                                               aaaaGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGc           <  2:4012312/93‑1 (MQ=255)
                                                                                aaaGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCt          >  1:5074858/1‑93 (MQ=255)
                                                                                aaaGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCt          >  2:5671098/1‑93 (MQ=255)
                                                                                aaaGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCt          >  2:260180/1‑93 (MQ=255)
                                                                                 aaGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCtt         <  1:6712635/93‑1 (MQ=255)
                                                                                  aGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCttt        >  1:6931758/1‑92 (MQ=255)
                                                                                  aGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCttt        >  1:10016350/1‑92 (MQ=255)
                                                                                  aGATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTg        <  1:8743458/93‑1 (MQ=255)
                                                                                   gATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTgt       >  1:2489695/1‑92 (MQ=255)
                                                                                   gATTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGc       <  2:7684177/93‑1 (MQ=255)
                                                                                    aTTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCt      <  1:5139267/93‑1 (MQ=255)
                                                                                    aTTATGTCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCt      >  2:3134602/1‑93 (MQ=255)
                                                                                        tGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGAcc  <  2:6877881/93‑1 (MQ=255)
                                                                                          |                                                                                          
GACTTCAGCTCGCAGACCACACCGAATGATGTGGTCACGGTGATTGCAACACAGCCGCTGACGGGCAATGAAACCTGGCAAAAGATTATGCCAGGCGAATGGCGCTTATTTTGCCTCGGGGAGCGTGTAGTTTGATGCCAGTTGTGGCTGCACAACTTCGTGGCTTAACGGCTTGCTGACC  >  NC_000913/244960‑245140

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