Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,460 C→A 12.5% intergenic (+499/+115) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,4600CA12.5% 341.9 / 9.3 152intergenic (+499/+115)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (61/61);  new base A (19/6);  total (84/68)
Fisher's exact test for biased strand distribution p-value = 2.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

AGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGT  >  CP009273/3263313‑3263608
                                                                                                                                                   |                                                                                                                                                    
aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                                                                                     >  1:3079486/1‑149 (MQ=38)
  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                                                                                   <  1:2820159/149‑1 (MQ=38)
  tAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATCAg                                                                                                                                                   >  1:2749121/1‑149 (MQ=17)
   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAgg                                                                                                                                                  <  1:2063418/149‑1 (MQ=37)
    caaaCCCGCCGTATCGGTTTATTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                                                                                 <  2:858639/148‑1 (MQ=11)
    aaaaCCCGCCGAAGCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATTGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACAAGGt                                                                                                                                                 <  2:331589/149‑1 (MQ=11)
      aaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTCTCGGGGAAAGTCATTCATAGTTCAGCCCCAACAAAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAAGAAGGTGa                                                                                                                                               >  2:616210/1‑149 (MQ=11)
        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCCTTCATCGTTCTGCGCCTCCCAAAGCCAAAACCCGCCCAAGCGGGGTTTTACGTACACCAGGTGaaa                                                                                                                                             >  2:3611319/1‑149 (MQ=11)
        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                                                                             <  2:3484216/149‑1 (MQ=37)
        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGaaa                                                                                                                                             >  2:2684051/1‑149 (MQ=37)
        cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCAACAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGaaa                                                                                                                                             >  2:1159515/1‑149 (MQ=11)
           gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTACAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAAACAGGTGAAACTg                                                                                                                                          >  2:1100576/1‑149 (MQ=255)
            ccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGa                                                                                                                                         <  2:2798622/149‑1 (MQ=17)
             cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                                                                        <  1:1722012/149‑1 (MQ=37)
               aaGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCAATAATAATTTTTCCCCACCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTTCGGAAAACAAGGGGAACAGACCg                                                                                                                                      >  1:2590/1‑149 (MQ=11)
                atcGGGTTATTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                                                                                                                     <  2:2496943/147‑1 (MQ=255)
                aGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                                                                                                                     <  1:3660660/149‑1 (MQ=37)
                 gCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                                                    >  1:2142414/1‑149 (MQ=255)
                 gCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATCTTTATTCCCCACCAAAATCAAAAAACCCCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGTCCgct                                                                                                                                    >  1:907899/1‑147 (MQ=11)
                  cGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATa                                                                                                                                   <  2:228927/149‑1 (MQ=37)
                   gggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATaa                                                                                                                                  >  1:2474590/1‑149 (MQ=21)
                   gggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAAAGTCATTCATCTTTGTGCCCCTACCAAATCAAAATAACGCCGAAGCGGGTTATTGAGTAAAACAGTTGAATCTGACCGTTaa                                                                                                                                  >  1:1947390/1‑149 (MQ=11)
                    ggTTTTTACGTAAAACAGGTGATACTGACCGATAAGCCGCTTTCTTTAGGGTATAGTGTCGAGGCCAGTCAATCAAGGTTGAGCCCCGCCAAAAGCACAAACCCGCCGAAGCCGGGTGTTACGTAAAGCACGTGAAACTGACCCATacg                                                                                                                                 >  2:1195718/1‑147 (MQ=11)
                     gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                                                                                <  1:1715904/149‑1 (MQ=37)
                     gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                                                                                <  2:428484/149‑1 (MQ=37)
                      tttttACGTAAATCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                                                                                               <  2:3368699/149‑1 (MQ=21)
                      ttattACGTAAAACAGGTTAAACTGACCGAGAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGcc                                                                                                                               <  2:2745714/146‑1 (MQ=255)
                       ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAAACCGCCGAAGCGGGTTTTTAAGTAAAACAGGTTAAACTTACAGATAATCCg                                                                                                                              >  1:3798797/1‑149 (MQ=11)
                       ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACGTACCGATAAGCCg                                                                                                                              >  1:3073837/1‑149 (MQ=11)
                            cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                                                                         <  1:1743823/149‑1 (MQ=40)
                             gTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                                                                        <  1:4188480/149‑1 (MQ=21)
                             gTACAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                                                                        <  2:3589797/149‑1 (MQ=17)
                             gTAAAACCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                                                                        <  2:3147521/149‑1 (MQ=255)
                                  aCAGTTGAAACTGACCGATAAGCCGCTTTCTTATGGTTATAGTGTCGTTGAAAGTCATACATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAATCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                                   <  1:1033766/149‑1 (MQ=11)
                                     ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCCGTCAATCATATTTATGCCCCGCCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAAACAGGGGAAACTGACCCATAATCCGCGTTCTTTTgggga                                                                                                                >  1:434869/1‑147 (MQ=11)
                                     ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCAATAATATTTATGCCCCACCCAAATCAAAAAAACCCCCAAGCGGGTTTTTAAGTAAAACAGGTGTAACAGACCGATAAAACGGGTTTTTTTggggg                                                                                                                >  1:771160/1‑147 (MQ=2)
                                     ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCAGTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTTCGTAAAACAGGTGTAACAGACCGATAATCCGGGTTCTTTTggggt                                                                                                                >  1:434856/1‑147 (MQ=11)
                                     ggTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAc                                                                                                                                        <  1:1721934/125‑1 (MQ=17)
                                     ggTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                                                                <  1:4234187/149‑1 (MQ=17)
                                       tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                                                                              <  2:1454431/149‑1 (MQ=35)
                                        taaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                                                             <  1:3273836/148‑1 (MQ=21)
                                        taaaCTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAg                                                                                                             <  2:3135598/148‑1 (MQ=21)
                                        gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCAATCATATTTGTGCCCCGCCAAAATCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAAACAGGTGAAAAAGACCGATAAGCCGCGTTATTTTTGGTTTAg                                                                                                             >  1:2171846/1‑149 (MQ=11)
                                        gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGGTCGATAAGCCGCTTTCTTTGGGGGATAg                                                                                                             >  2:2312235/1‑149 (MQ=11)
                                           aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                                                                          <  2:2923115/149‑1 (MQ=37)
                                              gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcgt                                                                                                       <  1:3004407/149‑1 (MQ=37)
                                                    tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAg                                                                                                 <  1:3664047/149‑1 (MQ=37)
                                                    tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAg                                                                                                 <  1:1023856/149‑1 (MQ=37)
                                                      aGCCGCATTCTTGTTCGTTTTGTGTCGTGGAAAGTCTCTACTATTTCGGCCCCTCAAAAAGCTTAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                                                               <  2:2332067/149‑1 (MQ=11)
                                                       gCCGCGTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCCATAAGCCGGTGTCTGTGGGGGAGAGAGACGAGGACAGtcc                                                                                              >  2:3314206/1‑148 (MQ=11)
                                                           ctttctttTGGGTAGAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                                                                          <  2:3658728/149‑1 (MQ=37)
                                                             ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAATCCGCTTTCTTTTGGGTATAGCGTCGTGGCAAGtcattaaac                                                                                        >  1:3570690/1‑145 (MQ=11)
                                                             ttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCGTTTGGGTATAGCGTCGTGTACAGtccataatc                                                                                        >  2:1829816/1‑142 (MQ=11)
                                                             ttctttTGGGTATAGTGTCGTGGACAGTCAATCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCATTCTTTTGGGTATATCGTCGGGGACAGTAATTAATc                                                                                        >  2:1829815/1‑149 (MQ=11)
                                                               ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                                                                                      <  1:40371/149‑1 (MQ=37)
                                                               ctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                                                                                      <  1:40331/149‑1 (MQ=37)
                                                                 tgtGGGGATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCACAAGAAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTc                                                                                    <  2:2859138/147‑1 (MQ=255)
                                                                  ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCt                                                                                   <  1:3611319/149‑1 (MQ=37)
                                                                    gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGc                                                                                 <  2:475601/149‑1 (MQ=37)
                                                                              tCGTTGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                       <  2:1949561/149‑1 (MQ=17)
                                                                              tCGTGGCCAGTCATTCATCTTTATGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                       <  2:2529169/149‑1 (MQ=255)
                                                                              tCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                                                       <  2:1950191/149‑1 (MQ=37)
                                                                                          attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                           >  2:4019616/1‑149 (MQ=37)
                                                                                              atCTTTCTGCCCCACAAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                       <  2:2610364/149‑1 (MQ=255)
                                                                                                cactctcCCCCTCCAACAGCTCAAACCCGCCGCAGCGGGTTATTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGc                                                     <  2:3370931/146‑1 (MQ=11)
                                                                                                 tttCTGACCCGCAACAATCAAAAACCCGCCGACGCGCGTTTTTTCTTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAggg                                                    <  2:611362/149‑3 (MQ=11)
                                                                                                        cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                             <  2:3079486/149‑1 (MQ=37)
                                                                                                        cccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtgtgt                                             >  1:2880872/1‑145 (MQ=37)
                                                                                                         ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTATGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                            >  2:2690703/1‑149 (MQ=37)
                                                                                                          cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                           >  1:687343/1‑149 (MQ=37)
                                                                                                             cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGTa                                        >  2:2014994/1‑149 (MQ=255)
                                                                                                               aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACCGTCATTAAACTTTATGCCCCGCCCAAAGCAAAAACCCGCCGAAGCGGGTGTTTACGGaaa                                      >  2:739780/1‑149 (MQ=11)
                                                                                                                 agtaaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATc                                    <  1:1878795/146‑1 (MQ=18)
                                                                                                                     aaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTACAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATCAGGt                                >  2:2229537/1‑149 (MQ=255)
                                                                                                                        aCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCGGTTGGGTATAGCGTCGTGGACAGTCATTCAACTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTGTTACGTAAATCAGGTGaa                             >  1:3842193/1‑149 (MQ=255)
                                                                                                                         cccGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGaaa                            <  1:1099781/149‑1 (MQ=21)
                                                                                                                          ccgccgTAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAAc                           <  2:3044162/149‑1 (MQ=12)
                                                                                                                             ccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACAGa                        >  2:2311768/1‑149 (MQ=37)
                                                                                                                              cgAATCGGTTTTTTACGTACACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAc                       <  2:3613759/149‑1 (MQ=255)
                                                                                                                                aaGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCg                     <  1:2230322/149‑1 (MQ=37)
                                                                                                                                 atcGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                    <  1:1871597/147‑1 (MQ=255)
                                                                                                                                 aGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                    >  2:3061881/1‑149 (MQ=37)
                                                                                                                                 aGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGAAAAAAACAGCCGAAGCGGGTTTTTACGTAAATAAGGTGAAACTGACCGa                    >  2:144555/1‑149 (MQ=11)
                                                                                                                                  gCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAATCAGGTGAAACTGACCGAt                   >  1:3432400/1‑149 (MQ=37)
                                                                                                                                   cGGGTTTTTACGTATAACAGCTGACTCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATa                  <  1:2994634/149‑1 (MQ=11)
                                                                                                                                    gggTTTTTTCGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                                           <  1:568564/107‑1 (MQ=25)
                                                                                                                                    gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATaa                 >  2:1641369/1‑149 (MQ=37)
                                                                                                                                    gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAAATGACCGATaa                 >  2:3815796/1‑149 (MQ=37)
                                                                                                                                       tttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTGTGGGTATAGCGTCGGGGCCAGGCATTCATCGTTCTGCCCCGCCAAAAGCAAAAACCCGCCGCAGCGGGGGTTTACGTAAATCAGGTGAAACTGACCGGTAAGcc              >  1:1950231/1‑149 (MQ=11)
                                                                                                                                        ttttACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATAAGGTGAAACTGACCGATAATCCg             >  1:1375080/1‑149 (MQ=255)
                                                                                                                                         tttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg            <  1:1170614/149‑1 (MQ=17)
                                                                                                                                              gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTCTAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAACAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt       <  2:2142414/149‑1 (MQ=255)
                                                                                                                                              gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATGTTTCTGCCCCTCCAAATTCAAAAACCCGCCGGAGCGGGTTTTTTCGTAAATCAGGTGAAACTGACCGATTAGCGGGGttgt       >  2:3587600/1‑147 (MQ=14)
                                                                                                                                              gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt       >  1:1744261/1‑149 (MQ=255)
                                                                                                                                              gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCt       >  1:1745599/1‑149 (MQ=255)
                                                                                                                                               tAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGGTTTACGTAAATCAGGTGAAACAGACCGATCAGCCGGGTTCTg      >  2:3523762/1‑149 (MQ=14)
                                                                                                                                                 aaCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTc    >  1:3513856/1‑149 (MQ=255)
                                                                                                                                                  aCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCg   >  2:1676942/1‑149 (MQ=255)
                                                                                                                                                  aCCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTAATAATTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAAAAAGGTGAAACTGACCGATAAGCCGGGTTCTGTCg   >  2:728779/1‑149 (MQ=14)
                                                                                                                                                   ccAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGt  >  1:4226996/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGT  >  CP009273/3263313‑3263608

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