Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,281,517 G→T 7.7% L71L (CTG→CTT narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,281,5170GT 7.7% 1171.6 / 40.3 324L71L (CTG→CTTnarInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base G (148/150);  new base T (17/9);  total (165/159)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

ATTGGTAGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACGCGTGCGTGCAA  >  CP009273/1281368‑1281665
                                                                                                                                                     |                                                                                                                                                    
  tGGTAGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAc                                                                                                                                                    <  1:3203688/150‑1 (MQ=255)
  tGGTAGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAc                                                                                                                                                    <  2:4879217/150‑1 (MQ=255)
   ggTAGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTAcg                                                                                                                                                   <  1:4409853/150‑1 (MQ=255)
     tAGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgcc                                                                                                                                                 >  1:955205/1‑150 (MQ=255)
      aGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                                                                                                                >  1:5007210/1‑150 (MQ=255)
      aGCTGGCTGCGTTATGACTACGGGCAGTACACCGGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGAGTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg                                                                                                                                                >  1:2788901/1‑150 (MQ=255)
         tggctgCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCAc                                                                                                                                             >  1:2540307/1‑150 (MQ=255)
                  tATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTAt                                                                                                                                    >  2:703084/1‑150 (MQ=255)
                   aTGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATg                                                                                                                                   <  1:3264721/150‑1 (MQ=255)
                   aTGACAACGGGCAGTACACCTGGCGCGCGGCGACAATCCAGAAGCTGGAACGCAAAGGGATGAACCTGGCTACGAACATGTTCCATATCGGGATTCTGGGGATTATTGTCGGTCACTACTTAGGTATGCTGACGCCTCACTGGATGTaag                                                                                                                                   >  2:2144346/1‑148 (MQ=255)
                     gACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa                                                                                                                                 >  1:4197017/1‑150 (MQ=255)
                     gACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa                                                                                                                                 >  1:2945360/1‑150 (MQ=255)
                       cTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGc                                                                                                                               <  2:445269/150‑1 (MQ=255)
                       cTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACGGGATGTATGAAGc                                                                                                                               >  2:922949/1‑150 (MQ=255)
                         aCGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCt                                                                                                                             <  2:1665354/150‑1 (MQ=255)
                         aCGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                                                                                                             >  2:3103516/1‑150 (MQ=255)
                         aCGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                                                                                                             <  1:3048151/150‑1 (MQ=255)
                         aCGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt                                                                                                                             <  2:1152157/150‑1 (MQ=255)
                          cGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg                                                                                                                            <  1:3241127/150‑1 (MQ=255)
                          cGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTg                                                                                                                            <  1:1617478/150‑1 (MQ=255)
                             gCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt                                                                                                                         <  1:733151/150‑1 (MQ=255)
                               aGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGc                                                                                                                       >  2:2321236/1‑150 (MQ=255)
                                  acacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGa                                                                                                                    >  1:1904426/1‑150 (MQ=255)
                                  acacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGAGTCTGGGGATTTTAGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCgc                                                                                                                    >  2:4172755/1‑149 (MQ=255)
                                   cacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAt                                                                                                                   <  1:1051302/150‑1 (MQ=255)
                                   cacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAt                                                                                                                   <  2:1638087/150‑1 (MQ=255)
                                      cTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa                                                                                                                >  2:424379/1‑150 (MQ=255)
                                       tGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGaa                                                                                                               <  1:4819051/150‑1 (MQ=255)
                                       tGGCGCGCGGCGTCCAGCCAGATGCGGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGGTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGaa                                                                                                               <  2:101122/150‑1 (MQ=255)
                                        ggCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAg                                                                                                              <  1:1503371/150‑1 (MQ=255)
                                         gcgcgcGGCGTCCAGCCAGATGCTGGATCGCAAAGCGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTNCCGATTGAAGt                                                                                                             <  1:2761833/150‑1 (MQ=255)
                                          cgcgcgGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTg                                                                                                            <  1:4471177/150‑1 (MQ=255)
                                             gcggcgTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa                                                                                                         >  2:1296364/1‑150 (MQ=255)
                                               ggcgTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACa                                                                                                       <  1:4001201/150‑1 (MQ=255)
                                                  gTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa                                                                                                    >  2:2743386/1‑150 (MQ=255)
                                                   tCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaa                                                                                                   >  2:3359060/1‑150 (MQ=255)
                                                     cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt                                                                                                 <  1:2133949/150‑1 (MQ=255)
                                                     cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt                                                                                                 <  1:1386243/150‑1 (MQ=255)
                                                     cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAt                                                                                                 <  1:2078833/150‑1 (MQ=255)
                                                        ccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGc                                                                                              <  1:2478864/150‑1 (MQ=255)
                                                        ccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGc                                                                                              <  1:1476686/150‑1 (MQ=255)
                                                         cagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCa                                                                                             >  2:1937342/1‑150 (MQ=255)
                                                           gATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACATAAAATGGCAAt                                                                                           >  2:3523470/1‑150 (MQ=255)
                                                            aTGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATg                                                                                          <  1:5230064/150‑1 (MQ=255)
                                                            aTGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATg                                                                                          <  1:3695823/150‑1 (MQ=255)
                                                               cTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGttt                                                                                       <  1:3174654/150‑1 (MQ=255)
                                                               cTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGttt                                                                                       <  2:3147999/150‑1 (MQ=255)
                                                                 ggATCGCAAAGGGATGAAACTGGCGTCTAACCTGTTCCATATCGGGATTCTGGTGATTTTTGTCGGTCAATTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGACTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGc                                                                                     <  1:3556746/150‑1 (MQ=255)
                                                                   aTCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAATTGGCAATGTTTGCtg                                                                                   >  2:2075151/1‑150 (MQ=255)
                                                                   aTCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg                                                                                   >  2:2930792/1‑150 (MQ=255)
                                                                      gCATAGGGATGAACCTGGCGTAGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtggtg                                                                                <  2:3579765/150‑1 (MQ=255)
                                                                      gCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtggtg                                                                                >  1:1404642/1‑150 (MQ=255)
                                                                      gCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtggtg                                                                                <  2:3076830/150‑1 (MQ=255)
                                                                       cAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtggtgg                                                                               >  2:848067/1‑150 (MQ=255)
                                                                        aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTggtggt                                                                              >  2:2654160/1‑150 (MQ=255)
                                                                        aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTggtggt                                                                              <  1:1187162/150‑1 (MQ=255)
                                                                        aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTggtggt                                                                              >  2:2944832/1‑150 (MQ=255)
                                                                        aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTggtggt                                                                              <  1:198461/150‑1 (MQ=255)
                                                                                gAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGAGTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAgcgg                                                                      >  1:1950944/1‑150 (MQ=255)
                                                                                 aaCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAgcggc                                                                     >  1:4776199/1‑150 (MQ=255)
                                                                                      ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCt                                                                <  1:2235984/150‑1 (MQ=255)
                                                                                        cGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCtgt                                                              >  2:3192919/1‑150 (MQ=255)
                                                                                         gTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCtgtg                                                             >  2:5276992/1‑150 (MQ=255)
                                                                                            gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTgtgtat                                                          <  1:2689975/150‑3 (MQ=255)
                                                                                            gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCt                                                          <  1:460766/150‑1 (MQ=255)
                                                                                             aaCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTg                                                         <  1:2205020/150‑1 (MQ=255)
                                                                                             aaCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTg                                                         >  2:2059256/1‑150 (MQ=255)
                                                                                             aaCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTg                                                         <  2:4009826/150‑1 (MQ=255)
                                                                                              aCCTGTTCCATATCGGTATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGa                                                        <  1:5130524/150‑1 (MQ=255)
                                                                                              aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGa                                                        >  2:3837733/1‑150 (MQ=255)
                                                                                                cTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGAtt                                                      >  1:3323529/1‑150 (MQ=255)
                                                                                                cTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGAtt                                                      >  2:2578141/1‑150 (MQ=255)
                                                                                                 tGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTg                                                     <  1:3835825/150‑1 (MQ=255)
                                                                                                       atatCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGgcggcg                                               >  1:2618413/1‑150 (MQ=255)
                                                                                                         atCGGGATTCTGGGGATTTTTGTCGGTCACTTCGTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTg                                             >  1:3621433/1‑150 (MQ=255)
                                                                                                          tCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGc                                            <  1:4888487/150‑1 (MQ=255)
                                                                                                          tCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGc                                            <  1:703084/150‑1 (MQ=255)
                                                                                                             ggATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGt                                         >  1:4425755/1‑150 (MQ=255)
                                                                                                                   tGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGa                                   >  1:2447782/1‑150 (MQ=255)
                                                                                                                    ggggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGaa                                  >  2:4413397/1‑150 (MQ=255)
                                                                                                                     gggATTTTTGTCGGTCACTTCTTCGGTATGCAGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGaca                                 >  2:4698470/1‑148 (MQ=255)
                                                                                                                      ggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAc                                <  1:830546/150‑1 (MQ=255)
                                                                                                                       gATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcg                               <  1:4547811/150‑1 (MQ=255)
                                                                                                                        aTTTTTGTCGGTCACTTCTTCGGTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgt                              <  2:1373386/150‑1 (MQ=255)
                                                                                                                         tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgtc                             <  1:4676881/150‑1 (MQ=255)
                                                                                                                          ttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgtcg                            <  1:180684/150‑1 (MQ=255)
                                                                                                                             tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCt                         >  2:4928013/1‑150 (MQ=255)
                                                                                                                             tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCt                         >  2:5154628/1‑150 (MQ=255)
                                                                                                                             tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCt                         >  1:792647/1‑150 (MQ=255)
                                                                                                                             tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCt                         >  1:587615/1‑150 (MQ=255)
                                                                                                                                cGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGtt                      <  1:307568/150‑1 (MQ=255)
                                                                                                                                    cACTACTTCGGTATGCTGACGCCGCACTGGAAGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAAAGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCAGTGACAGATTGGCGGCGTGCTGTTGCTGACACGACGTCTGGTCAGc                  >  2:2377250/1‑150 (MQ=255)
                                                                                                                                     aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCTGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGcc                 >  2:3391300/1‑150 (MQ=255)
                                                                                                                                       ttcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCa               >  2:5059839/1‑150 (MQ=255)
                                                                                                                                         cttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCAcg             >  1:1690294/1‑150 (MQ=255)
                                                                                                                                            cGGTATGCTTACGCCGCACTGGACGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACgcgt          >  1:609153/1‑150 (MQ=255)
                                                                                                                                              gTATGCTTACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACgcgtgc        >  2:250035/1‑150 (MQ=255)
                                                                                                                                               tATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACgcgtgcg       >  1:5079401/1‑150 (MQ=255)
                                                                                                                                                aTGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACgcgtgcgt      <  2:955205/150‑1 (MQ=255)
                                                                                                                                                   cTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACGCGTGCGTGCa   >  1:335853/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATTGGTAGCTGGCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACGCGTGCGTGCAA  >  CP009273/1281368‑1281665

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