Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,758,181 G→C 6.6% R470P (CGC→CCC)  ompQ → outer membrane pyoverdine efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,758,1810GC 6.6% 319.0 / 8.0 100R470P (CGC→CCC) ompQouter membrane pyoverdine efflux protein
Reads supporting (aligned to +/- strand):  ref base G (46/47);  new base C (2/5);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

CGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTT  >  NC_002947/4758035‑4758326
                                                                                                                                                  |                                                                                                                                                 
cgAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTa                                                                                                                                                >  2:1440047/1‑150 (MQ=255)
cgAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTa                                                                                                                                                <  1:1922653/150‑1 (MQ=255)
 gagaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAg                                                                                                                                               <  2:1939541/150‑1 (MQ=255)
  agaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTCAGACAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                              <  2:1068430/150‑1 (MQ=255)
  agaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                              <  1:978376/150‑1 (MQ=255)
  agaCCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                              >  2:1347047/1‑150 (MQ=255)
       tttcTGACCCTGCTCGACACCCAGCGCAGCCTGTACCAGGCACATGACCACCACCCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTc                                                                                                                                         >  2:1289953/4‑150 (MQ=255)
           tGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCACGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGc                                                                                                                                     <  1:2110895/150‑1 (MQ=255)
              ccGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCGGGCTAt                                                                                                                                  <  1:202336/150‑1 (MQ=255)
               cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCTAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATc                                                                                                                                 <  2:124504/150‑1 (MQ=255)
               cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCGGGCTATc                                                                                                                                 <  1:1339972/150‑1 (MQ=255)
               cGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTATCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATc                                                                                                                                 >  1:152333/1‑150 (MQ=255)
                 tGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                               <  1:883786/150‑1 (MQ=255)
                 tGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                               <  2:140858/150‑1 (MQ=255)
                    tCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtg                                                                                                                            >  1:734083/1‑150 (MQ=255)
                      gACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCggg                                                                                                                                      <  2:1896380/138‑1 (MQ=255)
                      gACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCTGACGCCCTAGCGCGTCggg                                                                                                                                      >  1:1896380/1‑138 (MQ=255)
                      gACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACgcgc                                                                                                                                                  >  1:15119/1‑126 (MQ=255)
                       acacCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTc                                                                                                                         <  2:1738777/150‑1 (MQ=255)
                        cacCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCa                                                                                                                        <  2:45448/150‑1 (MQ=255)
                           ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAg                                                                                                                                               >  1:1527029/1‑124 (MQ=255)
                           ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAg                                                                                                                                               <  2:1527029/124‑1 (MQ=255)
                           ccAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGa                                                                                                                     >  2:2238746/1‑150 (MQ=255)
                            cAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGAc                                                                                                                    >  2:1137191/1‑150 (MQ=255)
                             aGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGAcc                                                                                                                   <  1:1188020/150‑1 (MQ=255)
                                  aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                              <  2:1580278/150‑1 (MQ=255)
                                  aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                              <  1:1347047/150‑1 (MQ=255)
                                  aGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGc                                                                                                              <  1:699677/150‑1 (MQ=255)
                                   gCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGcc                                                                                                             <  1:1838795/150‑1 (MQ=255)
                                    ccTGTGTCAGGCCCAGGACCAGCACGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCt                                                                                                            <  1:593935/150‑1 (MQ=255)
                                     cTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGc                                                                                                                                            <  2:423999/117‑1 (MQ=255)
                                     cTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGc                                                                                                                                            >  1:423999/1‑117 (MQ=255)
                                     cTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTg                                                                                                           >  1:508352/1‑150 (MQ=255)
                                      tGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGa                                                                                                          >  2:656303/1‑150 (MQ=255)
                                        tACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAcc                                                                                                        >  1:1009107/1‑150 (MQ=255)
                                           cAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                     <  2:734083/150‑1 (MQ=255)
                                             ggCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAg                                                                                                   <  1:585338/150‑1 (MQ=255)
                                                ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAATCCCTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATATGCGGTTCACGACCGCTCACTGACCTTCAgcta                                                                                                >  2:1041328/1‑147 (MQ=255)
                                                ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAcc                                                                                                        >  1:2270164/1‑142 (MQ=255)
                                                ccAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGAcc                                                                                                        <  2:2270164/142‑1 (MQ=255)
                                                    gACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGa                                                                                            <  2:2267653/150‑1 (MQ=255)
                                                       cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                               <  1:405318/112‑1 (MQ=255)
                                                       cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTg                                                                                                                               >  2:405318/1‑112 (MQ=255)
                                                       cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCa                                                                                         >  2:433705/1‑150 (MQ=255)
                                                       cagcagGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCa                                                                                         <  2:422981/150‑1 (MQ=255)
                                                             gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                        >  2:1135864/1‑97 (MQ=255)
                                                             gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCa                                                                                                                                        <  1:1135864/97‑1 (MQ=255)
                                                             gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                     >  2:2113446/1‑132 (MQ=255)
                                                             gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTc                                                                                                     <  1:2113446/132‑1 (MQ=255)
                                                             gCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATGCCAg                                                                                   >  1:1322364/1‑150 (MQ=255)
                                                               gggCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAg                                                                                           >  2:2170632/3‑140 (MQ=255)
                                                               gggCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAg                                                                                           <  1:2170632/138‑1 (MQ=255)
                                                               ccGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAggg                                                                                 >  1:363250/1‑150 (MQ=255)
                                                                cGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAgggg                                                                                >  2:1536424/1‑150 (MQ=255)
                                                                cGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAgggg                                                                                <  2:508352/150‑1 (MQ=255)
                                                                      cAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGgcggcgc                                                                          <  1:95440/150‑1 (MQ=255)
                                                                          tAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGGTAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCATTGgcggcgcggcc                                                                      >  1:1697168/1‑149 (MQ=255)
                                                                             gCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAgg                                                                                                                                              >  1:96247/1‑75 (MQ=255)
                                                                             gCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtg                                                                                                                            >  1:1420111/1‑93 (MQ=255)
                                                                             gCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGgtg                                                                                                                            <  2:1420111/93‑1 (MQ=255)
                                                                              ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGGTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                  >  2:343577/1‑150 (MQ=255)
                                                                              ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCTCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                  >  1:412911/1‑150 (MQ=255)
                                                                              ccAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGACGTCAGGCTATCTGGTGTTCACAACCGCTGCCTGAACTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACttt                                                                  >  1:1385606/1‑150 (MQ=255)
                                                                                aaCTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAgg                                                                                                                                      >  2:2317880/1‑80 (MQ=255)
                                                                                aaCTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAgg                                                                                                                                      <  1:2317880/80‑1 (MQ=255)
                                                                                aaCTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGc                                                                >  2:1935546/1‑150 (MQ=255)
                                                                                      aGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCa                                                          >  1:1580079/1‑150 (MQ=255)
                                                                                          cAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATGCCAGGGGCGGCGCGGCACTTTGCCTACCAGccc                                                      <  2:1322364/150‑1 (MQ=255)
                                                                                                ggCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCaan                                                <  2:1296527/150‑2 (MQ=255)
                                                                                                 gCNCTTTACAAAGCCTTGGGCGGTGGCNGGCAAGCGGATANCCCGACNCNCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                               >  1:1296527/1‑150 (MQ=255)
                                                                                                 gCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                               >  1:2287430/1‑150 (MQ=255)
                                                                                                 gCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                               <  1:1525598/150‑1 (MQ=255)
                                                                                                 gCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAg                                               <  1:1935546/150‑1 (MQ=255)
                                                                                                    cTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAgcgc                                            >  2:669825/1‑150 (MQ=255)
                                                                                                     tttACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGct                                           >  1:99000/1‑150 (MQ=255)
                                                                                                         cAAAGCCTTGGCAGCTACCTGGTAAGCGGATAGAACGACACGCGACGGTGTCAGGCGATCTGGTGGTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCGTACCAGCCCGTCAAAGCGctcctc                                       >  2:1835208/1‑150 (MQ=255)
                                                                                                          aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGGGCATACCAGGGGCGGCGCGGCACTTTGcc                                                               <  1:1973521/125‑1 (MQ=255)
                                                                                                          aaaGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGGGCATACCAGGGGCGGCGCGGCACTTTGcc                                                               >  2:1973521/1‑125 (MQ=255)
                                                                                                              ccTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGc                                  <  2:152333/150‑1 (MQ=255)
                                                                                                              ccTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTAACCTTCAGGTGGCGAGCATACCAGGAGCGGCGCGCCTGTGGTCGTTCCAGCCCGTCAAAGCGCTCCTCGTCGc                                  >  1:539513/1‑150 (MQ=255)
                                                                                                                 tGGGCGGTGGCTGCCCAGCGGATAACCCTAGGCGCTAGGGCGTCAGGCTATCTTGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCtcgtcgtaat                               <  1:1289953/150‑5 (MQ=255)
                                                                                                                    gCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCCCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTaaa                            <  1:489039/150‑1 (MQ=255)
                                                                                                                      ggTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATc                          <  2:1580079/150‑1 (MQ=255)
                                                                                                                        tggctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAg                        <  1:2238746/150‑1 (MQ=255)
                                                                                                                          gctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgc                      <  1:433705/150‑1 (MQ=255)
                                                                                                                           ctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCa                                                          <  1:1228070/113‑1 (MQ=255)
                                                                                                                           ctggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCa                                                          >  2:1228070/1‑113 (MQ=255)
                                                                                                                            tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAggg                                                                                 >  2:268792/1‑89 (MQ=255)
                                                                                                                            tggcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgcgc                    >  1:1693080/1‑150 (MQ=255)
                                                                                                                            tggcAAGCGGATAGCCCAACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGCGGCGAGCATACCAGGGGCGGCGCGGCACTTTTCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTTTAAATCAGgcgc                    >  2:1858596/1‑150 (MQ=255)
                                                                                                                              gcAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCaa                  >  2:1517781/1‑150 (MQ=255)
                                                                                                                               cAAGCGGATAGCCCGCCGGCCTAGGGGGTGAGGGGAGCGGTGGTTCACGCCCGATGCCAGACGTGTGCCAGCCTAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAg                 <  1:656303/150‑1 (MQ=255)
                                                                                                                                aaGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTc                                                   <  2:1619888/115‑1 (MQ=255)
                                                                                                                                aaGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTc                                                   >  1:1619888/1‑115 (MQ=255)
                                                                                                                                aaGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAgg                >  2:2029895/1‑150 (MQ=255)
                                                                                                                                aaGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAgg                >  1:829606/1‑150 (MQ=255)
                                                                                                                                          cctgACGCCCTAGCGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTCAATCAGgc                      >  2:1342654/4‑134 (MQ=255)
                                                                                                                                          cctgACGCCCTAGCGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGgc                      <  1:1342654/131‑1 (MQ=255)
                                                                                                                                          cccGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCAt      >  2:2099202/1‑150 (MQ=255)
                                                                                                                                              aCGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGtt  <  1:1185913/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGACCAGCAGGCCCGCCTGCAGCTAAGCCAACTGCAGGCCAGCGTGGCCCTTTACAAAGCCTTGGGCGGTGGCTGGCAAGCGGATAGCCCGACGCGCTAGGGCGTCAGGCTATCTGGTGTTCACGACCGCTGCCTGACCTTCAGGTGGCGAGCATACCAGGGGCGGCGCGGCACTTTGCCTACCAGCCCGTCAAAGCGCTCCTCGTCGCTGTAAATCAGGCGCAAGGCCAGCTTCATGGTT  >  NC_002947/4758035‑4758326

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: