Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,596,243 T→A 94.4% intergenic (‑25/+243) leuQ ← / ← rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,596,2430TA94.4% 779.2 / 27.2 249intergenic (‑25/+243)leuQ/rsmCtRNA‑Leu/16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (5/9);  new base A (148/87);  total (153/96)
Fisher's exact test for biased strand distribution p-value = 5.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 251 total aligned reads displayed.

ACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACG  >  CP009273/4596096‑4596387
                                                                                                                                                   |                                                                                                                                                
aCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGAGTGGACCCCCACGGGCGGAAGGGCCCGAAAACCTGAAGGGAGCGCGTGAAACCATTACGGCACCGTGGCCAACCATCCATTGGTAAAAAGAAtt                                                                                                                                                  >  1:743412/1‑148 (MQ=255)
 ccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                                                                                                                                 <  2:798686/148‑1 (MQ=255)
     cGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTAc                                                                                                                                             >  2:1768488/1‑148 (MQ=255)
          aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACgg                                                                                                                                        <  2:3259817/148‑1 (MQ=255)
            tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAg                                                                                                                                      >  2:1454180/1‑148 (MQ=255)
            tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAg                                                                                                                                      <  1:2979091/148‑1 (MQ=255)
                 gCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCa                                                                                                                                 >  2:1262514/1‑148 (MQ=255)
                   aaTATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCAtt                                                                                                                               <  1:2602605/148‑1 (MQ=255)
                     tatCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCt                                                                                                                             <  1:3182027/148‑1 (MQ=255)
                      atCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTa                                                                                                                            >  1:1423991/1‑148 (MQ=255)
                      atCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTa                                                                                                                            >  2:1693469/1‑148 (MQ=255)
                      atCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTa                                                                                                                            >  1:419712/1‑148 (MQ=255)
                        cAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGt                                                                                                                          >  1:3157759/1‑148 (MQ=255)
                           cGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGt                                                                                                                       >  2:773768/1‑148 (MQ=255)
                            gtggtgTTTGGTGCGCGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGtt                                                                                                                      <  2:368016/148‑1 (MQ=255)
                             tgttttttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGttt                                                                                                                     <  1:841173/145‑1 (MQ=255)
                              ggTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGtttt                                                                                                                    >  2:3539125/1‑148 (MQ=255)
                              ggTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGtttt                                                                                                                    <  1:848656/148‑1 (MQ=255)
                              ggTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGtttt                                                                                                                    >  1:601603/1‑148 (MQ=255)
                               gTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTc                                                                                                                   <  2:3009249/148‑1 (MQ=255)
                                 ttttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAg                                                                                                                 >  1:1568300/1‑148 (MQ=255)
                                 ttttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAg                                                                                                                 >  1:3131780/1‑148 (MQ=255)
                                    tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCtt                                                                                                              <  1:2693358/148‑1 (MQ=255)
                                     ggTGCGTGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCttt                                                                                                             <  1:483906/148‑1 (MQ=255)
                                      gTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCtttt                                                                                                            <  1:1976029/148‑1 (MQ=255)
                                           aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCa                                                                                                       >  1:3047560/1‑148 (MQ=255)
                                            ggggggggACTTGAAACCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCaa                                                                                                      >  2:1944427/1‑148 (MQ=255)
                                              ggggggACTTCAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                                                                                                    <  1:950775/148‑1 (MQ=255)
                                              ggggggATTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAg                                                                                                   <  1:2361858/148‑1 (MQ=255)
                                               gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAg                                                                                                   <  1:570522/148‑1 (MQ=255)
                                               gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAg                                                                                                   <  1:3372572/148‑1 (MQ=255)
                                                  ggACTTCACCCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTa                                                                                                <  1:2707941/148‑1 (MQ=255)
                                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAAt                                                                                              <  2:1531383/148‑1 (MQ=255)
                                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACAACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAAt                                                                                              >  2:672788/1‑148 (MQ=255)
                                                     cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCAtt                                                                                                                               <  1:3430089/114‑1 (MQ=255)
                                                     cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCAtt                                                                                                                               >  2:3430089/1‑114 (MQ=255)
                                                      ttGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTa                                                                                            <  1:305524/148‑1 (MQ=255)
                                                      ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTa                                                                                            >  2:3083349/1‑148 (MQ=255)
                                                      ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTa                                                                                            >  1:3328165/1‑148 (MQ=255)
                                                      ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTa                                                                                            >  2:1809556/1‑148 (MQ=255)
                                                       tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTAt                                                                                           <  1:1314272/148‑1 (MQ=255)
                                                          aCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                                                                                        >  1:2507939/1‑148 (MQ=255)
                                                             cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTGTCAGCTTTTCGTCAATAGTTAATTATCGACag                                                                                     >  2:2083947/1‑148 (MQ=255)
                                                                    ccGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTa                                                                              <  2:598989/148‑1 (MQ=255)
                                                                       tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAtt                                                                           >  2:911263/1‑148 (MQ=255)
                                                                        aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCCTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTg                                                                          <  1:1254681/148‑1 (MQ=255)
                                                                         aGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGc                                                                         <  1:2331711/148‑1 (MQ=255)
                                                                          ggACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCt                                                                        <  1:674507/148‑1 (MQ=255)
                                                                           gACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTg                                                                       >  1:3459968/1‑148 (MQ=255)
                                                                           gACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTg                                                                       >  1:1694262/1‑148 (MQ=255)
                                                                            acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGGAATTGCTgg                                                                      >  1:2627917/1‑148 (MQ=255)
                                                                               cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGaaa                                                                   >  2:2929013/1‑148 (MQ=255)
                                                                                tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGaaaa                                                                  <  2:326550/148‑1 (MQ=255)
                                                                                tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGaaaa                                                                  <  1:2634578/148‑1 (MQ=255)
                                                                                  acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCAtt                                                                                                                               >  1:1523721/1‑85 (MQ=255)
                                                                                   cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATg                                                               >  1:2210974/1‑148 (MQ=255)
                                                                                   cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATg                                                               >  1:121039/1‑148 (MQ=255)
                                                                                   cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATg                                                               >  2:2004772/1‑148 (MQ=255)
                                                                                    acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGt                                                              >  1:2979498/1‑148 (MQ=255)
                                                                                     ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTc                                                             >  1:190908/1‑148 (MQ=255)
                                                                                     ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTc                                                             >  1:2832174/1‑148 (MQ=255)
                                                                                      cTGAAGCTAGCGCGTCTACCAATTCCACCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTcc                                                            >  2:1624602/1‑148 (MQ=255)
                                                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCAt                                                          <  2:2507939/148‑1 (MQ=255)
                                                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCAt                                                          >  1:1112357/1‑148 (MQ=255)
                                                                                            cTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAgg                                                      <  2:2899286/148‑1 (MQ=255)
                                                                                            cTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAgg                                                      <  2:513476/148‑1 (MQ=255)
                                                                                              aGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGaa                                                    >  2:1554711/1‑148 (MQ=255)
                                                                                              aGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGaa                                                    >  2:3324794/1‑148 (MQ=255)
                                                                                               gcgcGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGaaa                                                   >  1:2102730/1‑148 (MQ=255)
                                                                                                 gcgTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACt                                                 <  2:2675606/148‑1 (MQ=255)
                                                                                                  cgTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTa                                                >  2:1039600/1‑148 (MQ=255)
                                                                                                     cTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCg                                             >  2:178559/1‑148 (MQ=255)
                                                                                                     cTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCg                                             >  1:2456557/1‑148 (MQ=255)
                                                                                                      tACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGt                                            <  1:639783/148‑1 (MQ=255)
                                                                                                        ccAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGc                                          >  2:332108/1‑148 (MQ=255)
                                                                                                        ccAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGc                                          <  1:3336838/148‑1 (MQ=255)
                                                                                                         cAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCa                                         <  1:1372721/148‑1 (MQ=255)
                                                                                                                 ccACCTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGt                                 <  1:54227/148‑1 (MQ=255)
                                                                                                                   aCCGTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATAAGTACACTAGCGTGCGGGTTTGGGat                               >  2:246607/1‑148 (MQ=255)
                                                                                                                     cTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTatgc                             <  2:1094579/148‑1 (MQ=255)
                                                                                                                     cTTCGCATACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGGCAATAGTTACTTATGGACAGAGGTGTCATTGCTGGCCAAATGTCCATCAGGACACTAGCGTGCCGGGTTGGTgtgc                             >  2:2732854/1‑148 (MQ=255)
                                                                                                                           aTACCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCgg                       <  1:1262514/148‑1 (MQ=255)
                                                                                                                             aCCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCgggg                     <  1:911263/148‑1 (MQ=255)
                                                                                                                             aCCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCgggg                     >  2:3115878/1‑148 (MQ=255)
                                                                                                                             aCCATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTGTGGTATGCATGCgggg                     >  1:1132750/1‑148 (MQ=255)
                                                                                                                               cATCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGc                   >  2:2423252/1‑148 (MQ=255)
                                                                                                                                 tCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCCGGTTTGGTATGCATGCGGGGGcag                 >  2:3499434/1‑148 (MQ=255)
                                                                                                                                 tCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcag                 <  2:1655035/148‑1 (MQ=255)
                                                                                                                                 tCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcag                 >  2:1691872/1‑148 (MQ=255)
                                                                                                                                 tCAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcag                 >  1:1211474/1‑148 (MQ=255)
                                                                                                                                  cAATTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcaga                <  2:2781410/148‑1 (MQ=255)
                                                                                                                                    aTTCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatg              <  2:2887347/148‑1 (MQ=255)
                                                                                                                                     ttCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc             <  1:358013/148‑1 (MQ=255)
                                                                                                                                     ttCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc             >  2:855021/1‑148 (MQ=255)
                                                                                                                                     ttCTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgc             >  2:1320135/1‑148 (MQ=255)
                                                                                                                                       cTTAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgcca           >  1:1848153/1‑148 (MQ=255)
                                                                                                                                        ttAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccag          <  2:1409348/148‑1 (MQ=255)
                                                                                                                                        ttAAAAAGAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGcagatgccag          >  1:276294/1‑148 (MQ=255)
                                                                                                                                               gAATAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGAc   >  1:2406804/1‑148 (MQ=255)
                                                                                                                                                aaTAGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACg  <  1:199116/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTGGAAAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACG  >  CP009273/4596096‑4596387

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