Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,335,651 C→A 8.6% S7* (TCG→TAG)  kdsC → 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,335,6510CA 8.6% 1524.1 / 73.6 407S7* (TCG→TAG) kdsC3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (187/185);  new base A (12/23);  total (199/208)
Fisher's exact test for biased strand distribution p-value = 7.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 413 total aligned reads displayed.

CACTGAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTC  >  CP009273/3335502‑3335800
                                                                                                                                                     |                                                                                                                                                     
   tGAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCt                                                                                                                                                    >  1:3038701/1‑150 (MQ=255)
   tGAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCt                                                                                                                                                    >  2:1060515/1‑150 (MQ=255)
     aaCTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTg                                                                                                                                                  <  1:1326808/150‑1 (MQ=255)
        ttAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGa                                                                                                                                               <  1:1235525/150‑1 (MQ=255)
            tGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTg                                                                                                                                           <  1:2059142/150‑1 (MQ=255)
              cAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGtt                                                                                                                                         <  1:3559592/150‑1 (MQ=255)
                 tCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACg                                                                                                                                      >  1:261330/1‑150 (MQ=255)
                 tCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACg                                                                                                                                      >  2:4777179/1‑150 (MQ=255)
                    cGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGcc                                                                                                                                   <  1:67796/150‑1 (MQ=255)
                      ccATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTTCTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCt                                                                                                                                 <  2:2444758/150‑1 (MQ=255)
                        atatCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGt                                                                                                                               >  1:693453/1‑150 (MQ=255)
                         tttctccTCCGTGATGGTTGCCGATGGAGACCTTTTACTCGTTGTGTTACATATGCTTGATTTTCTTCGTGCTGGCGTAGTGTAAAGATTCAAGGATTAACATCAATGAGCATAGCAGGTGCGTCTCTTGCGACCTGTTACGGCCCTGTc                                                                                                                              <  2:2843301/145‑1 (MQ=255)
                          atCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCa                                                                                                                             <  2:1544216/150‑1 (MQ=255)
                          atCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCa                                                                                                                             >  2:5020009/1‑150 (MQ=255)
                          atCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCa                                                                                                                             <  1:2887525/150‑1 (MQ=255)
                            cACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAgc                                                                                                                           <  2:3423094/150‑1 (MQ=255)
                            cACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAgc                                                                                                                           >  1:4924854/1‑150 (MQ=255)
                            cACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAgc                                                                                                                           >  1:2079033/1‑150 (MQ=255)
                               cTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGcc                                                                                                                        >  1:1459369/1‑150 (MQ=255)
                                tCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCg                                                                                                                       >  2:116376/1‑150 (MQ=255)
                                  cGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGAc                                                                                                                     <  1:5082030/150‑1 (MQ=255)
                                  cGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGAc                                                                                                                     >  1:4487442/1‑150 (MQ=255)
                                  cGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGAc                                                                                                                     <  2:7131/150‑1 (MQ=255)
                                   gTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACg                                                                                                                    >  1:2792570/1‑150 (MQ=255)
                                   gTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACg                                                                                                                    >  1:4666636/1‑150 (MQ=255)
                                   gTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACg                                                                                                                    >  2:1306923/1‑150 (MQ=255)
                                   gTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACg                                                                                                                    >  2:4826225/1‑150 (MQ=255)
                                     gATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGtt                                                                                                                  <  2:4989026/150‑1 (MQ=255)
                                        ggTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTtata                                                                                                               <  1:936208/150‑1 (MQ=255)
                                               gATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAg                                                                                                        <  2:4724830/150‑1 (MQ=255)
                                                aTGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGc                                                                                                       >  2:4535386/1‑150 (MQ=255)
                                                aTGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGc                                                                                                       >  2:1271034/1‑150 (MQ=255)
                                                 tGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCa                                                                                                      <  2:2329502/150‑1 (MQ=255)
                                                  ggCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCag                                                                                                     >  2:539339/1‑150 (MQ=255)
                                                  ggCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCACTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCag                                                                                                     >  1:1842642/1‑150 (MQ=255)
                                                    cGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACACTGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCagag                                                                                                   <  2:2066631/150‑1 (MQ=255)
                                                     gACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAgaga                                                                                                  >  1:4606889/1‑150 (MQ=255)
                                                      aCCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGaa                                                                                                 >  2:584745/1‑150 (MQ=255)
                                                      aCCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGGTATAGCAAAAGCAGAGaa                                                                                                 >  2:300731/1‑150 (MQ=255)
                                                          tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTTCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             <  1:3900892/150‑1 (MQ=255)
                                                          tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             <  2:1895816/150‑1 (MQ=255)
                                                          tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             <  2:922189/150‑1 (MQ=255)
                                                          tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             <  1:4624469/150‑1 (MQ=255)
                                                          tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             <  1:3409283/150‑1 (MQ=255)
                                                           ttACTCGGTGTGTTACATAGGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTc                                                                                            <  2:4681974/150‑1 (MQ=255)
                                                             aCTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGt                                                                                          <  2:656988/150‑1 (MQ=255)
                                                             aCTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGt                                                                                          >  2:3162864/1‑150 (MQ=255)
                                                              cTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTc                                                                                         >  1:1046490/1‑150 (MQ=255)
                                                              cTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTc                                                                                         >  2:1845678/1‑150 (MQ=255)
                                                                    gtgtTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGa                                                                                   <  1:1257761/150‑1 (MQ=255)
                                                                    gtgtTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGa                                                                                   >  2:423463/1‑150 (MQ=255)
                                                                     tgtTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGAt                                                                                  <  2:2867945/150‑1 (MQ=255)
                                                                       ttACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATcc                                                                                >  2:1248465/1‑150 (MQ=255)
                                                                             atGCATGATTTACTGCGTGCAGGTGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcgatg                                                                          >  1:4253011/1‑150 (MQ=255)
                                                                                  tGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcgatgtcgat                                                                     <  2:366749/150‑1 (MQ=255)
                                                                                  tGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcgatgtcgat                                                                     <  1:4317859/150‑1 (MQ=255)
                                                                                   gATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcgatgtcgatg                                                                    >  1:2246664/1‑150 (MQ=255)
                                                                                   gATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcgatgtcgatg                                                                    <  2:261330/150‑1 (MQ=255)
                                                                                    aTTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATgg                                                                   >  2:5029226/1‑150 (MQ=255)
                                                                                          tGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACt                                                             >  1:274778/1‑150 (MQ=255)
                                                                                            cGTGCAGGCGTAGTGTAAAGATTCAAGGATTAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGt                                                           <  2:2686513/150‑1 (MQ=255)
                                                                                            cGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGt                                                           >  1:2129004/1‑150 (MQ=255)
                                                                                             gTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTc                                                          >  1:540085/1‑150 (MQ=255)
                                                                                                 aGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGAt                                                      >  1:2038095/1‑150 (MQ=255)
                                                                                                 aGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGAt                                                      <  1:3463316/150‑1 (MQ=255)
                                                                                                 aGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGAt                                                      >  2:1601760/1‑150 (MQ=255)
                                                                                                  ggCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATg                                                     <  1:90685/150‑1 (MQ=255)
                                                                                                   gCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATgg                                                    >  2:2790225/1‑150 (MQ=255)
                                                                                                    cGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGc                                                   >  1:4688962/1‑150 (MQ=255)
                                                                                                      tAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCt                                                 <  2:405162/150‑1 (MQ=255)
                                                                                                      tAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCt                                                 >  2:1780588/1‑150 (MQ=255)
                                                                                                        gtgtAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGa                                               <  2:3866456/150‑1 (MQ=255)
                                                                                                        gtgtAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGa                                               <  1:965886/150‑1 (MQ=255)
                                                                                                          gtAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGAtt                                             <  2:3356834/150‑1 (MQ=255)
                                                                                                           tAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGAttt                                            >  1:4044409/1‑150 (MQ=255)
                                                                                                             aaGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTtat                                          <  2:362164/150‑1 (MQ=255)
                                                                                                              aGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTtata                                         <  2:2468876/150‑1 (MQ=255)
                                                                                                                 ttCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATgg                                      >  1:1091718/1‑150 (MQ=255)
                                                                                                                  tCAAGGATACACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATggg                                     <  1:3769459/150‑1 (MQ=255)
                                                                                                                   cAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGc                                    >  2:4420250/1‑150 (MQ=255)
                                                                                                                    aaGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCa                                   >  1:4820492/1‑150 (MQ=255)
                                                                                                                    aaGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCa                                   <  1:3673581/150‑1 (MQ=255)
                                                                                                                      ggATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCaat                                 <  1:1112429/150‑1 (MQ=255)
                                                                                                                        aTAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCaataa                               >  2:3726959/1‑150 (MQ=255)
                                                                                                                          aaaCAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATg                             <  1:3111999/150‑1 (MQ=255)
                                                                                                                            actactttGAGCTTAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTTCTGTCAGATGGCCTGATTTATATGGGCAATAATGGc                           <  1:2965081/143‑1 (MQ=255)
                                                                                                                            acaacaATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGc                           <  2:2129269/150‑1 (MQ=255)
                                                                                                                            acaacaATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGc                           >  2:3494398/1‑150 (MQ=255)
                                                                                                                              aacaaTGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGa                         <  1:3339356/150‑1 (MQ=255)
                                                                                                                                caaTGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAag                       <  1:2410909/150‑1 (MQ=255)
                                                                                                                                   tGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGc                    <  1:4535851/150‑1 (MQ=255)
                                                                                                                                    gAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCt                   <  2:5063190/150‑1 (MQ=255)
                                                                                                                                      gCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGa                 >  1:1164673/1‑150 (MQ=255)
                                                                                                                                      gCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTGATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGGACTGTCAGATGGCCGGATTTAGATGGGCACTCATGGCGAAGAGCTGa                 >  2:1987124/1‑150 (MQ=255)
                                                                                                                                        aaaGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGaaa               >  1:4169182/1‑150 (MQ=255)
                                                                                                                                        aaaGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGaaa               >  2:823568/1‑150 (MQ=255)
                                                                                                                                            cAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGt           >  2:1380823/1‑150 (MQ=255)
                                                                                                                                              ggTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTc         >  1:1670900/1‑150 (MQ=255)
                                                                                                                                                  cGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCACAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATg     <  2:3888851/150‑1 (MQ=255)
                                                                                                                                                     cGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTc  >  1:3261425/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACTGAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTC  >  CP009273/3335502‑3335800

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: