Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,653,668 G→T 8.7% intergenic (‑20/+8) iscU ← / ← iscS iron‑sulfur cluster assembly scaffold protein/cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,653,6680GT 8.7% 1212.7 / 64.4 346intergenic (‑20/+8)iscU/iscSiron‑sulfur cluster assembly scaffold protein/cysteine desulfurase (tRNA sulfurtransferase), PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (157/158);  new base T (16/15);  total (173/173)
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 348 total aligned reads displayed.

GCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTAAAACGA  >  CP009273/2653519‑2653817
                                                                                                                                                     |                                                                                                                                                     
    cTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATAc                                                                                                                                                   <  2:1398148/150‑1 (MQ=255)
    cTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATAc                                                                                                                                                   >  1:3712261/1‑150 (MQ=255)
     ttCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATaac                                                                                                                                                  >  1:685440/1‑148 (MQ=255)
     ttCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATAcc                                                                                                                                                  <  2:4533165/150‑1 (MQ=255)
     ttCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATAcc                                                                                                                                                  <  1:2589224/150‑1 (MQ=255)
       catcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGa                                                                                                                                                >  2:2160986/1‑150 (MQ=255)
       catcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGa                                                                                                                                                >  1:4921786/1‑150 (MQ=255)
       catcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGa                                                                                                                                                <  1:4275006/150‑1 (MQ=255)
       catcaCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGa                                                                                                                                                <  1:2808759/150‑1 (MQ=255)
          caCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCAGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTa                                                                                                                                             <  2:1164913/150‑1 (MQ=255)
           aCGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTaa                                                                                                                                            <  2:1442216/150‑1 (MQ=255)
            cGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAat                                                                                                                                           <  2:556125/150‑1 (MQ=255)
            cGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAat                                                                                                                                           <  1:338461/150‑1 (MQ=255)
               cGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAatgat                                                                                                                                        >  2:3923153/1‑150 (MQ=255)
                gCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAatgatg                                                                                                                                       >  2:68215/1‑150 (MQ=255)
                  cacaGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAg                                                                                                                                     <  1:69376/150‑1 (MQ=255)
                   acaGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGc                                                                                                                                    <  1:2991675/150‑1 (MQ=255)
                         cGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCTATACCGATTAATGATGAGCCCATTc                                                                                                                              <  2:3464448/150‑1 (MQ=255)
                         cGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTc                                                                                                                              <  2:481928/150‑1 (MQ=255)
                         cGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTc                                                                                                                              >  1:166724/1‑150 (MQ=255)
                         cGGTGCCCCCACCATGCCGCGGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTc                                                                                                                              <  2:3710817/150‑1 (MQ=255)
                              cccccACCATGCCGCTGCCGACGTTCTCGTCGTTGTAGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGc                                                                                                                         >  2:4680152/1‑150 (MQ=255)
                               ccccACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCt                                                                                                                        <  1:4795324/150‑1 (MQ=255)
                                       tGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGAt                                                                                                                <  1:1061071/150‑1 (MQ=255)
                                        gCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTAGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTAAGATc                                                                                                               >  1:394012/1‑150 (MQ=255)
                                           gCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCAc                                                                                                            <  2:4150360/150‑1 (MQ=255)
                                           gCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCAc                                                                                                            <  1:2054504/150‑1 (MQ=255)
                                                cGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCt                                                                                                       >  1:5166177/1‑150 (MQ=255)
                                                    gTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCtt                                                                                                   <  1:4450981/150‑1 (MQ=255)
                                                      tctcGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCTATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGt                                                                                                 >  1:4706080/1‑150 (MQ=255)
                                                        tcgtcgTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTAc                                                                                               <  2:3508084/150‑1 (MQ=255)
                                                          gtcgtTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACAt                                                                                             >  2:4526100/1‑150 (MQ=255)
                                                                gttgtCAATGGAACCCACGTTACGCGTATTCTCGTTATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCTATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCca                                                                                       <  2:3784717/150‑1 (MQ=255)
                                                                gttgtCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCca                                                                                       <  1:2912943/150‑1 (MQ=255)
                                                                   gtCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAg                                                                                    >  2:5059900/1‑150 (MQ=255)
                                                                    tCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGTTGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGc                                                                                   >  2:1303416/1‑150 (MQ=255)
                                                                       aaGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGga                                                                                <  2:4089509/150‑1 (MQ=255)
                                                                        aGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATTCCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGgag                                                                               >  1:3919234/1‑150 (MQ=255)
                                                                         ggAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCAGTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGgaga                                                                              >  2:4596526/1‑150 (MQ=255)
                                                                           aaCCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGgagaca                                                                            >  2:2643118/1‑148 (MQ=255)
                                                                              ccACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGt                                                                         >  2:4657570/1‑150 (MQ=255)
                                                                              ccACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGt                                                                         <  2:2747976/150‑1 (MQ=255)
                                                                               cACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACCGCGGAGAAAGGTc                                                                        >  1:3103157/1‑150 (MQ=255)
                                                                               cACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTc                                                                        <  2:1358478/150‑1 (MQ=255)
                                                                                 cGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCAc                                                                      <  2:3112774/150‑1 (MQ=255)
                                                                                      cgcgGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAACGGTCACGCAga                                                                 >  1:2787794/1‑150 (MQ=255)
                                                                                       gcgGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAgac                                                                >  2:2571319/1‑150 (MQ=255)
                                                                                       gcgGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAgac                                                                <  1:1336263/150‑1 (MQ=255)
                                                                                        cgGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAgacg                                                               >  2:3394419/1‑150 (MQ=255)
                                                                                        cgGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAgacg                                                               >  2:3183092/1‑150 (MQ=255)
                                                                                           aTTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGAcc                                                            <  1:3746552/150‑1 (MQ=255)
                                                                                           aTTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGAcc                                                            <  1:5156257/150‑1 (MQ=255)
                                                                                           aTTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTACCACAGCGGAGAAAGGTCAAGCAGAAGAcc                                                            >  1:1951229/1‑150 (MQ=255)
                                                                                             tctcGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGa                                                          >  1:4114092/1‑150 (MQ=255)
                                                                                               tcGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATg                                                        <  2:4947030/150‑1 (MQ=255)
                                                                                                 gTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGa                                                      <  2:1852001/150‑1 (MQ=255)
                                                                                                  tAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTACTGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGAt                                                     <  2:1663014/150‑1 (MQ=255)
                                                                                                  tAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGAt                                                     >  2:3597375/1‑150 (MQ=255)
                                                                                                   aaTGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGAtt                                                    >  2:3641998/1‑150 (MQ=255)
                                                                                                   aaTGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGAtt                                                    <  1:4347719/150‑1 (MQ=255)
                                                                                                    aTGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCTATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGAttt                                                   <  1:3961819/150‑1 (MQ=255)
                                                                                                       gTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACg                                                >  1:4592207/1‑150 (MQ=255)
                                                                                                        tCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGa                                               >  2:4426945/1‑150 (MQ=255)
                                                                                                          gATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGaac                                             >  1:1381558/1‑150 (MQ=255)
                                                                                                           aTAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGaact                                            <  1:1513487/150‑1 (MQ=255)
                                                                                                            tAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGaacta                                           >  1:2363544/1‑150 (MQ=255)
                                                                                                              aCTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGaactaac                                         <  2:1547489/150‑1 (MQ=255)
                                                                                                               cTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGaactaact                                        >  1:4302329/1‑150 (MQ=255)
                                                                                                                   ttCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGAt                                    >  1:4895454/1‑150 (MQ=255)
                                                                                                                   ttCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGAt                                    <  1:978648/150‑1 (MQ=255)
                                                                                                                   ttCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGAt                                    >  2:1988899/1‑150 (MQ=255)
                                                                                                                    tCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATg                                   >  1:159827/1‑150 (MQ=255)
                                                                                                                    tCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATg                                   >  1:400629/1‑150 (MQ=255)
                                                                                                                      gCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGgt                                 <  1:873946/150‑1 (MQ=255)
                                                                                                                       cTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGgtg                                <  1:318648/150‑1 (MQ=255)
                                                                                                                       cTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGgtg                                <  2:1059414/150‑1 (MQ=255)
                                                                                                                        tGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGgtgt                               <  1:1495811/150‑1 (MQ=255)
                                                                                                                          tAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAg                             <  2:2975290/150‑1 (MQ=255)
                                                                                                                            aGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTc                           <  1:2571319/150‑1 (MQ=255)
                                                                                                                             gCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCg                          <  1:3633076/150‑1 (MQ=255)
                                                                                                                               cATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGAt                        <  2:3958492/150‑1 (MQ=255)
                                                                                                                               cATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGAt                        <  2:1081597/150‑1 (MQ=255)
                                                                                                                                   ataAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATctct                    >  2:2176873/1‑150 (MQ=255)
                                                                                                                                   ataAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGATGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATctct                    <  1:5244194/150‑1 (MQ=255)
                                                                                                                                    taAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGAGTTACGAACTAACTCGATGGTGTAGTCGATCtctt                   >  2:4879847/1‑150 (MQ=255)
                                                                                                                                        ttCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCttcttc               >  2:2349995/1‑150 (MQ=255)
                                                                                                                                        ttCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCttcttc               <  1:2344073/150‑1 (MQ=255)
                                                                                                                                        ttCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGATCAGATCCACGCCCTGCTTGTATATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCttcttc               >  1:770609/1‑150 (MQ=255)
                                                                                                                                          ctcCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCag             >  2:4021697/1‑150 (MQ=255)
                                                                                                                                             cTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCagtag          >  2:5125695/1‑150 (MQ=255)
                                                                                                                                              tGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCagtagt         >  2:52118/1‑150 (MQ=255)
                                                                                                                                              tGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCagtagt         >  2:4251469/1‑150 (MQ=255)
                                                                                                                                               gATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAgtagta        >  1:1223440/1‑150 (MQ=255)
                                                                                                                                               gATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAgtagta        <  2:4038869/150‑1 (MQ=255)
                                                                                                                                                aTTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTaa       >  2:823511/1‑150 (MQ=255)
                                                                                                                                                  tCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTaaaa     >  2:2254447/1‑150 (MQ=255)
                                                                                                                                                   ccGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTAAAAc    <  2:2814716/150‑1 (MQ=255)
                                                                                                                                                   ccGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTAAAAc    <  2:169152/150‑1 (MQ=255)
                                                                                                                                                    cGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTAAAACg   >  1:1258170/1‑150 (MQ=255)
                                                                                                                                                     gATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTAAAACGa  >  2:4913099/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAACTTCATCACGTCGCCACAGGCCGGTGCCCCCACCATGCCGCTGCCGACGTTCTCGTCGTTGTTGTCAAAGGAACCCACGTTACGCGGATTCTCGTAATGGTCGATAACTTTTTCGCTGTAAGCCATTATAAATTCTCCTGATTCCGATACCGATTAATGATGAGCCCATTCGATGCTGTTCAGATCCACGCCCTGCTTGTACATTTCCCACAGCGGAGAAAGGTCACGCAGACGACCGATGGATTTACGAACTAACTCGATGGTGTAGTCGATCTCTTCTTCAGTAGTAAAACGA  >  CP009273/2653519‑2653817

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: