Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,284,205 T→C 7.8% intergenic (+138/‑204) 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,2050TC 7.8% 841.4 / 47.7 255intergenic (+138/‑204)nrfG/gltPheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base T (119/116);  new base C (8/12);  total (127/128)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 261 total aligned reads displayed.

AGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGT  >  CP009273/4284056‑4284354
                                                                                                                                                     |                                                                                                                                                     
 gaTTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATc                                                                                                                                                      <  1:3179637/150‑1 (MQ=255)
 gaTTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATc                                                                                                                                                      >  1:1925489/1‑150 (MQ=255)
  aTTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATAAAACAAATAAAAATCTTTCCATAACAAAAGGGTATTCATTAAtac                                                                                                                                                     >  2:1657823/1‑148 (MQ=255)
    tGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTg                                                                                                                                                   <  2:2516655/150‑1 (MQ=255)
    tGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTg                                                                                                                                                   >  2:4017162/1‑150 (MQ=255)
    tGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTg                                                                                                                                                   <  1:2186508/150‑1 (MQ=255)
     ggATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGc                                                                                                                                                  <  2:2976241/150‑1 (MQ=255)
           atcatcTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATaaa                                                                                                                                            >  1:2760357/1‑150 (MQ=255)
             catcTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACt                                                                                                                                          <  1:3057519/150‑1 (MQ=255)
               tcTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTcg                                                                                                                                        >  1:2977391/1‑150 (MQ=255)
               tcTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTcg                                                                                                                                        <  1:2243738/150‑1 (MQ=255)
                   gTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGtt                                                                                                                                    >  2:3677360/1‑150 (MQ=255)
                       tgaAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGttatta                                                                                                                                >  1:1874253/1‑150 (MQ=255)
                       tgaAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGttatta                                                                                                                                <  1:3334111/150‑1 (MQ=255)
                       tgaAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGttatta                                                                                                                                >  1:34848/1‑150 (MQ=255)
                       tgaAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCAATAATCAGTATAATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGttatta                                                                                                                                >  1:1874066/1‑150 (MQ=255)
                        gaAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTtattat                                                                                                                               <  2:2576988/150‑1 (MQ=255)
                         aaGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATg                                                                                                                              >  2:3948086/1‑150 (MQ=255)
                         aaGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATg                                                                                                                              >  1:1447709/1‑150 (MQ=255)
                            aTTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCAt                                                                                                                           >  2:1340772/1‑150 (MQ=255)
                            aTTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCAt                                                                                                                           <  2:2665934/150‑1 (MQ=255)
                            aTTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCAt                                                                                                                           <  2:2029142/150‑1 (MQ=255)
                              tATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTa                                                                                                                         <  1:388885/150‑1 (MQ=255)
                                   cATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCa                                                                                                                    >  2:907825/1‑150 (MQ=255)
                                    aTTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAg                                                                                                                   >  1:2588529/1‑150 (MQ=255)
                                    aTTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAg                                                                                                                   >  1:3062549/1‑150 (MQ=255)
                                    aTTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAg                                                                                                                   >  1:3316860/1‑150 (MQ=255)
                                       gAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGa                                                                                                                <  1:2990413/150‑1 (MQ=255)
                                        aGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGaa                                                                                                               <  1:351951/150‑1 (MQ=255)
                                           tGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAg                                                                                                            <  1:1271348/150‑1 (MQ=255)
                                                 aaaaaTAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTTCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCAATAATCCTTCCATAAACACGCCTAATTATGCATTAATGCAGAGAAAAGCTCtgt                                                                                                      >  1:3638537/1‑150 (MQ=255)
                                                 aaaaaTAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCtgt                                                                                                      >  1:1219625/1‑150 (MQ=255)
                                                  aaaaTAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCtgtt                                                                                                     <  2:3106105/150‑1 (MQ=255)
                                                    aataatCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCtgttgt                                                                                                   <  2:426571/150‑1 (MQ=255)
                                                    aataatCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCtgttgt                                                                                                   <  2:3417915/150‑1 (MQ=255)
                                                      taatCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTa                                                                                                 <  1:2031225/150‑1 (MQ=255)
                                                          cGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAgg                                                                                             >  1:508109/1‑150 (MQ=255)
                                                            aaTGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGt                                                                                           >  1:1188436/1‑150 (MQ=255)
                                                             aTGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGtt                                                                                          <  1:1703402/150‑1 (MQ=255)
                                                               gAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGCAAAGCTCTGTTGTTAAAGGGTTgc                                                                                        >  2:2937080/1‑150 (MQ=255)
                                                                aaCATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTgcg                                                                                       <  2:2080148/150‑1 (MQ=255)
                                                                 aCATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTgcgc                                                                                      >  1:2138374/1‑150 (MQ=255)
                                                                     atGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACa                                                                                  >  1:1094645/1‑150 (MQ=255)
                                                                     atGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATCCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACa                                                                                  <  2:3792887/150‑1 (MQ=255)
                                                                       gCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATa                                                                                <  1:3694876/150‑1 (MQ=255)
                                                                       gCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATa                                                                                <  2:1862248/150‑1 (MQ=255)
                                                                        ccAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATAc                                                                               <  2:2359527/150‑1 (MQ=255)
                                                                        ccAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATAc                                                                               <  2:3066861/150‑1 (MQ=255)
                                                                         cAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATAcc                                                                              <  1:1257562/150‑1 (MQ=255)
                                                                         cAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATAcc                                                                              <  1:458776/150‑1 (MQ=255)
                                                                          aaaaaTAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACcg                                                                             >  1:3709765/1‑150 (MQ=255)
                                                                          aaaaaTAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACcg                                                                             <  1:728362/150‑1 (MQ=255)
                                                                          aaaaaTAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACcg                                                                             <  1:1628996/150‑1 (MQ=255)
                                                                                 atCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAAt                                                                      <  1:3228314/150‑1 (MQ=255)
                                                                                    aCTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGAt                                                                   <  1:1914821/150‑1 (MQ=255)
                                                                                    aCTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGAt                                                                   >  2:1032632/1‑150 (MQ=255)
                                                                                    aCTAATCAGTATTAGTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGAATCGAAAAGCTCTGTTTTTAAAGGGTTGCGCAACATACCGCGCAAATGAt                                                                   >  2:1565555/1‑150 (MQ=255)
                                                                                     cTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATa                                                                  >  1:3752461/1‑150 (MQ=255)
                                                                                     cTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATa                                                                  <  2:87631/150‑1 (MQ=255)
                                                                                      tAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATAc                                                                 >  1:1438469/1‑150 (MQ=255)
                                                                                       aaTCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACt                                                                <  1:2904288/150‑1 (MQ=255)
                                                                                        aTCAGTATTATTGCAGATTAAACAAATAAAATTCTTTCCATAACAAATTGTTATTAATTAACACTGCCATAAACTCGCTTTATTATGCATTAATGAAGAGAAAAGATCTGTTGTTAAAGGGTTGCGCAACATACAGCGCAAATGATAcag                                                               >  1:2794710/1‑148 (MQ=255)
                                                                                        aTCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTg                                                               <  1:3652731/150‑1 (MQ=255)
                                                                                          cAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGAt                                                             <  2:544289/150‑1 (MQ=255)
                                                                                          cAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGAt                                                             <  2:544356/150‑1 (MQ=255)
                                                                                           aGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATc                                                            <  2:14948/150‑1 (MQ=255)
                                                                                           aGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATc                                                            >  2:3902606/1‑150 (MQ=255)
                                                                                           aGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATc                                                            >  2:3902632/1‑150 (MQ=255)
                                                                                              tttattTCAGATTATACATATAAAAATTTTTCGATTACATATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGATAAGCTCTGTTGTTTAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATa                                                         <  1:711087/149‑1 (MQ=255)
                                                                                                tattGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAg                                                       >  1:3486802/1‑150 (MQ=255)
                                                                                                   tGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGt                                                    >  1:2298690/1‑150 (MQ=255)
                                                                                                   tGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGt                                                    >  1:566883/1‑150 (MQ=255)
                                                                                                         ttAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTaaaaa                                              >  1:2166502/1‑150 (MQ=255)
                                                                                                           aaaCAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTaaaaaaa                                            >  1:3999984/1‑150 (MQ=255)
                                                                                                           aaaCAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGATaaaaaaa                                            >  1:3999947/1‑150 (MQ=255)
                                                                                                           aaaCAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTaaaaaaa                                            <  2:751236/150‑1 (MQ=255)
                                                                                                            aaCAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAAt                                           <  1:2101380/150‑1 (MQ=255)
                                                                                                                    aaaaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAAc                                   <  2:2487591/150‑1 (MQ=255)
                                                                                                                     aaaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAAcc                                  >  1:98144/1‑150 (MQ=255)
                                                                                                                      aaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCa                                 >  1:1765619/1‑150 (MQ=255)
                                                                                                                      aaTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCa                                 >  1:2276648/1‑150 (MQ=255)
                                                                                                                       aTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCaa                                <  2:1203381/150‑1 (MQ=255)
                                                                                                                       aTCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCaa                                <  2:2260601/150‑1 (MQ=255)
                                                                                                                        tCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAAc                               <  2:2877171/150‑1 (MQ=255)
                                                                                                                              cATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                         <  2:1078001/150‑1 (MQ=255)
                                                                                                                               aTAACAAATGGTTATTCATTAACCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAGGCGTTAAAAAAATCTACAAACCAACGCAacac                        >  1:916847/1‑150 (MQ=255)
                                                                                                                                tAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                       <  2:1152478/150‑1 (MQ=255)
                                                                                                                                      aTGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                 >  1:3776012/1‑150 (MQ=255)
                                                                                                                                      aTGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                 >  2:3321457/1‑150 (MQ=255)
                                                                                                                                       tGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATg                >  2:3231252/1‑150 (MQ=255)
                                                                                                                                         gTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc              >  1:2929780/1‑150 (MQ=255)
                                                                                                                                         gTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc              <  2:3474274/150‑1 (MQ=255)
                                                                                                                                            aTTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg           >  1:901775/1‑150 (MQ=255)
                                                                                                                                                  tAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa     >  2:387152/1‑150 (MQ=255)
                                                                                                                                                   aaTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt    <  2:1619237/150‑1 (MQ=255)
                                                                                                                                                   aaTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt    >  2:82751/1‑150 (MQ=255)
                                                                                                                                                    aTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg   <  2:3470037/150‑1 (MQ=255)
                                                                                                                                                    aTCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg   >  1:3574510/1‑150 (MQ=255)
                                                                                                                                                    aTCCTGCCATAAACTCGCGGTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg   >  2:2768943/1‑150 (MQ=255)
                                                                                                                                                     tCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTTCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGt  >  2:1507698/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGATTGGATTAATCATCTGGTGATGAAGATTATGTCATTGAGCTGCATAAAAAATAATCGAATGAACATATGCCAAAAATAATCACTAATCAGTATTATTGCAGATTAAACAAATAAAAATCTTTCCATAACAAATGGTTATTCATTAATCCTGCCATAAACTCGCGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGT  >  CP009273/4284056‑4284354

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: