Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,184,471 A→T 40.8% I519I (ATT→ATA)  thiC ← phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,184,4710AT40.8% 158.0 / 288.6 223I519I (ATT→ATA) thiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (61/71);  new base T (41/50);  total (102/121)
Fisher's exact test for biased strand distribution p-value = 8.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

TTCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCGG  >  CP009273/4184324‑4184617
                                                                                                                                                   |                                                                                                                                                  
  cATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTAt                                                                                                                                                  >  1:1291085/1‑148 (MQ=255)
   aTCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATc                                                                                                                                                 <  1:3098619/148‑1 (MQ=255)
     cGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCtt                                                                                                                                               >  2:2442863/1‑148 (MQ=255)
     cGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCtt                                                                                                                                               >  2:318561/1‑148 (MQ=255)
       tGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGATATTCGAAGCGGGATTTCGACATGGCGTTATCGCGAATTGGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGTCAGCAATCTtat                                                                                                                                             <  2:935820/148‑1 (MQ=255)
       tGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTtat                                                                                                                                             >  1:143070/1‑148 (MQ=255)
          tAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAgg                                                                                                                                          <  1:3152111/148‑1 (MQ=255)
           aaGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGgcagca                                                                                                                                                    <  1:1159516/137‑1 (MQ=255)
           aaGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGt                                                                                                                                         >  2:346257/1‑148 (MQ=255)
            aGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTg                                                                                                                                        <  2:1432587/148‑1 (MQ=255)
                cgggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATaa                                                                                                                                    <  1:3154457/148‑1 (MQ=255)
                cgggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATaa                                                                                                                                    <  1:213531/148‑1 (MQ=255)
                 gggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAg                                                                                                                                   <  2:2228845/148‑1 (MQ=255)
                  ggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGc                                                                                                                                  >  1:831548/1‑148 (MQ=255)
                       tAAAGGGGTGGAGTGCCTGATTAATCTGGTCTTCCCAGCGAAATTCGTAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTg                                                                                                                             <  1:2672655/148‑1 (MQ=255)
                       tAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTg                                                                                                                             <  1:1657019/148‑1 (MQ=255)
                         aaCGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCt                                                                                                                           >  1:2259222/1‑148 (MQ=255)
                           cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTa                                                                                                                         >  1:3024276/1‑148 (MQ=255)
                           cGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTa                                                                                                                         <  1:3258452/148‑1 (MQ=255)
                            gggTCGAGGGCCAGATTACACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTaa                                                                                                                        <  2:527219/148‑1 (MQ=255)
                            gggTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTaa                                                                                                                        <  2:710942/148‑1 (MQ=255)
                             ggTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAAc                                                                                                                       <  2:1099757/148‑1 (MQ=255)
                               tCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAt                                                                                                                     <  2:1595742/148‑1 (MQ=255)
                                 gAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtct                                                                                                                   <  2:2787884/148‑1 (MQ=255)
                                   gggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtcttc                                                                                                                 <  1:2907612/148‑1 (MQ=255)
                                    ggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtcttct                                                                                                                <  1:1436097/148‑1 (MQ=255)
                                     gCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATcttctt                                                                                                               >  1:1355278/1‑148 (MQ=255)
                                      ccAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCttt                                                                                                              <  2:2679475/148‑1 (MQ=255)
                                       cAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTa                                                                                                             <  1:1121811/148‑1 (MQ=255)
                                        aGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAt                                                                                                            <  1:1897842/148‑1 (MQ=255)
                                         gATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAtt                                                                                                           <  1:1054752/148‑1 (MQ=255)
                                          aTTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTg                                                                                                          >  2:431090/1‑148 (MQ=255)
                                           ttAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTgg                                                                                                         >  1:2870941/1‑148 (MQ=255)
                                            tAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTggg                                                                                                        <  1:1705591/148‑1 (MQ=255)
                                            tAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTggg                                                                                                        <  1:2772282/148‑1 (MQ=255)
                                            tAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTggg                                                                                                        >  1:1899368/1‑148 (MQ=255)
                                                 tGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAc                                                                                                   >  1:2833026/1‑148 (MQ=255)
                                                  ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAcc                                                                                                  <  2:1161101/148‑1 (MQ=255)
                                                  ggTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAcc                                                                                                  >  1:404594/1‑148 (MQ=255)
                                                     cTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAg                                                                                               >  2:1192348/1‑148 (MQ=255)
                                                     cTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAg                                                                                               <  2:2671931/148‑1 (MQ=255)
                                                     cTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAg                                                                                               >  2:1780578/1‑148 (MQ=255)
                                                     cTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAg                                                                                               >  1:2560550/1‑148 (MQ=255)
                                                       tCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGAt                                                                                             >  1:1739363/1‑148 (MQ=255)
                                                       tCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGAt                                                                                             >  1:1739266/1‑148 (MQ=255)
                                                          cAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGct                                                                                          >  2:154139/1‑148 (MQ=255)
                                                           aGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctc                                                                                         >  2:1361308/1‑148 (MQ=255)
                                                           aGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctc                                                                                         <  2:2456513/148‑1 (MQ=255)
                                                            gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                        <  1:2363999/148‑1 (MQ=255)
                                                            gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                        <  1:2333908/148‑1 (MQ=255)
                                                            gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                        <  1:1900455/148‑1 (MQ=255)
                                                               aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                     <  2:258834/148‑1 (MQ=255)
                                                                aaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTg                                                                                    <  1:1198983/148‑1 (MQ=255)
                                                                   tCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCg                                                                                 >  2:1826903/1‑148 (MQ=255)
                                                                    cGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGt                                                                                <  2:2114939/148‑1 (MQ=255)
                                                                    cGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGt                                                                                >  1:2468960/1‑148 (MQ=255)
                                                                     gAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGtt                                                                               <  2:2428292/148‑1 (MQ=255)
                                                                      aaGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTa                                                                              >  1:1052325/1‑148 (MQ=255)
                                                                      aaGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGCCGTTa                                                                              <  2:850684/148‑1 (MQ=255)
                                                                       aGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTAc                                                                             <  2:2832926/148‑1 (MQ=255)
                                                                       aGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTAc                                                                             >  1:1590428/1‑148 (MQ=255)
                                                                        gCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACg                                                                            >  2:2609304/1‑148 (MQ=255)
                                                                           ggCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                         <  1:683268/148‑1 (MQ=255)
                                                                           ggCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                         <  2:2961271/148‑1 (MQ=255)
                                                                               ttCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACaga                                                                     >  2:1154894/1‑148 (MQ=255)
                                                                                       ggCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcgcg                                                             >  2:1459564/1‑148 (MQ=255)
                                                                                        gCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCgcgc                                                            >  1:1339267/1‑148 (MQ=255)
                                                                                        gCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCgcgc                                                            <  2:1055566/148‑1 (MQ=255)
                                                                                         cGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCa                                                           >  1:981956/1‑148 (MQ=255)
                                                                                           ttATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGc                                                         >  1:495950/1‑148 (MQ=255)
                                                                                                gcgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAAc                                                    <  2:2825620/148‑1 (MQ=255)
                                                                                                gcgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAAc                                                    <  2:1353655/148‑1 (MQ=255)
                                                                                                 cgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAAcc                                                   <  2:466517/148‑1 (MQ=255)
                                                                                                     tttGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTTGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGcc                                               >  1:2944140/1‑148 (MQ=255)
                                                                                                      ttGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCa                                              >  2:3300238/1‑148 (MQ=255)
                                                                                                          gcgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatca                                          <  2:2052676/148‑1 (MQ=255)
                                                                                                          gcgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatca                                          <  2:2052578/148‑1 (MQ=255)
                                                                                                           cgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatcat                                         >  1:2956411/1‑148 (MQ=255)
                                                                                                           cgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatcat                                         >  1:3294827/1‑148 (MQ=255)
                                                                                                               cGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAAACATcgcc                                     >  1:3091056/1‑148 (MQ=255)
                                                                                                                  aTGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCa                                  >  2:218966/1‑148 (MQ=255)
                                                                                                                         ttCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATcc                           <  2:1815372/148‑1 (MQ=255)
                                                                                                                         ttCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCGTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATcc                           >  2:1135318/1‑148 (MQ=255)
                                                                                                                              cAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGAc                      >  1:2359173/1‑148 (MQ=255)
                                                                                                                              cAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGAc                      >  2:3023412/1‑148 (MQ=255)
                                                                                                                               aGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACg                     <  1:2376766/148‑1 (MQ=255)
                                                                                                                                  tCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGa                  >  2:2424906/1‑148 (MQ=255)
                                                                                                                                    aGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAg                <  1:1428305/148‑1 (MQ=255)
                                                                                                                                        gcgTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGt            <  1:986135/148‑1 (MQ=255)
                                                                                                                                            gggCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATa        <  2:1291078/148‑1 (MQ=255)
                                                                                                                                            gggCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATa        >  2:1296040/1‑148 (MQ=255)
                                                                                                                                            gggCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATa        <  2:394459/148‑1 (MQ=255)
                                                                                                                                             ggCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAg       <  1:156407/148‑1 (MQ=255)
                                                                                                                                              gcagcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGc      <  2:161725/148‑1 (MQ=255)
                                                                                                                                               cagctATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGcc     >  2:477113/1‑148 (MQ=255)
                                                                                                                                                 gcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCg   >  2:322732/1‑148 (MQ=255)
                                                                                                                                                 gcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCg   >  1:2031514/1‑148 (MQ=255)
                                                                                                                                                 gcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCg   >  1:364641/1‑148 (MQ=255)
                                                                                                                                                 gcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCg   >  2:2157235/1‑148 (MQ=255)
                                                                                                                                                  caATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCgg  <  1:1704558/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCCCGG  >  CP009273/4184324‑4184617

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: ‑