Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,563,276 G→T 93.8% H491N (CAC→AAC)  gadC ← glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,563,2760GT93.8% 706.4 / 34.9 226H491N (CAC→AAC) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (11/3);  new base T (98/114);  total (109/117)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 226 total aligned reads displayed.

CGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGG  >  CP009273/1563131‑1563422
                                                                                                                                                 |                                                                                                                                                  
 gTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtt                                                                                                                                               <  2:3430139/150‑1 (MQ=255)
 gTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtt                                                                                                                                               <  1:3706489/150‑1 (MQ=255)
  ttAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCttt                                                                                                                                              >  1:4120775/1‑150 (MQ=255)
  ttAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCttt                                                                                                                                              <  1:2981264/150‑1 (MQ=255)
            aCGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTc                                                                                                                                    >  1:536941/1‑150 (MQ=255)
            aCGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTc                                                                                                                                    >  1:1785863/1‑150 (MQ=255)
            aCGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTc                                                                                                                                    <  2:1409458/150‑1 (MQ=255)
             cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCt                                                                                                                                   <  1:171780/150‑1 (MQ=255)
             cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCt                                                                                                                                   <  1:2136045/150‑1 (MQ=255)
             cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                                                                                                                                   >  1:920906/1‑150 (MQ=255)
             cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                                                                                                                                   >  1:3362854/1‑150 (MQ=255)
              gACAGCAAGCACGAAAAATGTAGCGATGAATTATCGCTCACTTTTCTTATAACCATTCAGACATGGTTAGTGTTTATTGTCATTCATCACAATATAGTTTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGAGTTCTg                                                                                                                                  <  1:4260133/150‑1 (MQ=255)
                 agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACt                                                                                                                               <  2:3901801/150‑1 (MQ=255)
                        aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGACGTTACCTTTTGGTGCGTTCTGACTGTGGAtt                                                                                                                        >  2:45676/1‑150 (MQ=255)
                        aCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGAtt                                                                                                                        <  2:1103562/150‑1 (MQ=255)
                          gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTgg                                                                                                                      >  1:2584971/1‑150 (MQ=255)
                           aaaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGc                                                                                                                     >  2:1764310/1‑150 (MQ=255)
                            aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                    >  2:3023042/1‑150 (MQ=255)
                             aaaGGGAGCTATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc                                                                                                                   >  1:3650305/1‑150 (MQ=255)
                              aaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTcc                                                                                                                  <  1:2649149/150‑1 (MQ=255)
                               aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 <  2:835459/150‑1 (MQ=255)
                               aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 <  2:835401/150‑1 (MQ=255)
                                  gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag                                                                                                              >  1:2767616/1‑150 (MQ=255)
                                     cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa                                                                                                           <  1:3361579/150‑1 (MQ=255)
                                      gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc                                                                                                          <  1:1912745/150‑1 (MQ=255)
                                        tGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGc                                                                                                        >  2:3552258/1‑150 (MQ=255)
                                          aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg                                                                                                      <  1:662302/150‑1 (MQ=255)
                                          aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg                                                                                                      >  2:2409498/1‑150 (MQ=255)
                                            ttATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGt                                                                                                    <  1:1471920/150‑1 (MQ=255)
                                              aTCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAt                                                                                                  <  1:1614404/150‑1 (MQ=255)
                                               tCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTCtt                                                                                                 >  2:4274312/1‑150 (MQ=255)
                                                cGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAttt                                                                                                >  2:1393078/1‑150 (MQ=255)
                                                cGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAttt                                                                                                <  1:1023508/150‑1 (MQ=255)
                                                   tCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgct                                                                                             >  2:2525853/1‑150 (MQ=255)
                                                        tGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGcc                                                                                        >  1:310570/1‑150 (MQ=255)
                                                         gTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCt                                                                                       <  1:3877821/150‑1 (MQ=255)
                                                          tCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCtt                                                                                      <  1:842299/150‑1 (MQ=255)
                                                          tCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCtt                                                                                      <  1:2411742/150‑1 (MQ=255)
                                                           cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                     >  1:2259289/1‑150 (MQ=255)
                                                               taACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGa                                                                                 >  1:4111420/1‑150 (MQ=255)
                                                                aaCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGAt                                                                                <  1:3266457/150‑1 (MQ=255)
                                                                 aCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATc                                                                               >  2:642605/1‑150 (MQ=255)
                                                                  ccATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCa                                                                              >  1:3969441/1‑150 (MQ=255)
                                                                  ccATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCa                                                                              <  1:1515973/150‑1 (MQ=255)
                                                                  ccATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGAGTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCa                                                                              <  2:2743723/150‑1 (MQ=255)
                                                                   cATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCAt                                                                             >  1:1166029/1‑150 (MQ=255)
                                                                   cATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCAt                                                                             <  2:3371004/150‑1 (MQ=255)
                                                                   cATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCAt                                                                             <  1:1932117/150‑1 (MQ=255)
                                                                      tCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGaa                                                                          <  2:3897407/150‑1 (MQ=255)
                                                                        aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                        <  2:2558642/150‑1 (MQ=255)
                                                                        aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACAGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                        >  1:2945596/1‑150 (MQ=255)
                                                                          aCATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGc                                                                      >  1:4156639/1‑150 (MQ=255)
                                                                                ttAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATagag                                                                >  1:2011499/1‑150 (MQ=255)
                                                                                  aGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGaa                                                              <  1:3327578/150‑1 (MQ=255)
                                                                                   gtgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGTGTGACGCCGGTATTTGCATTGCCTTTACGATCATGAACAGCATAGAGAAt                                                             <  1:903200/150‑1 (MQ=255)
                                                                                      tttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaaa                                                          >  2:153608/1‑150 (MQ=255)
                                                                                         cTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAaggg                                                       <  1:497026/150‑1 (MQ=255)
                                                                                           tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggca                                                     <  1:1839843/150‑1 (MQ=255)
                                                                                           tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggca                                                     <  2:1481141/150‑1 (MQ=255)
                                                                                            gTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcag                                                    >  1:822557/1‑150 (MQ=255)
                                                                                                   atcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                             >  2:4320793/1‑150 (MQ=255)
                                                                                                   atcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                             >  2:2245720/1‑150 (MQ=255)
                                                                                                   atcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                             >  2:2128802/1‑150 (MQ=255)
                                                                                                          tataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTAccacca                                      <  1:2149468/150‑1 (MQ=255)
                                                                                                          tataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTAccacca                                      >  1:4039704/1‑150 (MQ=255)
                                                                                                           ataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                     <  2:2374865/150‑1 (MQ=255)
                                                                                                           ataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                     >  1:2835662/1‑150 (MQ=255)
                                                                                                            taGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAgg                                    >  1:2207687/1‑150 (MQ=255)
                                                                                                            taGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAgg                                    <  1:947716/150‑1 (MQ=255)
                                                                                                            taGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGAATGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAgg                                    <  1:3131949/150‑1 (MQ=255)
                                                                                                              gtgtGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGaa                                  <  2:1978062/150‑1 (MQ=255)
                                                                                                              gtgtGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGaa                                  <  2:1463884/150‑1 (MQ=255)
                                                                                                                gtggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAAc                                >  2:2175635/1‑150 (MQ=255)
                                                                                                                 tggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACt                               <  2:2620102/150‑1 (MQ=255)
                                                                                                                  ggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTa                              >  2:3164571/1‑150 (MQ=255)
                                                                                                                  ggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTa                              >  1:2902510/1‑150 (MQ=255)
                                                                                                                    tgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaac                            >  1:2878616/1‑150 (MQ=255)
                                                                                                                     gAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCAGAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaaca                           >  2:930364/1‑150 (MQ=255)
                                                                                                                      aaCGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaacaa                          <  2:481927/150‑1 (MQ=255)
                                                                                                                      aaCGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaacaa                          >  1:3634016/1‑150 (MQ=255)
                                                                                                                         gTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                       >  2:3757300/1‑150 (MQ=255)
                                                                                                                         gTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                       <  1:3605701/150‑1 (MQ=255)
                                                                                                                               gcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAtt                 <  1:431900/150‑1 (MQ=255)
                                                                                                                               gcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAtt                 <  1:2603960/150‑1 (MQ=255)
                                                                                                                                cgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTc                >  1:2512386/1‑150 (MQ=255)
                                                                                                                                 gggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa               <  1:4145583/150‑1 (MQ=255)
                                                                                                                                  ggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCaa              >  2:2850605/1‑150 (MQ=255)
                                                                                                                                    tGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAaca            >  2:485970/1‑150 (MQ=255)
                                                                                                                                      cAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacata          >  1:2724124/1‑150 (MQ=255)
                                                                                                                                      cAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacata          >  2:732932/1‑150 (MQ=255)
                                                                                                                                        ggAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacataca        <  2:137928/150‑1 (MQ=255)
                                                                                                                                        ggAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacataca        >  2:3970477/1‑150 (MQ=255)
                                                                                                                                        ggAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacataca        >  2:289764/1‑150 (MQ=255)
                                                                                                                                        ggAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTGTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacataca        >  2:3337192/1‑150 (MQ=255)
                                                                                                                                         gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatacat       >  1:9774/1‑150 (MQ=255)
                                                                                                                                          aagaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATTAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAAcatacata      <  1:3023042/150‑1 (MQ=255)
                                                                                                                                           agaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACatat     <  1:1974213/150‑1 (MQ=255)
                                                                                                                                            gaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATc    <  2:4210606/150‑1 (MQ=255)
                                                                                                                                            gaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATc    <  1:3905120/150‑1 (MQ=255)
                                                                                                                                             aagTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCg   >  2:2318086/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGG  >  CP009273/1563131‑1563422

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: