Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,376,532 C→T 5.6% G886R (GGG→AGG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,376,5320CT 5.6% 1097.7 / 32.5 305G886R (GGG→AGG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (144/144);  new base T (13/4);  total (157/148)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

GCGAATCGGTTTTTCGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCA  >  CP009273/4376383‑4376680
                                                                                                                                                     |                                                                                                                                                    
  gAATCGGTTTTTCGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGc                                                                                                                                                    >  2:3933778/1‑150 (MQ=255)
        gTTTTTCGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatca                                                                                                                                              >  1:485350/1‑150 (MQ=255)
         tttttCGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatcag                                                                                                                                             <  1:2113576/150‑1 (MQ=255)
              cGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCTGCCAatcagcatca                                                                                                                                        >  1:382633/1‑150 (MQ=255)
              cGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatcagcatca                                                                                                                                        >  1:1894679/1‑150 (MQ=255)
               gAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatcagcatcag                                                                                                                                       <  1:3095077/150‑1 (MQ=255)
                aaCAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatcagcatcagc                                                                                                                                      <  1:153830/150‑1 (MQ=255)
                aaCAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatcagcatcagc                                                                                                                                      <  1:3468483/150‑1 (MQ=255)
                aaCAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAatcagcatcagc                                                                                                                                      >  1:2831270/1‑150 (MQ=255)
                 aCAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCTGCCACtcagcatcagca                                                                                                                                     >  2:2904875/1‑150 (MQ=255)
                     gATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGata                                                                                                                                 <  2:3710424/150‑1 (MQ=255)
                       tAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATAtt                                                                                                                               <  2:2040770/150‑1 (MQ=255)
                        aaTCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATAttt                                                                                                                              >  2:4559931/1‑150 (MQ=255)
                         aTCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtg                                                                                                                             >  2:1701376/1‑150 (MQ=255)
                           cAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggt                                                                                                                           >  2:1076460/1‑150 (MQ=255)
                            aGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtg                                                                                                                          >  2:664210/1‑150 (MQ=255)
                            aGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtg                                                                                                                          >  2:4113275/1‑150 (MQ=255)
                              aCCAGCGATACAGCTGGCGAAAATTTCCTGCACACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgat                                                                                                                        <  2:3606748/150‑1 (MQ=255)
                               ccAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatg                                                                                                                       <  2:2262048/150‑1 (MQ=255)
                               ccAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatg                                                                                                                       >  1:4287917/1‑150 (MQ=255)
                                 agagATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggt                                                                                                                     <  1:1032612/150‑1 (MQ=255)
                                  gagaTAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtg                                                                                                                    >  2:4644360/1‑150 (MQ=255)
                                   agaTAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtgg                                                                                                                   <  1:2362988/150‑1 (MQ=255)
                                    gaTAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtggt                                                                                                                  <  1:454884/150‑1 (MQ=255)
                                     aTAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCTGCCAATCAGCATCAGCAGATATTtggtgatggtggtg                                                                                                                 >  2:727839/1‑150 (MQ=255)
                                     aTAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtggtg                                                                                                                 >  2:3676913/1‑150 (MQ=255)
                                      tAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCTGCCAATCAGCATCAGCAGATATTtggtgatggtggtga                                                                                                                >  2:2160152/1‑150 (MQ=255)
                                       aaaGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtggtgat                                                                                                               >  1:2411346/1‑150 (MQ=255)
                                       aaaGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtggtgat                                                                                                               >  2:4591379/1‑150 (MQ=255)
                                       aaaGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtggtgat                                                                                                               >  2:2085997/1‑150 (MQ=255)
                                        aaGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTtggtgatggtggtgatg                                                                                                              <  2:1516496/150‑1 (MQ=255)
                                           ttGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCg                                                                                                           <  2:2698191/150‑1 (MQ=255)
                                                   aaTTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGt                                                                                                   >  1:1164301/1‑150 (MQ=255)
                                                    aTTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTa                                                                                                  >  1:4273201/1‑150 (MQ=255)
                                                    aTTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTa                                                                                                  <  1:1500226/150‑1 (MQ=255)
                                                     tttCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTAc                                                                                                 <  2:3407275/150‑1 (MQ=255)
                                                              aaCCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTa                                                                                        >  1:2013207/1‑150 (MQ=255)
                                                               aCCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTaa                                                                                       >  1:326816/1‑150 (MQ=255)
                                                                ccAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTaaa                                                                                      >  2:3539111/1‑150 (MQ=255)
                                                                 cAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAAt                                                                                     <  2:4505335/150‑1 (MQ=255)
                                                                 cAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAAt                                                                                     <  2:2640990/150‑1 (MQ=255)
                                                                 cAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAAt                                                                                     >  2:4190804/1‑150 (MQ=255)
                                                                  aaaaCCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCTGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATc                                                                                    >  2:1124518/1‑150 (MQ=255)
                                                                  aaaaCCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATc                                                                                    <  1:1076460/150‑1 (MQ=255)
                                                                  aaaaCCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATc                                                                                    >  1:4041911/1‑150 (MQ=255)
                                                                   aaaCCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATcc                                                                                   >  2:4477123/1‑150 (MQ=255)
                                                                      ccGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAgg                                                                                >  2:2434568/1‑150 (MQ=255)
                                                                       cGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGt                                                                               <  1:2952307/150‑1 (MQ=255)
                                                                        gagaCCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGtg                                                                              >  2:3741099/1‑150 (MQ=255)
                                                                        gagaCCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGtg                                                                              <  2:488928/150‑1 (MQ=255)
                                                                         agaCCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGtgc                                                                             >  2:1669012/1‑150 (MQ=255)
                                                                            ccAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGtgctgc                                                                          >  1:2568216/1‑150 (MQ=255)
                                                                             cAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCTGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCa                                                                         <  1:2669440/150‑1 (MQ=255)
                                                                             cAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCa                                                                         >  1:1901666/1‑150 (MQ=255)
                                                                              aaCACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCaa                                                                        >  2:1397078/1‑150 (MQ=255)
                                                                              aaCACCGAGCGCGGCAACCAGCCACTGCAATATCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCaa                                                                        <  1:1629425/150‑1 (MQ=255)
                                                                                  ccGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATc                                                                    <  2:4702765/150‑1 (MQ=255)
                                                                                     aGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGcc                                                                 >  2:380063/1‑150 (MQ=255)
                                                                                      gcgcGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCa                                                                <  1:2476966/150‑1 (MQ=255)
                                                                                      gcgcGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCa                                                                <  1:2434568/150‑1 (MQ=255)
                                                                                       cgcgGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAg                                                               >  2:1221499/1‑150 (MQ=255)
                                                                                          ggCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTc                                                            >  2:3301926/1‑150 (MQ=255)
                                                                                           gCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCCTTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTcc                                                           <  2:3194723/150‑1 (MQ=255)
                                                                                           gCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTcc                                                           <  1:2146039/150‑1 (MQ=255)
                                                                                            cAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCca                                                          <  2:522899/150‑1 (MQ=255)
                                                                                            cAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCca                                                          <  2:523217/150‑1 (MQ=255)
                                                                                                 aGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCAGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCagcagc                                                     >  1:3325943/1‑150 (MQ=255)
                                                                                                   ccACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAgcgc                                                   >  1:228603/1‑150 (MQ=255)
                                                                                                   ccACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAgcgc                                                   >  1:3973600/1‑150 (MQ=255)
                                                                                                    cACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGcgcg                                                  >  2:2538691/1‑150 (MQ=255)
                                                                                                       tGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGc                                               >  1:1646589/1‑150 (MQ=255)
                                                                                                        gCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCa                                              <  2:3263618/150‑1 (MQ=255)
                                                                                                         cAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCaa                                             <  2:4644646/150‑1 (MQ=255)
                                                                                                          aaTTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTATCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAg                                            <  1:4224217/150‑1 (MQ=255)
                                                                                                          aaTTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAg                                            <  1:2197291/150‑1 (MQ=255)
                                                                                                          aaTTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAg                                            <  1:866732/150‑1 (MQ=255)
                                                                                                          aaTTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAg                                            <  2:1740172/150‑1 (MQ=255)
                                                                                                            tttCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGtt                                          >  1:1292471/1‑150 (MQ=255)
                                                                                                              tCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTgc                                        <  1:3604984/150‑1 (MQ=255)
                                                                                                                    aCTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCa                                  >  2:4093878/1‑150 (MQ=255)
                                                                                                                      tCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGc                                >  2:4043132/1‑150 (MQ=255)
                                                                                                                       cAATACCAATCATTGAGAAGCCGACCAGCCTGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCt                               <  1:1443494/150‑1 (MQ=255)
                                                                                                                       aaaTACCAATCATTGAGAAGCCGACCAGCCTGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTGCCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCt                               <  1:4512188/149‑1 (MQ=255)
                                                                                                                        aaTACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTg                              <  1:3933778/150‑1 (MQ=255)
                                                                                                                         aTACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGCGGTTAATACCGGGTGCTCCAAATACGCCCGTACCACCACCCCGGGCACGTTGCCCCCCACCtcc                             >  2:3589274/1‑148 (MQ=255)
                                                                                                                          tACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCg                            >  1:1859693/1‑150 (MQ=255)
                                                                                                                             cAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCgtgg                         <  1:563916/150‑1 (MQ=255)
                                                                                                                             cAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCgtgg                         >  2:3907548/1‑150 (MQ=255)
                                                                                                                               aTCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCgtggtg                       <  2:1894681/150‑1 (MQ=255)
                                                                                                                                tCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGtga                      <  1:408514/150‑1 (MQ=255)
                                                                                                                                tCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGtga                      <  2:586460/150‑1 (MQ=255)
                                                                                                                                 cATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGtgat                     >  1:3217789/1‑150 (MQ=255)
                                                                                                                                  aTTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGtgatg                    <  2:2111371/150‑1 (MQ=255)
                                                                                                                                    tGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTgatgat                  <  2:4290512/150‑1 (MQ=255)
                                                                                                                                              aCCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGa        <  2:4377025/150‑1 (MQ=255)
                                                                                                                                              aCCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGa        <  2:3852342/150‑1 (MQ=255)
                                                                                                                                               ccAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTTGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGaa       <  1:3134160/150‑1 (MQ=255)
                                                                                                                                                cAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAAt      <  2:1416930/150‑1 (MQ=255)
                                                                                                                                                 aGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATg     >  1:3674103/1‑150 (MQ=255)
                                                                                                                                                 aGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATg     >  2:4667106/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGAATCGGTTTTTCGAACAGGATAATCAGACCAGAGATAAAGTTGGCGAAAATTTCCTGCAAACCAAAACCGAGACCAACACCGAGCGCGGCAACCAGCCACTGCAATTTCGACCACTCAATACCAATCATTGAGAAGCCGACCAGCCCGCCAATCAGCATCAGCAGATATTTGGTGATGGTGGTGATGGCGTAACCCGTACCCGGCGTTAAATCCAGGTGCTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCA  >  CP009273/4376383‑4376680

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