Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,513 T→C 83.7% D64G (GAT→GGT)  hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,5130TC83.7% 331.9 / 46.2 135D64G (GAT→GGT) hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base T (8/14);  new base C (55/58);  total (63/72)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 135 total aligned reads displayed.

CAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCG  >  NC_000913/3520417‑3520609
                                                                                                |                                                                                                
cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTTCAGTTTGTAGTGCCTGCTTCCAGGGAACACCTg                                                                                               <  2:1558976/100‑1 (MQ=255)
cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTg                                                                                               <  1:362922/100‑1 (MQ=255)
cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTg                                                                                               <  1:398327/100‑1 (MQ=255)
cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTg                                                                                               <  2:1229592/100‑1 (MQ=255)
cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATAAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTg                                                                                               >  1:2473720/1‑100 (MQ=255)
 aGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGACCACCTGt                                                                                              >  2:542493/1‑100 (MQ=255)
 aGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGt                                                                                              >  1:2500604/1‑100 (MQ=255)
  ggCAATGGAATCCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGtt                                                                                             <  1:126163/100‑1 (MQ=255)
   gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTa                                                                                            <  2:1339993/100‑1 (MQ=255)
   gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTa                                                                                            <  1:1899722/100‑1 (MQ=255)
   gCAATGGAATGCACTGAGAGAATGTTGCCCTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTa                                                                                            <  2:1713707/100‑1 (MQ=255)
    ccagtgAATGCACCGCAAGCATGTTTCCTTCCTGCCGGGGTATCACGCCGGCGCTTTTCACTACAGTTTGTCGTGCCTGCTTCCACGGAACACCTGTTa                                                                                            <  2:2566394/94‑1 (MQ=255)
    cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAg                                                                                           <  1:2085536/100‑1 (MQ=255)
    cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAg                                                                                           >  1:1439189/1‑100 (MQ=255)
     aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGa                                                                                          <  2:2270080/100‑1 (MQ=255)
     aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTCGTAGTGCCTGCTTCCAGGGAACACCTGTTAGa                                                                                          <  1:2469360/100‑1 (MQ=255)
       tgtaaTGCACTGTGATAATGTTGCGTTCCTGCCGCGTTATCAGTCCGGCGTTTTTCACTATAGTTTGTAGTTCCTGCTTCCAGGGAACACCTGTTAg                                                                                           <  2:891389/94‑1 (MQ=255)
       tGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATg                                                                                        <  1:2549128/100‑1 (MQ=255)
       tGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATg                                                                                        <  1:1624059/100‑1 (MQ=255)
         gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTa                                                                                      <  1:965710/100‑1 (MQ=255)
          aaTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTaa                                                                                     <  1:1632617/100‑1 (MQ=255)
           aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAAc                                                                                    <  2:2549444/100‑1 (MQ=255)
           aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAAc                                                                                    <  2:853748/100‑1 (MQ=255)
           aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAAc                                                                                    >  1:1258975/1‑100 (MQ=255)
              cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGac                                                                                 >  1:1473130/1‑100 (MQ=255)
                cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGacac                                                                               >  2:1378094/1‑100 (MQ=255)
                cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGacac                                                                               >  1:2485588/1‑100 (MQ=255)
                cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGacac                                                                               >  1:2132511/1‑100 (MQ=255)
                cTGAGAGAATGTTGCCTTCCAGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGacac                                                                               >  2:1589520/1‑100 (MQ=255)
                 tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACAcc                                                                              <  1:10217/100‑1 (MQ=255)
                 tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACAcc                                                                              <  1:2579621/100‑1 (MQ=255)
                     agaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTAc                                                                          <  2:339259/100‑1 (MQ=255)
                          gTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGc                                                                       >  2:1179083/1‑98 (MQ=255)
                             gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACg                                                                  <  1:1686526/100‑1 (MQ=255)
                             gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGAc                                                                   >  2:892670/1‑99 (MQ=255)
                             gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACg                                                                  <  1:908175/100‑1 (MQ=255)
                              ccTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGt                                                                 >  1:2277350/1‑100 (MQ=255)
                              ccTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGt                                                                 >  2:2266468/1‑100 (MQ=255)
                                ttCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGtgt                                                               >  2:1577739/1‑98 (MQ=255)
                                 tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                              <  1:1504913/100‑1 (MQ=255)
                                 tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                              >  1:966383/1‑100 (MQ=255)
                                   cTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                            >  2:2312038/1‑100 (MQ=255)
                                    tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                           <  1:153129/100‑1 (MQ=255)
                                    tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                           >  1:1645628/1‑100 (MQ=255)
                                    tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                           >  1:869103/1‑100 (MQ=255)
                                     gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                          >  1:715714/1‑100 (MQ=255)
                                     gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                          >  1:548458/1‑100 (MQ=255)
                                      ccGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGac                                                         <  2:1966990/100‑1 (MQ=255)
                                       cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGaca                                                        <  1:1830193/100‑1 (MQ=255)
                                       cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGaca                                                        <  1:1586000/100‑1 (MQ=255)
                                        gcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGacac                                                       >  1:2056164/1‑100 (MQ=255)
                                         cgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACg                                                      <  1:1514353/100‑1 (MQ=255)
                                         cgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACg                                                      <  2:545909/100‑1 (MQ=255)
                                          gTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGa                                                     >  1:92638/1‑100 (MQ=255)
                                           ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                    <  1:2265648/100‑1 (MQ=255)
                                                 gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                              >  1:119663/1‑100 (MQ=255)
                                                 gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                              >  2:1903119/1‑100 (MQ=255)
                                                 gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                              <  1:440286/100‑1 (MQ=255)
                                                   ccGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCa                                            <  1:2456003/100‑1 (MQ=255)
                                                   ccGGCGCTTTTCACTACAGTTTGTAGGGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGttn                                             <  2:1540172/99‑2 (MQ=255)
                                                    cGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCaa                                           >  2:623055/1‑100 (MQ=255)
                                                     ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTGCCAGGGAACACCTGTTAGATGTAACGACATCGTACCGCTGACGTCTGGCGACACGACCAGGt                                               <  2:1271061/95‑1 (MQ=255)
                                                         cTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGCCGACACGACCAGGTTCAACTtct                                      >  2:110170/1‑100 (MQ=255)
                                                          ttttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                     <  2:926418/100‑1 (MQ=255)
                                                          ttttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                     <  2:704040/100‑1 (MQ=255)
                                                           tttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                    <  2:659982/100‑1 (MQ=255)
                                                           tttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                    <  1:1473951/100‑1 (MQ=255)
                                                             tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACCTCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAACAGGTTCAACTTCTCCTg                                  <  2:2074219/100‑1 (MQ=255)
                                                             tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTg                                  >  1:677107/1‑100 (MQ=255)
                                                              cACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGt                                 <  2:1454205/100‑1 (MQ=255)
                                                               aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATTTAACGACACGTTACGGCTAACGTCTGGCGACACGACCAGGTTCAATTtctctgtt                                 >  2:1776658/1‑95 (MQ=255)
                                                               aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGtt                                >  1:2461255/1‑100 (MQ=255)
                                                               aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGtt                                >  1:20108/1‑100 (MQ=255)
                                                                  aCAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                             <  2:1007937/100‑1 (MQ=255)
                                                                    aGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTTCTGTTcagc                            >  2:1125957/1‑99 (MQ=255)
                                                                    aGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagcc                           >  1:4968/1‑100 (MQ=255)
                                                                         gTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagccag                         >  2:1551343/1‑97 (MQ=255)
                                                                         gTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAgcgc                      >  2:1543914/1‑100 (MQ=255)
                                                                          tAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                     <  2:677103/100‑1 (MQ=255)
                                                                           aGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCt                    >  2:1376738/1‑100 (MQ=255)
                                                                            gTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTg                   <  1:225503/100‑1 (MQ=255)
                                                                            gTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTg                   <  2:1937542/100‑1 (MQ=255)
                                                                             tGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                  <  2:973217/100‑1 (MQ=255)
                                                                              gCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCATCTTCTCCTGTTCAGCCAGCGCCTGCa                 <  1:1585503/100‑1 (MQ=255)
                                                                               ccTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCaa                >  2:1151016/1‑100 (MQ=255)
                                                                               ccTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCaa                >  1:1274986/1‑100 (MQ=255)
                                                                                cTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac               <  1:481022/100‑1 (MQ=255)
                                                                                  gCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAacac             >  2:2485806/1‑100 (MQ=255)
                                                                                      ccAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACAccaga         >  1:733668/1‑97 (MQ=255)
                                                                                       cAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg        >  1:1634598/1‑100 (MQ=255)
                                                                                       cAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg        <  2:2739/100‑1 (MQ=255)
                                                                                        aGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGc       <  2:182152/100‑1 (MQ=255)
                                                                                        aGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGc       <  2:2604935/100‑1 (MQ=255)
                                                                                         gggAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt      >  1:688398/1‑100 (MQ=255)
                                                                                         gggAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt      >  1:470975/1‑100 (MQ=255)
                                                                                         gggAACACCTGTTAGATGTAACGACACCGTACCACTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt      >  1:587456/1‑100 (MQ=255)
                                                                                          ggAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTa     >  2:1181457/1‑100 (MQ=255)
                                                                                            aaCACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAcc   <  2:2104232/100‑1 (MQ=255)
                                                                                             aCATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCg  <  1:1006820/100‑1 (MQ=255)
                                                                                             aCATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCg  <  1:2164277/100‑1 (MQ=255)
                                                                                                |                                                                                                
CAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCG  >  NC_000913/3520417‑3520609

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