Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,941,341 G→A 100% L343L (CTC→CTT amiC ← N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,941,3410GA96.6% 764.2 / 12.8 237L343L (CTC→CTTamiCN‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (120/109);  total (123/114)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 237 total aligned reads displayed.

GGAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTTGT  >  CP009273/2941196‑2941482
                                                                                                                                                 |                                                                                                                                             
ggAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAg                                                                                                                                             <  1:1572018/148‑1 (MQ=255)
 gAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGc                                                                                                                                            <  1:779070/148‑1 (MQ=255)
    aGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGACTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTa                                                                                                                                         >  2:1166475/1‑148 (MQ=255)
     gTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTAt                                                                                                                                        <  2:906868/148‑1 (MQ=255)
       cGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTc                                                                                                                                      <  2:1717078/148‑1 (MQ=255)
          cGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGt                                                                                                                                   <  1:3176178/148‑1 (MQ=255)
           gTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTa                                                                                                                                  <  2:1599299/148‑1 (MQ=255)
           gTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTa                                                                                                                                  >  1:421447/1‑148 (MQ=255)
            ttttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTAc                                                                                                                                 >  1:3508081/1‑148 (MQ=255)
             tttCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTAcc                                                                                                                                <  2:3463897/148‑1 (MQ=255)
               tCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGc                                                                                                                              >  2:3412343/1‑148 (MQ=255)
                cAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCt                                                                                                                             <  1:2512251/148‑1 (MQ=255)
                   tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                          >  1:1314076/1‑148 (MQ=255)
                   tttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTa                                                                                                                          >  1:346578/1‑148 (MQ=255)
                    ttACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTAc                                                                                                                         <  2:1107259/148‑1 (MQ=255)
                        gCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCaaa                                                                                                                     <  1:797650/148‑1 (MQ=255)
                         ctctTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAAc                                                                                                                    >  1:1690565/1‑148 (MQ=255)
                                 cAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCt                                                                                                            >  1:2864978/1‑148 (MQ=255)
                                  aaCGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTg                                                                                                           >  2:1503112/1‑148 (MQ=255)
                                    cGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGACTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTc                                                                                                         >  2:3266086/1‑148 (MQ=255)
                                      ttACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCgg                                                                                                       >  1:665395/1‑148 (MQ=255)
                                      ttACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCgg                                                                                                       <  2:2140601/148‑1 (MQ=255)
                                       tACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGc                                                                                                      >  1:384579/1‑148 (MQ=255)
                                          tGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAAt                                                                                                   >  1:2249019/1‑148 (MQ=255)
                                           gATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATg                                                                                                  <  1:473888/148‑1 (MQ=255)
                                            aTAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATgg                                                                                                 >  1:1402594/1‑148 (MQ=255)
                                             tAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGt                                                                                                >  1:3597782/1‑148 (MQ=255)
                                               aaCGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCa                                                                                              >  2:366606/1‑148 (MQ=255)
                                               aaCGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTGTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTAACGCTTTACCAAACATCAGGCTGTCGGCAATGGTCa                                                                                              >  2:1337584/1‑148 (MQ=255)
                                                   ccGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGa                                                                                          <  2:1609305/148‑1 (MQ=255)
                                                      tttCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTTATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTg                                                                                       <  1:1223554/148‑1 (MQ=255)
                                                      tttCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTg                                                                                       >  2:841385/1‑148 (MQ=255)
                                                       ttCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGt                                                                                      >  1:2828574/1‑148 (MQ=255)
                                                             cAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCata                                                                                <  2:2059088/148‑1 (MQ=255)
                                                              aGAATGGAGGGAATATGTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTTATTTTTATGCTGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCatat                                                                               <  2:1366292/148‑1 (MQ=255)
                                                                aaTGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCg                                                                             >  1:2692358/1‑148 (MQ=255)
                                                                 aTGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCTCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGa                                                                            >  2:3098870/1‑148 (MQ=255)
                                                                  tGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGaa                                                                           <  1:3030275/148‑1 (MQ=255)
                                                                  tGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGaa                                                                           >  2:655048/1‑148 (MQ=255)
                                                                   ggAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAAc                                                                          <  2:418402/148‑1 (MQ=255)
                                                                     aGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACAt                                                                        <  1:1602358/148‑1 (MQ=255)
                                                                     aGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACAt                                                                        <  1:2194093/148‑1 (MQ=255)
                                                                     aGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACAt                                                                        <  1:1269785/148‑1 (MQ=255)
                                                                      gggAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATg                                                                       >  1:204439/1‑148 (MQ=255)
                                                                        gAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGgt                                                                     >  1:887643/1‑148 (MQ=255)
                                                                        gAATATCTGGTGCCTTTAGTACGGCAAACCCGACCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCCGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGgt                                                                     <  1:2704113/148‑1 (MQ=255)
                                                                        gAATATCTGGTGCCTTTAGTACGGCAAACCCGACCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGgt                                                                     <  1:3507137/148‑1 (MQ=255)
                                                                         aaTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGgtg                                                                    <  2:1846311/148‑1 (MQ=255)
                                                                         aaTATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTAAGCGACTCTACCATATCGAACATGttg                                                                    >  2:2869009/1‑148 (MQ=255)
                                                                            atCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTgg                                                                 >  2:269035/1‑148 (MQ=255)
                                                                             tCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGt                                                                <  1:2403302/148‑1 (MQ=255)
                                                                             tCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGt                                                                <  2:1649056/148‑1 (MQ=255)
                                                                              cTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTc                                                               >  1:2851265/1‑148 (MQ=255)
                                                                                  tGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACa                                                           >  2:3063130/1‑148 (MQ=255)
                                                                                   gCCTTTTGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACAt                                                          <  1:1884080/148‑1 (MQ=255)
                                                                                           tACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCa                                                  >  1:3177894/1‑148 (MQ=255)
                                                                                            aCGCCAAACCCGTCATGTACAACTTGTTTTTTTTGCAGCTTGTTGTTTTTACCAAGATTATTCAGTACCGCTTTACCAAACTACCGGCTTTGGGAAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCAc                                                 <  1:2591524/148‑1 (MQ=255)
                                                                                              ggCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCg                                               <  2:3556979/148‑1 (MQ=255)
                                                                                                cAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCt                                             <  1:1683712/148‑1 (MQ=255)
                                                                                                   aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAATCTTATTCAGTACCGCTTTACCAAACTTAAGTAAGTCGGCAATGGTCAGCGACTGTAACATATCGAACATGTTTTGGTCGACATAGCGGTCACCGCtttt                                          >  1:3545929/1‑148 (MQ=255)
                                                                                                   aCCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCtttt                                          >  1:2949968/1‑148 (MQ=255)
                                                                                                    cccGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTg                                         >  2:323150/1‑148 (MQ=255)
                                                                                                     ccGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGc                                        >  2:1604692/1‑148 (MQ=255)
                                                                                                        gCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCa                                     >  1:2080484/1‑148 (MQ=255)
                                                                                                        gCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCa                                     <  2:2301792/148‑1 (MQ=255)
                                                                                                         ccTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCAc                                    >  1:1003732/1‑148 (MQ=255)
                                                                                                         ccTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCAc                                    >  1:589033/1‑148 (MQ=255)
                                                                                                              tCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccac                               <  1:2963487/148‑1 (MQ=255)
                                                                                                              tCAACTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGccac                               >  2:1213157/1‑148 (MQ=255)
                                                                                                                  cTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAAt                           >  2:503994/1‑148 (MQ=255)
                                                                                                                  cTTGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAAt                           >  1:3310502/1‑148 (MQ=255)
                                                                                                                    tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                         >  1:1836971/1‑148 (MQ=255)
                                                                                                                    tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                         >  2:1894152/1‑148 (MQ=255)
                                                                                                                    tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                         >  1:3329093/1‑148 (MQ=255)
                                                                                                                    tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                         >  1:3425124/1‑148 (MQ=255)
                                                                                                                    tGATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCa                         <  1:1470650/148‑1 (MQ=255)
                                                                                                                     gATTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCaa                        <  1:825860/148‑1 (MQ=255)
                                                                                                                      aTTTTTATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAg                       >  2:2006643/1‑148 (MQ=255)
                                                                                                                       tttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGt                      <  1:1472412/148‑1 (MQ=255)
                                                                                                                        ttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTc                     >  2:1440284/1‑148 (MQ=255)
                                                                                                                        ttttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTc                     >  1:248787/1‑148 (MQ=255)
                                                                                                                          ttATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCg                   <  2:2193714/148‑1 (MQ=255)
                                                                                                                           tATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGa                  <  1:3266086/148‑1 (MQ=255)
                                                                                                                           tATGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGa                  <  1:1363337/148‑1 (MQ=255)
                                                                                                                            aTGCAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAg                 <  1:1166475/148‑1 (MQ=255)
                                                                                                                               cAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCg              <  2:971095/148‑1 (MQ=255)
                                                                                                                               cAGCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCg              <  2:1572916/148‑1 (MQ=255)
                                                                                                                                 gCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGtt            <  2:1374039/148‑1 (MQ=255)
                                                                                                                                 gCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGtt            <  1:3081068/148‑1 (MQ=255)
                                                                                                                                 gCTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGtt            <  1:1917877/148‑1 (MQ=255)
                                                                                                                                  cTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTc           <  1:655048/148‑1 (MQ=255)
                                                                                                                                  cTTGTTGATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTc           <  1:3255029/148‑1 (MQ=255)
                                                                                                                                    tgttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTg         <  2:1551254/148‑1 (MQ=255)
                                                                                                                                    tgttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTg         >  2:882655/1‑148 (MQ=255)
                                                                                                                                     gttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTgg        <  1:1656296/148‑1 (MQ=255)
                                                                                                                                      ttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTggg       >  2:408201/1‑148 (MQ=255)
                                                                                                                                      ttgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTggg       <  2:617091/148‑1 (MQ=255)
                                                                                                                                       tgATTTTACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGt      <  2:3597782/148‑1 (MQ=255)
                                                                                                                                          ttttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtg   >  2:1482085/1‑148 (MQ=255)
                                                                                                                                           tttACCAAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTtgt  >  2:729265/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GGAAAGTCGCCGTTTTCAGTTTACGCTCTTCCTCAACGTTACTGATAAACGCCGTTTCGACCAGAATGGAGGGAATATCTGGTGCCTTTAGTACGGCAAACCCGGCCTGTTCAACTTGATTTTTATGCAGCTTGTTGATTTTACCGAGCTTATTCAGTACCGCTTTACCAAACTTCAGGCTGTCGGCAATGGTCAGCGACTGTACCATATCGAACATGGTGTGGTCGACATAGCGGTCACCGCTTTTGCTCACGCCACCAATCAAGTCCGAGGCGTTCTGGGTTTGT  >  CP009273/2941196‑2941482

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: