Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,289,287 C→A 7.8% P171P (CCG→CCT)  ccmC ← heme exporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,289,2870CA 7.8% 1274.6 / 64.8 340P171P (CCG→CCT) ccmCheme exporter subunit
Reads supporting (aligned to +/- strand):  ref base C (165/148);  new base A (10/17);  total (175/165)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 343 total aligned reads displayed.

ATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACC  >  CP009273/2289138‑2289435
                                                                                                                                                     |                                                                                                                                                    
    gcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGTATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCa                                                                                                                                                  <  2:1733801/150‑1 (MQ=255)
     catcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAg                                                                                                                                                 >  1:4076265/1‑150 (MQ=255)
       tcagcgtcagcgTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGAt                                                                                                                                               <  2:2161354/150‑1 (MQ=255)
       tcagcgtcagcgTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGGTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGAt                                                                                                                                               >  2:3704377/1‑150 (MQ=255)
        cagcgtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAATGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGAtt                                                                                                                                              >  2:864634/1‑150 (MQ=255)
        cagcgtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAATGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGGTCCACCACTCCACGGAGTAATGACTAATCGGCAGAtt                                                                                                                                              >  2:4031147/1‑150 (MQ=255)
        cagcgtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGAtt                                                                                                                                              >  1:3524891/1‑150 (MQ=255)
        cagcgtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGAtt                                                                                                                                              <  1:4092863/150‑1 (MQ=255)
         agcgtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTc                                                                                                                                             <  2:585828/150‑1 (MQ=255)
         agcgtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTc                                                                                                                                             <  1:1557038/150‑1 (MQ=255)
           cgtcagcgtTGCAGACAGGTGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcac                                                                                                                                           <  2:5067741/150‑1 (MQ=255)
             tcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcacca                                                                                                                                         >  2:3484602/1‑150 (MQ=255)
               agcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACg                                                                                                                                       >  1:1705393/1‑150 (MQ=255)
               agcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACg                                                                                                                                       <  1:4465576/150‑1 (MQ=255)
                 cgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATCCTTTGCTGCATCCGCGGTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGcc                                                                                                                                     >  1:2658288/1‑150 (MQ=255)
                   tGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCaa                                                                                                                                   >  1:3282592/1‑150 (MQ=255)
                    ggCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAAt                                                                                                                                  <  2:4661376/150‑1 (MQ=255)
                    ggCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAAt                                                                                                                                  <  1:4646480/150‑1 (MQ=255)
                     gCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATc                                                                                                                                 >  2:3188804/1‑150 (MQ=255)
                       aGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                                                                                                                               >  1:2229486/1‑150 (MQ=255)
                       aGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                                                                                                                               >  2:4986690/1‑150 (MQ=255)
                        gACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGc                                                                                                                              >  2:2486540/1‑150 (MQ=255)
                        gACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGGAATGAATAATCGGCAGATTCACCACGCCAATCAGc                                                                                                                              >  1:1316943/1‑150 (MQ=255)
                         aCAGGAGCAGGAAGCCAAACATCGTCCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCa                                                                                                                             <  2:3665684/150‑1 (MQ=255)
                          cAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAAGCAGCAc                                                                                                                            >  2:4461810/1‑150 (MQ=255)
                           aGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCAcc                                                                                                                           >  1:2902394/1‑150 (MQ=255)
                             gAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAg                                                                                                                         <  1:3716388/150‑1 (MQ=255)
                                 aGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATa                                                                                                                     >  2:3126311/1‑150 (MQ=255)
                                   gAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATAcc                                                                                                                   >  1:2781168/1‑150 (MQ=255)
                                          aaaaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCa                                                                                                            <  2:4055447/150‑1 (MQ=255)
                                            aaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACg                                                                                                          <  1:1150139/150‑1 (MQ=255)
                                            aaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACg                                                                                                          <  1:2546083/150‑1 (MQ=255)
                                             aTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACgg                                                                                                         >  1:5173526/1‑150 (MQ=255)
                                               cGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGcc                                                                                                       >  2:4133334/1‑150 (MQ=255)
                                                gACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGccc                                                                                                      >  1:4235984/1‑150 (MQ=255)
                                                  ccAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGc                                                                                                    >  2:2621570/1‑150 (MQ=255)
                                                    agcgcagcgGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCa                                                                                                  >  2:4071410/1‑150 (MQ=255)
                                                         agcgGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcgg                                                                                             >  1:4879841/1‑150 (MQ=255)
                                                          gcggcgAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcggc                                                                                            <  1:5464624/150‑1 (MQ=255)
                                                            ggcgAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcggcgg                                                                                          >  1:2404560/1‑150 (MQ=255)
                                                            ggcgAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGAcggcgg                                                                                          >  1:2069968/1‑150 (MQ=255)
                                                              cgAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGgtc                                                                                        <  1:2090240/150‑1 (MQ=255)
                                                                aaCGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGgtcgt                                                                                      >  1:4950637/1‑150 (MQ=255)
                                                                 aCGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGgtcgtc                                                                                     <  1:3125357/150‑1 (MQ=255)
                                                                  cGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGtcgtcg                                                                                    <  1:4328585/150‑1 (MQ=255)
                                                                   gCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGa                                                                                   >  1:3434610/1‑150 (MQ=255)
                                                                    cATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGaa                                                                                  <  1:2721605/150‑1 (MQ=255)
                                                                      tCGCCGGATCTATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAgg                                                                                <  2:849934/150‑1 (MQ=255)
                                                                         ccGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGt                                                                             <  2:4197527/150‑1 (MQ=255)
                                                                            gATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGcc                                                                          >  2:3400056/1‑150 (MQ=255)
                                                                            gATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGcc                                                                          <  1:1098508/150‑1 (MQ=255)
                                                                               cGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCcaca                                                                       <  2:1553677/150‑1 (MQ=255)
                                                                                gATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAg                                                                      >  2:2263710/1‑150 (MQ=255)
                                                                                  tACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACTCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAggg                                                                    >  1:2944496/1‑150 (MQ=255)
                                                                                   aCTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGc                                                                   <  1:4265639/150‑1 (MQ=255)
                                                                                         ctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATcac                                                             <  2:1369313/150‑1 (MQ=255)
                                                                                         ctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATcac                                                             <  1:3188082/150‑1 (MQ=255)
                                                                                          tgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATcaca                                                            <  2:295217/150‑1 (MQ=255)
                                                                                                  cgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATa                                                    <  2:431659/150‑1 (MQ=255)
                                                                                                  cgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATa                                                    <  2:2665055/150‑1 (MQ=255)
                                                                                                   gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCTCAGGGCAATCACACCCACATAc                                                   <  2:4419988/150‑1 (MQ=255)
                                                                                                   gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATAc                                                   >  1:1477636/1‑150 (MQ=255)
                                                                                                   gTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATAc                                                   <  1:995330/150‑1 (MQ=255)
                                                                                                    ttGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACa                                                  >  2:1930883/1‑150 (MQ=255)
                                                                                                       aTCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaa                                               >  2:1825803/1‑150 (MQ=255)
                                                                                                          ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAca                                            <  2:2073752/150‑1 (MQ=255)
                                                                                                            tGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcagc                                          >  1:3583413/1‑150 (MQ=255)
                                                                                                             gATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcagca                                         <  2:2789574/150‑1 (MQ=255)
                                                                                                              aTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcagcag                                        <  2:3775954/150‑1 (MQ=255)
                                                                                                              aTGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcagcag                                        <  1:1439098/150‑1 (MQ=255)
                                                                                                                 cAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAc                                     <  2:5111271/150‑1 (MQ=255)
                                                                                                                 cAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAc                                     <  1:3921671/150‑1 (MQ=255)
                                                                                                                  aGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAcc                                    >  1:175633/1‑150 (MQ=255)
                                                                                                                  aGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCAcc                                    >  1:806409/1‑150 (MQ=255)
                                                                                                                    ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAg                                  >  1:5293778/1‑150 (MQ=255)
                                                                                                                    ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAg                                  >  1:5293786/1‑150 (MQ=255)
                                                                                                                    ggTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAg                                  <  2:3282631/150‑1 (MQ=255)
                                                                                                                    ggTGTTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAg                                  <  2:4282230/150‑1 (MQ=255)
                                                                                                                      tgtTCCACCACTCCACGGAGTAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt                                <  2:3864152/150‑1 (MQ=255)
                                                                                                                        ttCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCa                              <  1:5386395/150‑1 (MQ=255)
                                                                                                                          ccaccaCTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATCCCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCACTCACACCCACATACAAAAACAGCAGCACCAGTTCAGa                            >  2:528867/1‑150 (MQ=255)
                                                                                                                           caccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGATGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGaa                           >  2:5389944/1‑150 (MQ=255)
                                                                                                                            accacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAg                          <  2:230598/150‑1 (MQ=255)
                                                                                                                               acTCCTCGGAGTAATGAATAATCGTCAGATTCACCACGCCAATCAGCACCAGGTTACCTGCCGCACGGCCCTCCAGACGGCGGTCGTCGAAGGGGTGCCTCTGGGCAATCTCACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCa                       <  1:4129034/150‑1 (MQ=255)
                                                                                                                                   cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACg                   <  1:3400044/150‑1 (MQ=255)
                                                                                                                                   cACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACg                   <  2:1238632/150‑1 (MQ=255)
                                                                                                                                    aCGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGt                  >  2:5028064/1‑150 (MQ=255)
                                                                                                                                     cGGAGTAATGAATACTCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTg                 <  2:1514625/150‑1 (MQ=255)
                                                                                                                                      ggAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGc                >  1:1955058/1‑150 (MQ=255)
                                                                                                                                      ggAGTAATGAATAATCGGCAGATTCACCACGCCAAACAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAAACACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGc                >  1:5481826/1‑150 (MQ=255)
                                                                                                                                       gAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCa               <  1:4593291/150‑1 (MQ=255)
                                                                                                                                          tAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc            <  1:2480693/150‑1 (MQ=255)
                                                                                                                                          tAATGAATAATAGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATcc            >  1:2120699/1‑150 (MQ=255)
                                                                                                                                           aaTGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATccc           >  2:526681/1‑150 (MQ=255)
                                                                                                                                            aTGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCa          <  2:2069968/150‑1 (MQ=255)
                                                                                                                                               aataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAc       <  1:1730119/150‑1 (MQ=255)
                                                                                                                                                ataatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAcc      >  2:4928762/1‑150 (MQ=255)
                                                                                                                                                ataatCGGCAGAGTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAcc      <  2:4990150/150‑1 (MQ=255)
                                                                                                                                                 taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAccc     >  2:2522389/1‑150 (MQ=255)
                                                                                                                                                 taatCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATAccc     <  2:4539322/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAAGTCAGACGTGCATCCCATACCCACC  >  CP009273/2289138‑2289435

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: ‑