Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,110,751 C→A 93.3% R574L (CGT→CTT)  yghJ ← DUF4092 family putative lipoprotein peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,110,7510CA93.3% 816.7 / 34.4 284R574L (CGT→CTT) yghJDUF4092 family putative lipoprotein peptidase
Reads supporting (aligned to +/- strand):  ref base C (6/15);  new base A (134/127);  total (140/144)
Fisher's exact test for biased strand distribution p-value = 6.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

CTGGCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAACGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTAATGTACG  >  CP009273/3110604‑3110898
                                                                                                                                                   |                                                                                                                                                    
cTGGCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGc                                                                                                                                                    >  2:3619234/1‑150 (MQ=255)
 tGGCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCa                                                                                                                                                   >  2:918489/1‑150 (MQ=255)
 tGGCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCa                                                                                                                                                   >  2:1423124/1‑150 (MQ=255)
  ggCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCaa                                                                                                                                                  >  1:1041905/1‑150 (MQ=255)
  ggCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCaa                                                                                                                                                  <  1:4730281/150‑1 (MQ=255)
   gCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGgcaaagcaaa                                                                                                                                                 <  1:3148137/150‑1 (MQ=255)
           gTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCtgtt                                                                                                                                        <  2:3429338/150‑1 (MQ=255)
           gTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCtgtt                                                                                                                                        <  2:4165324/150‑1 (MQ=255)
            tCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCtgttg                                                                                                                                       >  1:2354650/1‑150 (MQ=255)
              cGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTa                                                                                                                                     >  1:4548470/1‑150 (MQ=255)
               ggCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAg                                                                                                                                    <  1:439075/150‑1 (MQ=255)
               ggCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAg                                                                                                                                    >  2:1064110/1‑150 (MQ=255)
               ggCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACAGCCCATATTTATTAACCCCCCCCTTCCAGCTGTAACCCATCGGGCAAAGCAAAATGCTGTTGTat                                                                                                                                    >  1:3494466/1‑149 (MQ=255)
               ggCGTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAACATGCTGTTGTAg                                                                                                                                    <  2:2509420/150‑1 (MQ=255)
                   tCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGtgtg                                                                                                                                >  2:1034447/1‑150 (MQ=255)
                   tCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGtgtg                                                                                                                                >  1:3486508/1‑150 (MQ=255)
                    ccATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTgg                                                                                                                               <  1:1698913/150‑1 (MQ=255)
                        ttATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTa                                                                                                                           >  2:3859942/1‑150 (MQ=255)
                          aTCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTAcc                                                                                                                         <  1:1956412/150‑1 (MQ=255)
                           tcgtcgGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTTCCCATATTTATTAACGCCCCCGTTCCAGCTGTAAACGTTCGGGCAAAGCACAATGCTGTTGTAGTGTGGGTTACCg                                                                                                                        >  1:265893/1‑150 (MQ=255)
                              tcgGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATa                                                                                                                     >  1:4176574/1‑150 (MQ=255)
                               cgGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATaa                                                                                                                    >  1:1324629/1‑150 (MQ=255)
                                    aGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGGAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATc                                                                                                               <  2:2595730/150‑1 (MQ=255)
                                     gATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCa                                                                                                              <  1:766244/150‑1 (MQ=255)
                                      aTAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAg                                                                                                             >  1:2535413/1‑150 (MQ=255)
                                       tAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGt                                                                                                            <  2:4466663/150‑1 (MQ=255)
                                         gcgcAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGttt                                                                                                          <  2:2089656/150‑1 (MQ=255)
                                         gcgcAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGttt                                                                                                          <  2:4429563/150‑1 (MQ=255)
                                          cgcAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAAAGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTg                                                                                                         >  1:4417985/1‑150 (MQ=255)
                                             agcaCGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGccc                                                                                                      <  1:3048675/150‑1 (MQ=255)
                                                 cGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCAc                                                                                                  >  1:4698567/1‑150 (MQ=255)
                                                  gTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCAcc                                                                                                 <  1:1423080/150‑1 (MQ=255)
                                                  gTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCAcc                                                                                                 <  1:1423236/150‑1 (MQ=255)
                                                     ctcCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGAc                                                                                              <  1:309163/150‑1 (MQ=255)
                                                        cATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTg                                                                                           >  2:2054631/1‑150 (MQ=255)
                                                         aTGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGc                                                                                          >  2:1751059/1‑150 (MQ=255)
                                                           gAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGGGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGccc                                                                                        >  2:1137767/1‑150 (MQ=255)
                                                            aaGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGcccc                                                                                       >  2:2287982/1‑150 (MQ=255)
                                                            aaGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGcccc                                                                                       >  2:715590/1‑150 (MQ=255)
                                                             aGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCa                                                                                      <  1:1146361/150‑1 (MQ=255)
                                                              gTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAg                                                                                     >  1:3177963/1‑150 (MQ=255)
                                                                tcttcATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAACGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCg                                                                                   >  2:8271/1‑150 (MQ=255)
                                                                     gtGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATa                                                                              <  2:391338/149‑1 (MQ=255)
                                                                           tCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGc                                                                        <  2:3103201/150‑1 (MQ=255)
                                                                           tCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGc                                                                        <  2:3102462/150‑1 (MQ=255)
                                                                             cGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAg                                                                      >  1:754851/1‑150 (MQ=255)
                                                                               ggTCGGTGTTGAGCGTACACTGCCCATTTTTATTAACGCCCCCGTTCCTGCTGTAACGGTTGGGGCAAATCAAATTGCTGTTGTTGTGTGGGTTACCTATAACCTTAATTTTGCCCTCACCGAGTTGCCCCAGCGAAATAAACGGCAGGt                                                                    <  1:1527290/150‑1 (MQ=255)
                                                                                gTCGGTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCTTTCGGGCAAAGCAAAATGCTGTTGTTGTGTGGGTTACCGATAACCATCTTTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGtt                                                                   <  2:585919/150‑1 (MQ=255)
                                                                                  cGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTaa                                                                 <  2:1559323/150‑1 (MQ=255)
                                                                                  cGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGAAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTaa                                                                 <  2:3964921/150‑1 (MQ=255)
                                                                                     tGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTACCCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAgg                                                              <  2:4140349/150‑1 (MQ=255)
                                                                                     tGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAgg                                                              >  1:599786/1‑150 (MQ=255)
                                                                                     ggTTGAGCGTACACGGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAATGGCTGTTGTTGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAgg                                                              <  2:200787/149‑1 (MQ=255)
                                                                                        tGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTgg                                                           >  1:4440749/1‑150 (MQ=255)
                                                                                        tGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTgg                                                           >  1:3357154/1‑150 (MQ=255)
                                                                                            cGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGt                                                       <  2:3342502/150‑1 (MQ=255)
                                                                                            cGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGt                                                       >  2:2400797/1‑150 (MQ=255)
                                                                                                cacTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATcg                                                   <  2:2787363/150‑1 (MQ=255)
                                                                                                   tGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGAAAAATGCTGTTGTAGTGTGGGTTACCTATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATcgcgc                                                <  1:2749776/150‑1 (MQ=255)
                                                                                                    gCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATcgcgcg                                               >  2:562330/1‑150 (MQ=255)
                                                                                                     cccATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGt                                              >  1:4281946/1‑150 (MQ=255)
                                                                                                      ccATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTa                                             >  1:406275/1‑150 (MQ=255)
                                                                                                      caaTCTTTATTACCGCCCCCGTTCCAGCTGTCCCCGTGGGGGAAATGAAATTTGTTGTTGTTGTGTGGGTTACCTATACCATTCTGTTGGCCCACACCGAGTTCCCCCAGCGAAAAACACGGCAGGTTAAAGGTGGCGGTATCGCGCGTa                                             <  1:905836/148‑1 (MQ=255)
                                                                                                        aTCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAAc                                           >  1:1018802/1‑150 (MQ=255)
                                                                                                        aTCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAAc                                           >  2:3564556/1‑150 (MQ=255)
                                                                                                        aTCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTGTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAAc                                           <  2:1922470/150‑1 (MQ=255)
                                                                                                         tCTTTATTCACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGAAAACTGCTTTTGTAGTGTGGGTTACCGATAACAATCAGTTTGCCCCCCCCGACTTGCCCCCGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACg                                          <  1:4682669/150‑1 (MQ=255)
                                                                                                           tttATTACCGCCCCCGTTCCAGCGGTACCGGTTCGGGAAACCAAAATGGCTGTTGTTGTGTGGGTTACCTATAACCAACTGTTGGCCCACACCGACTTGCCCCAGCTAAATACACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGtt                                        <  1:2069848/150‑1 (MQ=255)
                                                                                                               ttaACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGAAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTcc                                    <  2:4614168/150‑1 (MQ=255)
                                                                                                                 aaCGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCgg                                  <  2:3446439/150‑1 (MQ=255)
                                                                                                                   cGCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCt                                >  2:4014695/1‑150 (MQ=255)
                                                                                                                    gCCCCCGTTCCAGCTGTAACCGTTCGGGCAACGCAAAATGCTGTTGTTGTGTGGGTTACCTATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTc                               <  1:1170656/150‑1 (MQ=255)
                                                                                                                    gCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTc                               <  2:2719288/150‑1 (MQ=255)
                                                                                                                    gCCCCCGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTc                               >  2:284896/1‑150 (MQ=255)
                                                                                                                     cccccGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTca                              >  2:1852676/1‑150 (MQ=255)
                                                                                                                     cccccGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTca                              >  1:3617102/1‑150 (MQ=255)
                                                                                                                      ccccGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTcaa                             >  1:1906960/1‑150 (MQ=255)
                                                                                                                        ccGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTcaaca                           <  2:2659984/150‑1 (MQ=255)
                                                                                                                        ccGTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTcaaca                           >  1:270717/1‑150 (MQ=255)
                                                                                                                          gTTCCAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAag                         >  2:4318553/1‑150 (MQ=255)
                                                                                                                          gTTCCAGCTGTAACCGTTCGGGAAAAGAAAATTGCTGTTGTTGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAag                         <  1:2763643/150‑1 (MQ=255)
                                                                                                                             ccAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAgaga                      >  1:710782/1‑150 (MQ=255)
                                                                                                                              ctggtGTACCCGTTGGGGAAAAGAAAATTGGTTTTGTTGTGTGGTTTACCGATAACAATCTTTTTGCCCTCACCGAGTTGCCCCTGCGAAATACACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGaa                     <  1:2454335/146‑1 (MQ=255)
                                                                                                                              ctggtGTAACGGTTGGGGAAAAGCAAATTGGTGTTGTTGTGTGGGTTACGGATAACCATCTGTTTGCCCTCACGGAGTTGCCCCTGCTAAATACACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGaa                     <  1:2601053/146‑1 (MQ=255)
                                                                                                                              cAGCTGTAACCGTTCGGGCAAAGCAAATTGCTGTTGTTGTGTGGGTTACCGATAACCATCTGTTTGCCCTCACCGACTTGCCCCTGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGaa                     <  1:3958441/150‑1 (MQ=255)
                                                                                                                              cAGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGaa                     >  1:1401033/1‑150 (MQ=255)
                                                                                                                              cAGCTGTAACCGTGCGGGCAAAGCAAATTGCTGTTGTAGTGTGGGTTACCTATAACCATCAGTGTGCCCTCACCGACTTGCCCCCGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGaa                     <  2:1185476/150‑1 (MQ=255)
                                                                                                                               aGCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAg                    >  2:2752208/1‑150 (MQ=255)
                                                                                                                                gCTGTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAgg                   >  2:1337700/1‑150 (MQ=255)
                                                                                                                                  tGTAACCGTTGGGGAAACGCAAAATGCTGTTGTTGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATACACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGgcg                 <  1:145382/150‑1 (MQ=255)
                                                                                                                                  tGTAACCGTTCGGGCAACGCAAAATGCTTTTGTAGTGTGGGTTACCGATAACCATCTGTTTGCCCTCACCGATTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGgcg                 <  1:3990289/150‑1 (MQ=255)
                                                                                                                                   gTAACCGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGgcgc                >  1:1253683/1‑150 (MQ=255)
                                                                                                                                       ccGTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTcc            >  1:2203535/1‑150 (MQ=255)
                                                                                                                                        cGTTCGGGCAACGCAAAATGCTGTTGTAGGGTGGGTTACCGATAACCATCGGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCg           <  2:2449500/150‑1 (MQ=255)
                                                                                                                                         gTTCGGGCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGt          <  2:4779670/150‑1 (MQ=255)
                                                                                                                                          ttCGGGAAACGAAAAATGCTGTTGTAGTGTGGTTTACCTATAACAAACTTTTGGCCCTCCCCGACTTGCCCCAGCTAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTa         <  1:1003058/150‑1 (MQ=255)
                                                                                                                                             gggCAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTAATg      >  2:1472580/1‑150 (MQ=255)
                                                                                                                                                cAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTAATGTAc   >  2:4607562/1‑150 (MQ=255)
                                                                                                                                                cAAAGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTAATGTAc   >  2:2907607/1‑150 (MQ=255)
                                                                                                                                                 aaCGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTAATGTACg  >  1:978849/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGGCTTTCGCGTCCGGCTTCCATTTATCGTCGGACAGATAGCGCAGCACGTTCTCCATGAAGTTCTTCATGTCATCCGGGTCGCTGTTGAGCGTACACTGCCCATCTTTATTAACGCCCCCGTTCCAGCTGTAACCGTTCGGGCAACGCAAAATGCTGTTGTAGTGTGGGTTACCGATAACCATCAGTTTGCCCTCACCGACTTGCCCCAGCGAAATAAACGGCAGGTTAAAGGTGGCGGTATCGCGCGTAACGTTTTCCGGCTCAACAAGAGAAGGCGCTTCCGTAATGTACG  >  CP009273/3110604‑3110898

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