Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,488,966 G→C 90.5% A533P (GCC→CCC)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,488,9660GC90.5% 722.2 / 53.2 203A533P (GCC→CCC) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base C (105/78);  total (116/87)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

CCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTG  >  CP009273/3488821‑3489112
                                                                                                                                                 |                                                                                                                                                  
ccTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcctc                                                                                                                                                <  1:1305379/150‑1 (MQ=255)
ccTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcctc                                                                                                                                                >  1:1770351/1‑150 (MQ=255)
ccTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcctc                                                                                                                                                <  1:4194668/150‑1 (MQ=255)
  tgtgCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcctcct                                                                                                                                              <  1:3347044/150‑1 (MQ=255)
  tgtgCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcctcct                                                                                                                                              >  1:4151268/1‑150 (MQ=255)
   gtgCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCctg                                                                                                                                             <  1:97222/150‑1 (MQ=255)
     gCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCctgct                                                                                                                                           <  2:1401180/150‑1 (MQ=255)
       gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                         >  2:4275845/1‑150 (MQ=255)
             tttGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGaa                                                                                                                                   >  2:4342302/1‑150 (MQ=255)
                  ttATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAt                                                                                                                              <  2:4337159/150‑1 (MQ=255)
                    aTTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTc                                                                                                                            <  2:2855020/150‑1 (MQ=255)
                      tCGCGTTAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTg                                                                                                                          <  2:680986/150‑1 (MQ=255)
                       cgcgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGa                                                                                                                         >  2:3395605/1‑150 (MQ=255)
                         cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                       >  1:2664309/1‑150 (MQ=255)
                         cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                       >  2:2618058/1‑150 (MQ=255)
                           tCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGcc                                                                                                                     <  1:2010850/150‑1 (MQ=255)
                           tCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGcc                                                                                                                     <  2:2985173/150‑1 (MQ=255)
                           tCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGcc                                                                                                                     <  1:675351/150‑1 (MQ=255)
                              ggAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAACGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGa                                                                                                                  <  2:754523/150‑1 (MQ=255)
                               gAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGaa                                                                                                                 >  1:4011448/1‑150 (MQ=255)
                                aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                <  1:359647/150‑1 (MQ=255)
                                  cTGTTCCATTTGATCCGCGTTGAAGGCATTAAATCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACAt                                                                                                              <  2:112078/150‑1 (MQ=255)
                                   tGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATa                                                                                                             <  1:3501442/150‑1 (MQ=255)
                                     ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                           >  1:3247612/1‑150 (MQ=255)
                                       ccATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTcc                                                                                                         <  2:423767/150‑1 (MQ=255)
                                       ccATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTcc                                                                                                         <  2:3698753/150‑1 (MQ=255)
                                        cATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCg                                                                                                        >  1:3244081/1‑150 (MQ=255)
                                        cATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCg                                                                                                        >  2:2889761/1‑150 (MQ=255)
                                         aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGACAATACTCCGc                                                                                                       >  1:197030/1‑150 (MQ=255)
                                              aTCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAg                                                                                                  >  1:4282726/1‑150 (MQ=255)
                                              aTCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAg                                                                                                  >  1:4284947/1‑150 (MQ=255)
                                                      tGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGGTGTGAAGTGTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGACCGTATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTaa                                                                                          >  2:4311963/1‑150 (MQ=255)
                                                        aaGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACg                                                                                        >  1:492644/1‑150 (MQ=255)
                                                        aaGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACg                                                                                        <  2:2550833/150‑1 (MQ=255)
                                                        aaGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACg                                                                                        >  1:1676226/1‑150 (MQ=255)
                                                            cATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACaa                                                                                    <  2:767677/150‑1 (MQ=255)
                                                              ttAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACt                                                                                  >  1:2665286/1‑150 (MQ=255)
                                                              ttAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTAGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACt                                                                                  <  1:2341741/150‑1 (MQ=255)
                                                               tAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAGTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACtt                                                                                 >  2:859741/1‑150 (MQ=255)
                                                                aaaaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTc                                                                                >  2:4532316/1‑150 (MQ=255)
                                                                 aaaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTcc                                                                               <  1:940319/150‑1 (MQ=255)
                                                                   aCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTc                                                                             >  1:3162315/1‑150 (MQ=255)
                                                                   aCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTc                                                                             <  2:4281619/150‑1 (MQ=255)
                                                                     cTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGc                                                                           >  2:1669452/1‑150 (MQ=255)
                                                                     cTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGc                                                                           >  1:4408852/1‑150 (MQ=255)
                                                                     cTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGc                                                                           <  2:1820571/150‑1 (MQ=255)
                                                                       tCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTa                                                                         <  2:1389706/150‑1 (MQ=255)
                                                                         gaagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAAc                                                                       <  2:4043250/150‑1 (MQ=255)
                                                                         gaagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTACTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAAc                                                                       >  2:4479551/1‑150 (MQ=255)
                                                                                  ctgGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGaaa                                                              >  2:2135472/1‑150 (MQ=255)
                                                                                      cGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACg                                                          >  2:2486098/1‑150 (MQ=255)
                                                                                       gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACgg                                                         >  2:1052620/1‑150 (MQ=255)
                                                                                       gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACATCCTCGCTAACATCCAGAAAGACgg                                                         >  2:2882851/1‑150 (MQ=255)
                                                                                        aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        <  2:1607003/150‑1 (MQ=255)
                                                                                        aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGc                                                        >  1:4357611/1‑150 (MQ=255)
                                                                                           cacGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCAcc                                                     <  1:1658407/150‑1 (MQ=255)
                                                                                            acGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCt                                                    <  1:3901139/150‑1 (MQ=255)
                                                                                            acGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCt                                                    >  2:2825509/1‑150 (MQ=255)
                                                                                            acGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCt                                                    >  2:3192265/1‑150 (MQ=255)
                                                                                             cGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTa                                                   <  1:1665811/150‑1 (MQ=255)
                                                                                                 aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                               >  2:2020405/1‑150 (MQ=255)
                                                                                                 aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                               >  2:2020283/1‑150 (MQ=255)
                                                                                                 aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                               <  1:4034522/150‑1 (MQ=255)
                                                                                                 aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACACg                                               >  2:2967027/1‑150 (MQ=255)
                                                                                                    ggCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTg                                            >  2:2308666/1‑150 (MQ=255)
                                                                                                       tACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATc                                         <  1:2889178/150‑1 (MQ=255)
                                                                                                       tACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCAGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATc                                         >  1:690004/1‑150 (MQ=255)
                                                                                                          ggTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCcg                                      >  2:1279944/1‑150 (MQ=255)
                                                                                                           gTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCcgc                                     >  1:2402597/1‑150 (MQ=255)
                                                                                                           gTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCcgc                                     >  1:138228/1‑150 (MQ=255)
                                                                                                           gTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCcgc                                     >  2:1683414/1‑150 (MQ=255)
                                                                                                           gTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCcgc                                     >  1:3233826/1‑150 (MQ=255)
                                                                                                             tgtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGt                                   <  1:4509851/150‑1 (MQ=255)
                                                                                                              gtgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGtt                                  <  2:3818969/150‑1 (MQ=255)
                                                                                                                    tttGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCgg                            <  1:2235635/150‑1 (MQ=255)
                                                                                                                       gTAAACCACCCGTCGGTTCGCAGCTGGCCTCCGGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCg                         <  2:4223716/150‑1 (MQ=255)
                                                                                                                       gTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCg                         <  1:339355/150‑1 (MQ=255)
                                                                                                                       gTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCg                         >  1:2467765/1‑150 (MQ=255)
                                                                                                                         aaaCCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGt                       <  1:2594280/150‑1 (MQ=255)
                                                                                                                          aaccaaccGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTg                      <  2:4392776/150‑1 (MQ=255)
                                                                                                                           accaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGa                     >  2:2174600/1‑150 (MQ=255)
                                                                                                                            ccaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGaa                    >  2:3538818/1‑150 (MQ=255)
                                                                                                                             caaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGaaa                   >  2:1827676/1‑150 (MQ=255)
                                                                                                                              aaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAAt                  >  1:1019921/1‑150 (MQ=255)
                                                                                                                               accGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATc                 >  1:1733471/1‑150 (MQ=255)
                                                                                                                                 cGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc               <  2:1853052/150‑1 (MQ=255)
                                                                                                                                 cGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc               <  2:1852937/150‑1 (MQ=255)
                                                                                                                                 cGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc               >  1:3731941/1‑150 (MQ=255)
                                                                                                                                  gTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcc              <  2:2282117/150‑1 (MQ=255)
                                                                                                                                    cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            <  2:2588384/150‑1 (MQ=255)
                                                                                                                                    cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            <  2:2964051/150‑1 (MQ=255)
                                                                                                                                    cGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc            >  2:713412/1‑150 (MQ=255)
                                                                                                                                     ggTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCg           >  2:4427085/1‑150 (MQ=255)
                                                                                                                                     ggTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCg           >  2:323559/1‑150 (MQ=255)
                                                                                                                                       ttCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGa         >  1:4071421/1‑150 (MQ=255)
                                                                                                                                           ctgctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAggg     <  1:1688102/150‑1 (MQ=255)
                                                                                                                                           ctgctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAggg     >  1:1103311/1‑150 (MQ=255)
                                                                                                                                            tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc    <  2:4519917/150‑1 (MQ=255)
                                                                                                                                              ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg  <  1:3714510/150‑1 (MQ=255)
                                                                                                                                              ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGAGTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg  <  2:2591689/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTG  >  CP009273/3488821‑3489112

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: