Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 827,319 C→T 10.2% A169A (GCC→GCT)  ycaC‑I → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947827,3190CT10.2% 482.8 / 31.3 156A169A (GCC→GCT) ycaC‑Iputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (62/78);  new base T (4/12);  total (66/90)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 157 total aligned reads displayed.

CCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCA  >  NC_002947/827171‑827468
                                                                                                                                                    |                                                                                                                                                     
 cTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTg                                                                                                                                                     >  1:5105484/1‑150 (MQ=255)
  ttCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGc                                                                                                                                                    <  2:4861692/150‑1 (MQ=255)
     ccGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGaa                                                                                                                                                 <  2:860011/150‑1 (MQ=255)
      cGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCGCCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAAc                                                                                                                                                <  2:1415530/150‑1 (MQ=255)
      cGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAAc                                                                                                                                                >  2:3028781/1‑150 (MQ=255)
       ggCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACt                                                                                                                                               <  2:3784217/150‑1 (MQ=255)
         cccTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGc                                                                                                                                             >  1:2122097/1‑150 (MQ=255)
             gTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCg                                                                                                                                         >  1:3926793/1‑150 (MQ=255)
                  cATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACt                                                                                                                                    <  1:4591656/150‑1 (MQ=255)
                     tGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGgc                                                                                                                                 <  2:1098133/150‑1 (MQ=255)
                       gaggagGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGgcgc                                                                                                                               >  2:1315342/1‑150 (MQ=255)
                        aggaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCa                                                                                                                              >  1:4249130/1‑150 (MQ=255)
                         ggaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCaa                                                                                                                             >  1:1433119/1‑150 (MQ=255)
                         ggaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCaa                                                                                                                             >  1:1433210/1‑150 (MQ=255)
                          gaggGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTCCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAAc                                                                                                                            <  1:16854/150‑1 (MQ=255)
                             ggCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGAc                                                                                                                         <  1:1267305/150‑1 (MQ=255)
                              gtttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGata                                                                                                                        <  1:2247073/148‑3 (MQ=2)
                               tttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGatat                                                                                                                       <  1:58776/149‑1 (MQ=2)
                               tttCGATGTCTTTGTCGTCACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGatat                                                                                                                       <  2:3616919/149‑1 (MQ=2)
                                ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAGCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                                                                                                      >  1:293613/1‑150 (MQ=255)
                                ttCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATc                                                                                                                      <  2:1245524/150‑1 (MQ=255)
                                   gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                                                                                                   <  1:2220625/150‑1 (MQ=255)
                                   gAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGaa                                                                                                                   <  1:4059641/150‑1 (MQ=255)
                                    aGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                  <  1:785373/150‑1 (MQ=255)
                                     ggTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAgg                                                                                                                 >  1:2102493/1‑150 (MQ=255)
                                        aTTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCt                                                                                                              <  1:191784/150‑1 (MQ=255)
                                         tttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTg                                                                                                             >  1:4616580/4‑150 (MQ=2)
                                         ttCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTg                                                                                                             <  2:3025898/150‑1 (MQ=255)
                                          ttgtcgtcACAGATGCGTCGGGAACCTTCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTgg                                                                                                            <  2:177167/148‑1 (MQ=2)
                                          tcgtcgTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTgg                                                                                                            <  1:3961710/150‑1 (MQ=255)
                                          tcgtcgTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTgg                                                                                                            <  2:2622799/150‑1 (MQ=255)
                                            gtcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGt                                                                                                          >  2:1575548/1‑150 (MQ=255)
                                             tcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTa                                                                                                         >  1:2302036/1‑150 (MQ=255)
                                             tcgtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTa                                                                                                         >  2:1945089/1‑150 (MQ=255)
                                               gtcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAcc                                                                                                       <  1:4426018/150‑1 (MQ=255)
                                                tcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTAccc                                                                                                      <  1:3227054/150‑1 (MQ=255)
                                                tcACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGAGACTGGCGCAACGACATCGAAGGGCTGGGTAccc                                                                                                      >  2:4849537/1‑150 (MQ=255)
                                                  aCCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTg                                                                                                    <  1:4939395/150‑1 (MQ=255)
                                                   ccGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGt                                                                                                   <  1:498594/150‑1 (MQ=255)
                                                    cGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGtt                                                                                                  <  2:772831/150‑1 (MQ=255)
                                                      aCGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTtct                                                                                                <  2:4020464/150‑1 (MQ=255)
                                                        gCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTcc                                                                                              <  2:4271239/150‑1 (MQ=255)
                                                           tctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAc                                                                                           <  2:3611939/150‑1 (MQ=255)
                                                            ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                          <  1:4852602/150‑1 (MQ=255)
                                                            ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                          <  1:2719304/150‑1 (MQ=255)
                                                            ctGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAAcc                                                                                          <  1:2719414/150‑1 (MQ=255)
                                                               gCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCata                                                                                       >  1:4846493/1‑150 (MQ=255)
                                                                  ccttCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCacattc                                                                                    >  2:1993294/3‑148 (MQ=2)
                                                                     tCAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCgg                                                                                 >  1:3792496/1‑150 (MQ=2)
                                                                      cAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGa                                                                                <  1:1764462/150‑1 (MQ=2)
                                                                      cAATGAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGa                                                                                <  2:4170473/150‑1 (MQ=2)
                                                                      cAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGa                                                                                >  1:869204/1‑150 (MQ=255)
                                                                       aatgAAGTGACCCGTCATTCGGCATGGGACCGCATGAGTGCAGCCGGTGCGCAGCTGATGACCTGGTTCGGCGTTGCTTGTGAACTGCATCGCGACTGGCGCAACGATATTGAAGGTCTGGGTACGCTGTTTGCAAACCACATTCCGGat                                                                               <  1:2996529/147‑2 (MQ=2)
                                                                       aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                                                                               >  1:5105800/1‑150 (MQ=255)
                                                                       aaCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGAc                                                                               >  1:1358966/1‑150 (MQ=255)
                                                                        aCGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACt                                                                              >  2:208469/1‑150 (MQ=255)
                                                                          gAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTAc                                                                            >  2:3385532/1‑150 (MQ=255)
                                                                           aaCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTAcc                                                                           <  1:4501505/150‑1 (MQ=255)
                                                                            aCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                  <  2:1117567/110‑1 (MQ=255)
                                                                            aCTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAg                                                                                                                  >  1:1117567/1‑110 (MQ=255)
                                                                               gACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTACACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCACaa                                                                       >  1:545994/1‑150 (MQ=255)
                                                                                aCCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAAc                                                                      >  2:406948/1‑150 (MQ=255)
                                                                                 cccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAAcc                                                                     >  2:4536301/1‑150 (MQ=255)
                                                                                  ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt                                                                    >  2:4145747/1‑150 (MQ=255)
                                                                                  ccGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCt                                                                    >  1:1350917/1‑150 (MQ=255)
                                                                                    gTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtga                                                                  <  2:3827102/150‑1 (MQ=255)
                                                                                    gTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtga                                                                  <  2:577744/150‑1 (MQ=255)
                                                                                     tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                 <  2:1474045/150‑1 (MQ=255)
                                                                                     tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                 <  2:5013918/150‑1 (MQ=255)
                                                                                     tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                 <  1:1765973/150‑1 (MQ=255)
                                                                                     tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                 <  1:1839265/150‑1 (MQ=255)
                                                                                     tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                 >  2:1375630/1‑150 (MQ=255)
                                                                                     tGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCtgat                                                                 >  2:1840083/1‑150 (MQ=255)
                                                                                         tCGGCCTGGCGACGCATGGAGGCTGCCGGGGCACAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGAcc                                                             >  2:4269846/1‑150 (MQ=255)
                                                                                          cGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCa                                                            <  1:4034069/150‑1 (MQ=255)
                                                                                           ggCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAg                                                           <  2:2603396/150‑1 (MQ=255)
                                                                                              cTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTa                                                        <  2:2374112/150‑1 (MQ=255)
                                                                                                ggCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTaca                                                      >  1:1842447/1‑150 (MQ=255)
                                                                                                  cGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTacaac                                                    >  1:5040482/1‑150 (MQ=255)
                                                                                                   gACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTacaaca                                                   >  2:4266242/1‑150 (MQ=255)
                                                                                                      gCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGc                                                >  2:2480328/1‑150 (MQ=255)
                                                                                                      gCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATAGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGc                                                >  2:4304505/1‑150 (MQ=255)
                                                                                                       cATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCt                                               <  2:4805599/150‑1 (MQ=255)
                                                                                                             ggCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaa                                         <  1:1988901/150‑1 (MQ=255)
                                                                                                              gCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCaaa                                        <  1:2356584/150‑1 (MQ=255)
                                                                                                                tGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATa                                      <  1:4299568/150‑1 (MQ=255)
                                                                                                                           aaCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgg                           >  2:1663223/1‑150 (MQ=255)
                                                                                                                            aCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggg                          <  1:1139431/150‑1 (MQ=255)
                                                                                                                            aCTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggg                          <  2:2122097/150‑1 (MQ=255)
                                                                                                                             cTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTgggg                         <  1:4515825/150‑1 (MQ=255)
                                                                                                                              tgatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggggg                        >  1:701835/1‑150 (MQ=255)
                                                                                                                              tgatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggggg                        >  1:3596790/1‑150 (MQ=255)
                                                                                                                              tgatgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTggggg                        <  1:2869567/150‑1 (MQ=255)
                                                                                                                                atgaCCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCt                      <  1:3028781/150‑1 (MQ=255)
                                                                                                                                    ccTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCAc                  <  2:2243626/150‑1 (MQ=255)
                                                                                                                                     cTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACg                 >  2:1249645/1‑150 (MQ=255)
                                                                                                                                     cTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACg                 >  1:1324203/1‑150 (MQ=255)
                                                                                                                                        gTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACa              >  2:432210/1‑150 (MQ=255)
                                                                                                                                             gCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGcccc         >  2:4505902/1‑150 (MQ=255)
                                                                                                                                                tGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCcac      >  2:4929935/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTTCCCGGCCCTGTCGGCATTGGAGGAGGGCTTCGAGGTATTCGTCGTCACCGACGCCTCTGGCACGTTCAACGAACTGACCCGTGACTCGGCCTGGCGACGCATGGAGGCTGCCGGGGCGCAACTGATGACCTGGTTCGGCGTGGCCTGCGAACTGCACCGTGACTGGCGCAACGACATCGAAGGGCTGGGTACCCTGTTCTCCAACCATATCCCGGACTACCGCAACCTGATGACCAGCTACAACAAGCTGGCCAAATAGTGCCTTCTGGGGGCTGCACGACAGCCCCCACATCA  >  NC_002947/827171‑827468

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