Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,488,966 G→C 89.7% 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,9660GC89.7% 1081.1 / 102.4 291A533P (GCC→CCC) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base G (19/11);  new base C (159/102);  total (178/113)
Fisher's exact test for biased strand distribution p-value = 8.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

CAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGA  >  CP009273/3488818‑3489113
                                                                                                                                                    |                                                                                                                                                   
     tgtgCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCTCTGCTGCcctcct                                                                                                                                               <  1:3680650/150‑1 (MQ=255)
     tgtgCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCcctcct                                                                                                                                               <  1:4378641/150‑1 (MQ=255)
        gCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCctgct                                                                                                                                            >  2:4740972/1‑150 (MQ=255)
          gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                          >  2:1869679/1‑150 (MQ=255)
          gAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTgg                                                                                                                                          >  2:3949144/1‑150 (MQ=255)
           aaCACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGa                                                                                                                                         >  1:1792721/1‑150 (MQ=255)
            acacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGaa                                                                                                                                        >  2:1423987/1‑150 (MQ=255)
             cacTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAAc                                                                                                                                       >  2:62305/1‑150 (MQ=255)
                tttGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGaa                                                                                                                                    >  2:1478499/1‑150 (MQ=255)
                 ttGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAAt                                                                                                                                   >  2:2524997/1‑150 (MQ=255)
                   gCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAc                                                                                                                                 >  2:3910571/1‑150 (MQ=255)
                   gCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAc                                                                                                                                 >  2:5130806/1‑150 (MQ=255)
                   gCTTAGTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATAc                                                                                                                                 <  2:2796581/150‑1 (MQ=255)
                     ttATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACAt                                                                                                                               >  2:4161048/1‑150 (MQ=255)
                       aTTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTc                                                                                                                             <  1:2971617/150‑1 (MQ=255)
                        ttCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCt                                                                                                                            >  1:86637/1‑150 (MQ=255)
                          cgcgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGa                                                                                                                          >  1:3910847/1‑150 (MQ=255)
                          cgcgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGa                                                                                                                          >  1:1636347/1‑150 (MQ=255)
                           gcgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGaa                                                                                                                         <  2:1788815/150‑1 (MQ=255)
                            cgTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAg                                                                                                                        >  2:4659014/1‑150 (MQ=255)
                             gTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGc                                                                                                                       <  2:667916/150‑1 (MQ=255)
                              tCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGcc                                                                                                                      <  1:2705763/150‑1 (MQ=255)
                                aGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCgg                                                                                                                    >  2:4324579/1‑150 (MQ=255)
                                aGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCgg                                                                                                                    >  2:2036820/1‑150 (MQ=255)
                                aGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCgg                                                                                                                    >  2:4433444/1‑150 (MQ=255)
                                  gAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCAGCTGGAACGAATACATTCTGAAGCCGGaa                                                                                                                  >  2:4354440/1‑150 (MQ=255)
                                   aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                 >  1:287242/1‑150 (MQ=255)
                                   aaCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAAc                                                                                                                 >  1:2572423/1‑150 (MQ=255)
                                    aCTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACa                                                                                                                <  1:661391/150‑1 (MQ=255)
                                      tGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATa                                                                                                              <  1:4034707/150‑1 (MQ=255)
                                        ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                            >  1:4562458/1‑150 (MQ=255)
                                        ttCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACt                                                                                                            <  1:100095/150‑1 (MQ=255)
                                          ccATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTcc                                                                                                          >  2:3412225/1‑150 (MQ=255)
                                           cATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCg                                                                                                         <  1:2008807/150‑1 (MQ=255)
                                            aTTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGc                                                                                                        >  2:4611097/1‑150 (MQ=255)
                                             tttGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCt                                                                                                       <  1:1070614/150‑1 (MQ=255)
                                              ttGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTg                                                                                                      >  2:858438/1‑150 (MQ=255)
                                                gATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCa                                                                                                    >  1:3899771/1‑150 (MQ=255)
                                                gATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCa                                                                                                    <  2:3094161/150‑1 (MQ=255)
                                                      cgTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTc                                                                                              >  1:1468628/1‑150 (MQ=255)
                                                      cgTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTc                                                                                              >  2:4609050/1‑150 (MQ=255)
                                                        ttGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTa                                                                                            >  1:2430622/1‑150 (MQ=255)
                                                         tGAAGGCAGTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGACGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTaa                                                                                           >  2:3480405/1‑150 (MQ=255)
                                                          gAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAAc                                                                                          >  1:2210911/1‑150 (MQ=255)
                                                            aGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGa                                                                                        >  1:5141404/1‑150 (MQ=255)
                                                                aTTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAAc                                                                                    >  1:3997972/1‑150 (MQ=255)
                                                                  tAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACtt                                                                                  <  2:2294452/150‑1 (MQ=255)
                                                                    anaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTcc                                                                                >  1:1220999/1‑150 (MQ=255)
                                                                     aaCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCt                                                                               >  1:4998443/1‑150 (MQ=255)
                                                                      aCCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTc                                                                              <  2:1069177/150‑1 (MQ=255)
                                                                       ccTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCg                                                                             <  1:3764297/150‑1 (MQ=255)
                                                                       ccTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCg                                                                             <  1:3558230/150‑1 (MQ=255)
                                                                           cGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTaa                                                                         >  2:4005993/1‑150 (MQ=255)
                                                                             aagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACa                                                                       <  2:2179445/150‑1 (MQ=255)
                                                                             aagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACa                                                                       >  1:3323096/1‑150 (MQ=255)
                                                                             aagaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACa                                                                       >  1:745534/1‑150 (MQ=255)
                                                                              agaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACAt                                                                      <  2:1537500/150‑1 (MQ=255)
                                                                              agaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACAt                                                                      >  2:2312440/1‑150 (MQ=255)
                                                                              agaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACAt                                                                      >  1:2477897/1‑150 (MQ=255)
                                                                              agaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACAt                                                                      <  2:4869032/150‑1 (MQ=255)
                                                                               gaaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATc                                                                     >  2:1634635/1‑150 (MQ=255)
                                                                                aaCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATcc                                                                    >  2:751277/1‑150 (MQ=255)
                                                                                 aCTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCa                                                                   >  2:4457906/1‑150 (MQ=255)
                                                                                   tgctgGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGa                                                                 >  1:3637701/1‑150 (MQ=255)
                                                                                   tgctgGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGa                                                                 <  1:2155320/150‑1 (MQ=255)
                                                                                      tgGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAg                                                              <  2:1381147/150‑1 (MQ=255)
                                                                                        gCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGAc                                                            >  2:4811880/1‑150 (MQ=255)
                                                                                        gCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGAc                                                            <  1:4725572/150‑1 (MQ=255)
                                                                                          gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACgg                                                          >  2:3290227/1‑150 (MQ=255)
                                                                                          gAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACgg                                                          >  2:4629813/1‑150 (MQ=255)
                                                                                           aaaCACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCAGCTGCCCACCTGCTGGAACGAATACAATCTGAAGCCGGAAAATACTCCGCAGCAGGATTCTAACGACAACTTCCTCGCTCACATCCAGAAAGACGGc                                                         >  1:3096291/1‑150 (MQ=255)
                                                                                              cacGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCAcc                                                      <  2:318449/150‑1 (MQ=255)
                                                                                                cGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTa                                                    <  2:3660661/150‑1 (MQ=255)
                                                                                                 ggCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTAc                                                   >  2:3725789/1‑150 (MQ=255)
                                                                                                    aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                                >  2:3084531/1‑150 (MQ=255)
                                                                                                    aaaGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCg                                                >  1:3986177/1‑150 (MQ=255)
                                                                                                     aaGGCTACGGTTGTCAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCgg                                               >  1:1582438/1‑150 (MQ=255)
                                                                                                                  tgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                  >  2:2732312/1‑150 (MQ=255)
                                                                                                                  tgAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                  >  2:1257776/1‑150 (MQ=255)
                                                                                                                  tgAAGTTAGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtc                                  <  2:932268/150‑1 (MQ=255)
                                                                                                                   gAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTtct                                 <  1:3068711/150‑1 (MQ=255)
                                                                                                                     aGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTcc                               >  1:5057841/1‑150 (MQ=255)
                                                                                                                       tttGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCgg                             >  1:3203501/1‑150 (MQ=255)
                                                                                                                       tttGTAAACCAACCGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCgg                             <  1:4933125/150‑1 (MQ=255)
                                                                                                                         tGTAACCCAACCGTCGGTTCGCTGCGGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGc                           <  1:1350973/150‑1 (MQ=255)
                                                                                                                             aaccaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTg                       <  1:4211913/150‑1 (MQ=255)
                                                                                                                               ccaaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGaa                     >  1:694877/1‑150 (MQ=255)
                                                                                                                                 aaccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAAt                   >  1:4125860/1‑150 (MQ=255)
                                                                                                                                  accGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATc                  >  2:562272/1‑150 (MQ=255)
                                                                                                                                   ccGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCa                 >  2:4943823/1‑150 (MQ=255)
                                                                                                                                    cGTCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAc                <  1:1847898/150‑1 (MQ=255)
                                                                                                                                      tCGGTTCGCTGCTGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAccc              <  1:4462475/150‑1 (MQ=255)
                                                                                                                                       cGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACGTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCAcccc             <  1:1142747/150‑1 (MQ=255)
                                                                                                                                          ttCGCTGCTGCCCTCCTGCTGGAACGAATACAGTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGa          <  2:2996414/150‑1 (MQ=255)
                                                                                                                                             gctgctGCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAgg       >  2:2613876/1‑150 (MQ=255)
                                                                                                                                               tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc     <  1:1202723/150‑1 (MQ=255)
                                                                                                                                               tgctgcCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGc     <  2:2507781/150‑1 (MQ=255)
                                                                                                                                                 ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg   >  2:5023440/1‑150 (MQ=255)
                                                                                                                                                 ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg   >  1:3067451/1‑150 (MQ=255)
                                                                                                                                                 ctgCCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTg   >  2:3557689/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACCAACCGTCGGTTCGCTGCTGGCCTCCTGCTGGAACGAATACATTCTGAAGCCGGAACATACTCCGCTGCAGGATTCTAACGACAACTTCCTCGCTAACATCCAGAAAGACGGCACCTACTCGGTGATCCCGCGTTCTCCGGGCGGTGAAATCACCCCGGAAGGGCTGA  >  CP009273/3488818‑3489113

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: