Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,284,083 G→A 5.5% intergenic (+16/‑326) nrfG → / → gltP heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,284,0830GA 5.5% 1103.4 / 31.8 308intergenic (+16/‑326)nrfG/gltPheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base G (138/152);  new base A (8/9);  total (147/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

ATGCAGGCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCA  >  CP009273/4283936‑4284232
                                                                                                                                                   |                                                                                                                                                     
 tGCAGGCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGAtt                                                                                                                                                    <  1:3708149/150‑1 (MQ=255)
   cAGGCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTAt                                                                                                                                                  <  1:3439953/150‑1 (MQ=255)
    aGGCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATg                                                                                                                                                 <  2:1872037/150‑1 (MQ=255)
    aGGCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAAATTATg                                                                                                                                                 >  2:3796036/1‑150 (MQ=255)
     ggCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGt                                                                                                                                                <  2:123578/150‑1 (MQ=255)
      gCGAACTACGCGCAAGCCATAGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGAGAAATTGTTTCAGCGGAGATTGGATTAATCAAATGGAGATGACTATAATGTc                                                                                                                                               >  1:785752/1‑150 (MQ=255)
        gAATTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCAt                                                                                                                                             <  1:4417452/150‑1 (MQ=255)
        gAACTACGCTCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGATCACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATAATCTGGTGATGAAGATTATGTCAt                                                                                                                                             <  2:1078867/150‑1 (MQ=255)
        gAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCAt                                                                                                                                             >  2:1197842/1‑150 (MQ=255)
         aaCTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCAtt                                                                                                                                            <  1:4119148/150‑1 (MQ=255)
            tACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAg                                                                                                                                         >  1:2201876/1‑150 (MQ=255)
              cgcgCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCt                                                                                                                                       >  1:1017745/1‑150 (MQ=255)
                cgcAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGc                                                                                                                                     <  2:3302210/150‑1 (MQ=255)
                  cAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAAATTATGTCATTGAGCTGCAt                                                                                                                                   <  2:3190054/150‑1 (MQ=255)
                   aaGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATTGCGAAATTGTTGCATCGGAGATTGGATAAATCATATTGTGATGAAGATAATTTCATTGAGCTGCATa                                                                                                                                  >  1:1296238/1‑150 (MQ=255)
                   aaGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATa                                                                                                                                  <  1:4629883/150‑1 (MQ=255)
                   aaGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATa                                                                                                                                  >  1:1224575/1‑150 (MQ=255)
                   aaGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATa                                                                                                                                  >  1:2904822/1‑150 (MQ=255)
                      ccATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATaaaa                                                                                                                               <  1:3493502/150‑1 (MQ=255)
                       cATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAAATTATGTCATTGAGCTGCATaaaaa                                                                                                                              >  2:369810/1‑150 (MQ=255)
                         tCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAaat                                                                                                                            >  2:2488629/1‑150 (MQ=255)
                          cGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAaata                                                                                                                           <  2:4097277/150‑1 (MQ=255)
                          cGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAaata                                                                                                                           <  2:5228253/150‑1 (MQ=255)
                             aCTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATc                                                                                                                        >  1:1370842/1‑150 (MQ=255)
                             aCTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATc                                                                                                                        >  2:4319879/1‑150 (MQ=255)
                             aCTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGATGCATAAAAAATAATc                                                                                                                        >  2:2068045/1‑150 (MQ=255)
                               tATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGa                                                                                                                      >  1:922581/1‑150 (MQ=255)
                                aTGGCAAAAAGTGATGGATCTCAACTCACCGCGATTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGaa                                                                                                                     <  1:4127831/150‑1 (MQ=255)
                                 tGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAAt                                                                                                                    <  1:1096159/150‑1 (MQ=255)
                                 tGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAAt                                                                                                                    <  1:2085265/150‑1 (MQ=255)
                                  ggCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATAGAATg                                                                                                                   >  2:4123764/1‑150 (MQ=255)
                                     aaaaaGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACATCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAAc                                                                                                                <  2:2986074/150‑1 (MQ=255)
                                      aaaaGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATAATGTCATTGAGCTGCATAAAAAATAATCGAATGAACa                                                                                                               >  1:2626066/1‑150 (MQ=255)
                                       aaaGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACat                                                                                                              <  1:2255036/150‑1 (MQ=255)
                                       aaaGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACat                                                                                                              <  2:3977432/150‑1 (MQ=255)
                                        aaGTGTTGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTAGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACata                                                                                                             <  1:4854373/150‑1 (MQ=255)
                                           tGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGc                                                                                                          <  2:3834719/150‑1 (MQ=255)
                                             aTGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCa                                                                                                        <  1:3856901/150‑1 (MQ=255)
                                                gATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAAATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCaaaa                                                                                                     <  1:2692697/150‑1 (MQ=255)
                                                   ctcAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAaata                                                                                                  >  1:3294032/1‑150 (MQ=255)
                                                    tcAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAaataa                                                                                                 >  2:1502775/1‑150 (MQ=255)
                                                        cTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCAc                                                                                             <  2:47415/150‑1 (MQ=255)
                                                        cTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAAATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCAc                                                                                             >  1:3750782/1‑150 (MQ=255)
                                                         tCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACt                                                                                            >  2:494472/1‑150 (MQ=255)
                                                         tCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACt                                                                                            <  2:838985/150‑1 (MQ=255)
                                                           aCCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTaa                                                                                          <  1:1546043/150‑1 (MQ=255)
                                                            ccGCGAGTTAACCTAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTTGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTaag                                                                                         <  2:2355963/150‑2 (MQ=255)
                                                            ccGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCTGAGATTGGATAAATCATCTGGTTATGAAAATTATTTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAAt                                                                                         <  2:4215695/150‑1 (MQ=255)
                                                            ccGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAAt                                                                                         >  2:2452179/1‑150 (MQ=255)
                                                             cgcgAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATc                                                                                        >  1:4764761/1‑150 (MQ=255)
                                                             cgcgAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATc                                                                                        >  2:5182258/1‑150 (MQ=255)
                                                                gAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGt                                                                                     >  2:3612380/1‑150 (MQ=255)
                                                                    tAACCGAACACAGCTTGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGtatta                                                                                 <  1:3353596/150‑1 (MQ=255)
                                                                    tAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGtatta                                                                                 <  1:2663600/150‑1 (MQ=255)
                                                                    tAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGtatta                                                                                 >  2:484617/1‑150 (MQ=255)
                                                                        cGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGc                                                                             >  1:529960/1‑150 (MQ=255)
                                                                          aaCACCGCTGGTTGAGTAGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAg                                                                           <  2:2560492/150‑1 (MQ=255)
                                                                          aaCACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAg                                                                           >  2:5203897/1‑150 (MQ=255)
                                                                          aaCACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAg                                                                           <  1:5169959/150‑1 (MQ=255)
                                                                            cacaGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGAt                                                                         <  1:4774882/150‑1 (MQ=255)
                                                                            cacaGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGAt                                                                         <  1:3165460/150‑1 (MQ=255)
                                                                             acaGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGAtt                                                                        >  1:2210110/1‑150 (MQ=255)
                                                                               aGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTaa                                                                      >  2:2103561/1‑150 (MQ=255)
                                                                                   ggTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACaa                                                                  <  2:2201585/150‑1 (MQ=255)
                                                                                    gTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACaaa                                                                 <  1:4689228/150‑1 (MQ=255)
                                                                                    gTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATAACTAATCAGTATTATTGCAGATTAAACaaa                                                                 >  2:1170881/1‑150 (MQ=255)
                                                                                     ttGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAAt                                                                <  1:5066575/150‑1 (MQ=255)
                                                                                     ttGAGTCGAATAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAAt                                                                <  2:4436041/150‑1 (MQ=255)
                                                                                         gTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATaaaa                                                            <  2:3467326/150‑1 (MQ=255)
                                                                                             aTTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTTCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCt                                                        <  2:2628060/150‑1 (MQ=255)
                                                                                               tAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCttt                                                      <  2:3010381/150‑1 (MQ=255)
                                                                                               tAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCttt                                                      >  1:661550/1‑150 (MQ=255)
                                                                                                aaTATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAATTCTTTc                                                     >  1:3666564/1‑150 (MQ=255)
                                                                                                aaTATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTc                                                     >  1:1753993/1‑150 (MQ=255)
                                                                                                    tGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATa                                                 <  1:3110180/150‑1 (MQ=255)
                                                                                                     ggCGAAATTGTTGCATCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATaa                                                <  2:2707029/150‑1 (MQ=255)
                                                                                                     ggCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCTGATTAAACAAATAAAAATCTTTCCATaa                                                >  2:4583714/1‑150 (MQ=255)
                                                                                                      gagAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAAc                                               <  2:1742382/148‑1 (MQ=255)
                                                                                                       cGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACa                                              >  2:2208994/1‑150 (MQ=255)
                                                                                                         aaaTTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACaaa                                            >  2:558079/1‑150 (MQ=255)
                                                                                                           aTTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATg                                          >  1:355058/1‑150 (MQ=255)
                                                                                                           aTTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATg                                          <  1:2523100/150‑1 (MQ=255)
                                                                                                             tgttgCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATAAGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGt                                        >  2:1480236/1‑150 (MQ=255)
                                                                                                                tgCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTAt                                     >  1:4160047/1‑150 (MQ=255)
                                                                                                                 gCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTAtt                                    <  2:1360785/150‑1 (MQ=255)
                                                                                                                   aGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCa                                  >  1:169189/1‑150 (MQ=255)
                                                                                                                    gCGGATATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCAt                                 <  1:3984882/150‑1 (MQ=255)
                                                                                                                          aTTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATcc                           <  1:1755200/150‑1 (MQ=255)
                                                                                                                           ttGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCt                          <  1:4904365/150‑1 (MQ=255)
                                                                                                                              gATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGcc                       <  1:3886388/150‑1 (MQ=255)
                                                                                                                              gATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGcc                       <  1:5094669/150‑1 (MQ=255)
                                                                                                                                 tAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATa                    <  2:1799212/150‑1 (MQ=255)
                                                                                                                                      atcTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTc               >  1:3788109/1‑150 (MQ=255)
                                                                                                                                         tGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTcgcg            <  1:2595386/150‑1 (MQ=255)
                                                                                                                                           gTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGtt          <  1:2142118/150‑1 (MQ=255)
                                                                                                                                           gTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGtt          <  2:3140250/150‑1 (MQ=255)
                                                                                                                                            tgatgaAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGtta         >  1:1828628/1‑150 (MQ=255)
                                                                                                                                                 aaGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTGATTATg    >  1:1945018/1‑150 (MQ=255)
                                                                                                                                                   gATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCa  >  1:1886238/1‑150 (MQ=255)
                                                                                                                                                   gATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCa  >  1:3149663/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATGCAGGCGAACTACGCGCAAGCCATCGAACTATGGCAAAAAGTGATGGATCTCAACTCACCGCGAGTTAACCGAACACAGCTGGTTGAGTCGATTAATATGGCGAAATTGTTGCAGCGGAGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCA  >  CP009273/4283936‑4284232

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