Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 121,004 T→C 5.5% Q183R (CAG→CGG)  aroP ← aromatic amino acid:H(+) symporter AroP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913121,0040TC 5.5% 502.6 / 12.2 161Q183R (CAG→CGG) aroParomatic amino acid:H(+) symporter AroP
Reads supporting (aligned to +/- strand):  ref base T (79/73);  new base C (6/3);  total (85/76)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 165 total aligned reads displayed.

CTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAAT  >  NC_000913/120854‑121153
                                                                                                                                                      |                                                                                                                                                     
cTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCTGATCCCACAGGTTGCTAACGGTCGCCt                                                                                                                                                       >  1:921141/1‑150 (MQ=255)
    gCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCg                                                                                                                                                    >  1:1145828/1‑150 (MQ=255)
     cTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCgg                                                                                                                                                   <  1:2243075/150‑1 (MQ=255)
     cTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCgg                                                                                                                                                   >  1:1287419/1‑150 (MQ=255)
      tCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGACGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCggg                                                                                                                                                  >  1:696001/1‑150 (MQ=255)
       ccGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGgc                                                                                                                                                 >  2:65264/1‑150 (MQ=255)
           gTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGgccgcc                                                                                                                                             <  1:937177/150‑1 (MQ=255)
              aTCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGCACGAGAACATGATAATCGCCATCCTCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtg                                                                                                                                          >  1:1749094/1‑149 (MQ=255)
              aTCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGtt                                                                                                                                          >  1:2107326/1‑150 (MQ=255)
               tCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTg                                                                                                                                         >  2:1626266/1‑150 (MQ=255)
                cAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGc                                                                                                                                        <  1:981624/150‑1 (MQ=255)
                 aGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGcc                                                                                                                                       <  2:2221932/150‑1 (MQ=255)
                 aGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGCGGCCCGCAGGCGACCGttagc                                                                                                                                       <  2:2261328/150‑4 (MQ=255)
                  gCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCCCCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCa                                                                                                                                      >  2:2221419/1‑150 (MQ=255)
                   cTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCAc                                                                                                                                     >  2:8582/1‑150 (MQ=255)
                       tgctgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCAc                                                                                                                                     <  2:1558855/146‑1 (MQ=255)
                       tgctgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCAc                                                                                                                                     >  1:1558855/1‑146 (MQ=255)
                       tgctgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGaa                                                                                                                                 <  2:1194450/150‑1 (MQ=255)
                          tgcGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATTCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATag                                                                                                                              >  1:2209082/1‑150 (MQ=255)
                               tGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagcc                                                                                                                         <  1:2313082/150‑1 (MQ=255)
                                gATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagcca                                                                                                                        >  1:800548/1‑150 (MQ=255)
                                gATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACAGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagcca                                                                                                                        >  1:1374823/1‑150 (MQ=255)
                                  tCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGcagccagc                                                                                                                      >  2:25457/1‑150 (MQ=255)
                                   ccccACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCCCCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCagccagcc                                                                                                                     >  2:941012/1‑150 (MQ=255)
                                    cccACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTACCGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccg                                                                                                                    >  2:2451637/1‑150 (MQ=255)
                                      caccaGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAgccgcc                                                                                                                  <  2:653898/150‑1 (MQ=255)
                                           gTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAga                                                                                                             <  2:702069/150‑1 (MQ=255)
                                             tCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAgatg                                                                                                           <  1:678357/150‑1 (MQ=255)
                                                   ccaccGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGatcatc                                                                                                     >  1:372752/1‑150 (MQ=255)
                                                      ccGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCt                                                                                                  >  2:1155660/1‑150 (MQ=255)
                                                       cGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTa                                                                                                 >  2:1353416/1‑150 (MQ=255)
                                                         aaCGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCt                                                                                                  <  1:1313270/147‑1 (MQ=255)
                                                         aaCGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCt                                                                                                  >  2:1313270/1‑147 (MQ=255)
                                                           cGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTaccac                                                                                             >  1:1635564/1‑150 (MQ=255)
                                                             agaACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACcg                                                                                           >  1:202167/1‑150 (MQ=255)
                                                                  aTGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATa                                                                                      <  2:456162/150‑1 (MQ=255)
                                                                  aTGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATa                                                                                      <  1:347041/150‑1 (MQ=255)
                                                                   tGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCTGAAGATGATCATCGCTACCACCGCGATaa                                                                                     >  2:1092742/1‑150 (MQ=255)
                                                                      tAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACtt                                                                                  >  1:837130/1‑150 (MQ=255)
                                                                        aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTg                                                                                >  1:783905/1‑150 (MQ=255)
                                                                        aTCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTg                                                                                <  1:240342/150‑1 (MQ=255)
                                                                         tCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGa                                                                               >  1:1966158/1‑150 (MQ=255)
                                                                           gCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATa                                                                             <  1:25457/150‑1 (MQ=255)
                                                                           gCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATa                                                                             <  1:183424/150‑1 (MQ=255)
                                                                           gCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATa                                                                             <  1:1765681/150‑1 (MQ=255)
                                                                              atcatcatcACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATg                                                                          <  2:1145828/150‑1 (MQ=255)
                                                                                 atcatcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCa                                                                       <  2:2044210/150‑1 (MQ=255)
                                                                                    atcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAAc                                                                    <  1:2379481/150‑1 (MQ=255)
                                                                                     tcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACTGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAAcc                                                                   <  2:1374823/150‑1 (MQ=255)
                                                                                     tcaCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAAcc                                                                   >  2:150235/1‑150 (MQ=255)
                                                                                      caCCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCa                                                                  >  2:406125/1‑150 (MQ=255)
                                                                                          aGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAAc                                                              >  2:141881/1‑150 (MQ=255)
                                                                                            cccGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTc                                                            <  2:800548/150‑1 (MQ=255)
                                                                                             ccGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTcc                                                           >  1:1611840/1‑150 (MQ=255)
                                                                                              cGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCa                                                          <  2:202167/150‑1 (MQ=255)
                                                                                               ggTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCAt                                                         <  1:1548689/150‑1 (MQ=255)
                                                                                                 tGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGATGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCAtct                                                       <  1:1637428/150‑1 (MQ=255)
                                                                                                    aGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGTTAACTTTGATAATGGCAAACCAGAACTCCATCTCGc                                                    <  2:2083390/150‑1 (MQ=255)
                                                                                                    aGCCGTGCGGCAGGAAACAACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGc                                                    >  1:530509/1‑150 (MQ=255)
                                                                                                      cggTGCGGCAGGAAACCACCCTGATCCCACAGGTTTCTAACGGTCGCCTGCGTGCGGCCGTTGCCACTGAATAGCAGCCAGCCTGCGAAGATGATCATCTCTACCACCGCGATAACTTTGTTATTGGCAAACCAGAACTCAATCTCGCCa                                                  <  2:30990/148‑1 (MQ=255)
                                                                                                      ccGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCa                                                  <  1:711106/150‑1 (MQ=255)
                                                                                                      ccGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCaa                                                 >  1:1005756/1‑150 (MQ=255)
                                                                                                       cGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCaa                                                 <  2:818003/150‑1 (MQ=255)
                                                                                                         ggCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAac                                               <  2:512689/149‑1 (MQ=255)
                                                                                                           cGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGTCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAacac                                             <  1:595159/150‑1 (MQ=255)
                                                                                                            ggCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACt                                            <  2:2107326/150‑1 (MQ=255)
                                                                                                             gCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACtt                                           >  1:1349893/1‑150 (MQ=255)
                                                                                                              cAGGAAACCACCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTa                                         <  2:921141/150‑1 (MQ=255)
                                                                                                                      caccCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTgg                                  <  2:1729225/150‑1 (MQ=255)
                                                                                                                      caccCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCACGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTgg                                  <  1:2036696/150‑1 (MQ=255)
                                                                                                                          cTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAg                              >  1:558482/1‑150 (MQ=255)
                                                                                                                           tGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAgg                             <  2:545424/150‑1 (MQ=255)
                                                                                                                           tGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAgg                             <  2:963695/150‑1 (MQ=255)
                                                                                                                            gATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGt                            <  2:1966158/150‑1 (MQ=255)
                                                                                                                             aTCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGtt                           <  1:872072/150‑1 (MQ=255)
                                                                                                                               cccACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGa                         <  2:2452834/150‑1 (MQ=255)
                                                                                                                                ccACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGAt                        >  2:189691/1‑150 (MQ=255)
                                                                                                                                  acaGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCTGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATgg                      <  1:1092742/150‑1 (MQ=255)
                                                                                                                                   caGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGc                     <  1:131656/150‑1 (MQ=255)
                                                                                                                                   caGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGc                     <  1:2079259/150‑1 (MQ=255)
                                                                                                                                     ggTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGTTTGATGGCGt                   <  1:1155660/150‑1 (MQ=255)
                                                                                                                                      gTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAGCACTTTAACGTTGGTCAGGTTGATGGCGtt                  <  2:2209082/150‑1 (MQ=255)
                                                                                                                                       gtggcaaCGGCGGCCCGCAGGCGACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTa                 >  2:1334344/6‑150 (MQ=255)
                                                                                                                                        tggcaaCGGCGGCCCGCAGGCGACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTaa                >  1:988264/5‑150 (MQ=255)
                                                                                                                                         gCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAAt               <  2:1633227/150‑1 (MQ=255)
                                                                                                                                          cTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATc              <  2:507551/150‑1 (MQ=255)
                                                                                                                                           tAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGtt                  >  2:1744024/1‑145 (MQ=255)
                                                                                                                                           tAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGtt                  <  1:1744024/145‑1 (MQ=255)
                                                                                                                                           tAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATca             >  1:1484636/1‑150 (MQ=255)
                                                                                                                                            aaCGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcac            >  1:2104645/1‑150 (MQ=255)
                                                                                                                                            aaCGGCGGCCCGCAGGCGACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcac            >  1:1445443/1‑150 (MQ=255)
                                                                                                                                            aaCGGCGGCCCGCAGGCGACCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcac            >  1:1455153/1‑150 (MQ=255)
                                                                                                                                             aCGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcacc           <  1:150235/150‑1 (MQ=255)
                                                                                                                                              cGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcacca          <  2:372752/150‑1 (MQ=255)
                                                                                                                                              cGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcacca          >  2:409940/1‑150 (MQ=255)
                                                                                                                                              cGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGGCAGGTTGATGGCGTTAATcacca          >  2:1643800/1‑150 (MQ=255)
                                                                                                                                              cGGTCGCCTGAGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATcacca          >  1:2244939/1‑150 (MQ=255)
                                                                                                                                                   gCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAg     >  1:1702576/1‑150 (MQ=255)
                                                                                                                                                      tGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAAt  >  1:1370049/1‑150 (MQ=255)
                                                                                                                                                      tGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAAt  <  2:1262038/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTTTGCTCCGGGTTATCAGCTTCTGCTGCGGTGATCCCCACCAGTTCCAGACCACCGAACGAGAACATGATAATCGCCATCATCATCACCAGCCCGGTGAAGCCGTGCGGCAGGAAACCACCCTGATCCCACAGGTTGCTAACGGTCGCCTGCGGGCCGCCGTTGCCACTGAATAGCAGCCAGCCGCCGAAGATGATCATCGCTACCACCGCGATAACTTTGATAATGGCAAACCAGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAAT  >  NC_000913/120854‑121153

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