Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,563,276 G→T 100% 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,2760GT96.9% 805.2 / 13.6 245H491N (CAC→AAC) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base T (148/89);  total (150/95)
Fisher's exact test for biased strand distribution p-value = 5.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 245 total aligned reads displayed.

CCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGG  >  CP009273/1563129‑1563422
                                                                                                                                                   |                                                                                                                                                  
cccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTAc                                                                                                                                                  >  1:3524161/1‑150 (MQ=255)
     tAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCTGACATGGTTCGTGTTTCTTGTCATTCATCAGAATATAGTGTGGTGAACGTGCACGCGGGTGCTGGAAGACGTTAccgtgt                                                                                                                                             >  2:2540611/1‑146 (MQ=255)
     tAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtttt                                                                                                                                             <  2:2766673/150‑1 (MQ=255)
        gTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGt                                                                                                                                          <  2:4085742/150‑1 (MQ=255)
         tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATCACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTg                                                                                                                                         <  1:400752/150‑1 (MQ=255)
            aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGt                                                                                                                                      >  2:1641023/1‑150 (MQ=255)
            aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGt                                                                                                                                      <  1:998005/150‑1 (MQ=255)
               cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGCAGAAGTTACCGTTTGGTGCGTTCt                                                                                                                                   >  2:337124/1‑150 (MQ=255)
               cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCt                                                                                                                                   >  1:3479469/1‑150 (MQ=255)
                      aagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGtt                                                                                                                            >  1:2424805/1‑150 (MQ=255)
                      aagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGGGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTCCCTTTTGGTGCGTTCTGACTGtt                                                                                                                            >  2:199621/1‑150 (MQ=255)
                      aagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGGTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGtt                                                                                                                            >  2:1376177/1‑150 (MQ=255)
                       agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTg                                                                                                                           >  1:3846178/1‑150 (MQ=255)
                       agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTg                                                                                                                           >  1:2693082/1‑150 (MQ=255)
                         caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGAt                                                                                                                         >  1:4612766/1‑150 (MQ=255)
                         caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                                                                                                                         <  2:2872460/150‑1 (MQ=255)
                         caCGAAAAAGGGAGCGATGAATGATCGCTCCCTTGTCTTATCACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGAt                                                                                                                         <  2:3096059/150‑1 (MQ=255)
                            gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTgg                                                                                                                      >  1:3214525/1‑150 (MQ=255)
                              aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                    >  2:251992/1‑150 (MQ=255)
                                aaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTcc                                                                                                                  >  2:2070451/1‑150 (MQ=255)
                                 aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 <  1:2597928/150‑1 (MQ=255)
                                 aGGGAGCGATGAATTATCGCTCACTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTTGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 <  1:3870394/150‑1 (MQ=255)
                                   ggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga                                                                                                               >  1:1499142/1‑150 (MQ=255)
                                    gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag                                                                                                              <  2:4503106/150‑1 (MQ=255)
                                    gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag                                                                                                              >  1:3251297/1‑150 (MQ=255)
                                    gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag                                                                                                              <  2:2318940/150‑1 (MQ=255)
                                     tgCGATGATTTTTCGCTCCCTTGTCTTATCCCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGTTGAACGTGCACGCGGGTGCAGTTAGAAGTTACCTTTTGGTGAGTTCTGTCTGTTGATTGGCTCCAGAGt                                                                                                             <  2:263325/149‑1 (MQ=255)
                                       cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCGGAGTGa                                                                                                           >  2:3540701/1‑150 (MQ=255)
                                        gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc                                                                                                          >  2:1677764/1‑150 (MQ=255)
                                        gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc                                                                                                          <  1:1416112/150‑1 (MQ=255)
                                           gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc                                                                                                       <  1:3688400/150‑1 (MQ=255)
                                           gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc                                                                                                       >  2:4172240/1‑150 (MQ=255)
                                            aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg                                                                                                      >  2:4276012/1‑150 (MQ=255)
                                            aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg                                                                                                      >  1:4499976/1‑150 (MQ=255)
                                             aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                     >  1:2009200/1‑150 (MQ=255)
                                             aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                     >  2:2523879/1‑150 (MQ=255)
                                                 tCGCTCCCTTGTATTATAACCATTCAGACAAGGTTAGTGTTTCTTGTCATTCATCACAATATAGAGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTTACTGTTGATTGGCTCCAGAGTGACTCCGGTAtt                                                                                                 >  2:467459/1‑150 (MQ=255)
                                                  cGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAttt                                                                                                <  1:2138076/150‑1 (MQ=255)
                                                     tCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgct                                                                                             >  1:3540948/1‑150 (MQ=255)
                                                     tCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgct                                                                                             >  1:4368004/1‑150 (MQ=255)
                                                        cTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgctttg                                                                                          >  1:722622/1‑150 (MQ=255)
                                                        cTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGGTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgctttg                                                                                          <  1:1243593/150‑1 (MQ=255)
                                                             cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                     >  1:2205032/1‑150 (MQ=255)
                                                             cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                     <  1:4416692/150‑1 (MQ=255)
                                                             cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                     >  2:3836813/1‑150 (MQ=255)
                                                              ttATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                    <  2:2151639/150‑1 (MQ=255)
                                                              ttATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                    >  1:3230393/1‑150 (MQ=255)
                                                              ttATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                    <  2:4198568/150‑1 (MQ=255)
                                                               tataACCATTCAGACATTGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTAc                                                                                   >  1:180938/1‑150 (MQ=255)
                                                               tataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTAc                                                                                   >  1:756357/1‑150 (MQ=255)
                                                                  aaCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGAt                                                                                >  2:4465635/1‑150 (MQ=255)
                                                                  aaCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGAt                                                                                >  2:3312721/1‑150 (MQ=255)
                                                                  aaCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGAt                                                                                <  2:901329/150‑1 (MQ=255)
                                                                      aTTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATg                                                                            >  1:2794859/1‑150 (MQ=255)
                                                                       ttCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGa                                                                           >  2:2238435/1‑150 (MQ=255)
                                                                       ttCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGa                                                                           >  2:2238329/1‑150 (MQ=255)
                                                                          aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                        >  2:3076440/1‑150 (MQ=255)
                                                                            aCATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGc                                                                      >  2:641973/1‑150 (MQ=255)
                                                                               tGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATa                                                                   <  1:4522365/150‑1 (MQ=255)
                                                                                ggTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATag                                                                  <  2:1573738/150‑1 (MQ=255)
                                                                                   tAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAgaga                                                               <  2:1772756/150‑1 (MQ=255)
                                                                                   tAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAgaga                                                               >  1:608110/1‑150 (MQ=255)
                                                                                   tAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAgaga                                                               >  1:2404593/1‑150 (MQ=255)
                                                                                     gtgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAAt                                                             >  2:2082147/1‑150 (MQ=255)
                                                                                           cTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAaggg                                                       <  2:4127621/150‑1 (MQ=255)
                                                                                             tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggca                                                     <  1:1117485/150‑1 (MQ=255)
                                                                                             tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggca                                                     >  1:4130412/1‑150 (MQ=255)
                                                                                                 attcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcagggc                                                 >  2:506643/1‑150 (MQ=255)
                                                                                                  ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAgggcagggca                                                <  2:3160010/150‑1 (MQ=255)
                                                                                                  ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAgggcagggca                                                <  1:2797591/150‑1 (MQ=255)
                                                                                                       cacaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTac                                           <  1:4245735/150‑1 (MQ=255)
                                                                                                         caATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTacca                                         <  2:16133/150‑1 (MQ=255)
                                                                                                          aaTATAGTGTGGTGAACGTGCACGAGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTaccac                                        <  2:375655/150‑1 (MQ=255)
                                                                                                              taGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAgg                                    >  2:3112081/1‑150 (MQ=255)
                                                                                                              taGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAgg                                    <  1:4401062/150‑1 (MQ=255)
                                                                                                               aGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGa                                   >  1:4348719/1‑150 (MQ=255)
                                                                                                               aGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGa                                   >  1:4349231/1‑150 (MQ=255)
                                                                                                                gtgtGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGaa                                  <  2:4379332/150‑1 (MQ=255)
                                                                                                                   tggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACt                               >  2:2426497/1‑150 (MQ=255)
                                                                                                                    ggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTa                              <  1:3015493/150‑1 (MQ=255)
                                                                                                                       gAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaaca                           >  2:2407330/1‑150 (MQ=255)
                                                                                                                        aaCGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaacaa                          >  1:542923/1‑150 (MQ=255)
                                                                                                                          cGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAAcc                        <  2:2969121/150‑1 (MQ=255)
                                                                                                                           gTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                       >  2:1785885/1‑150 (MQ=255)
                                                                                                                             gCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGt                     >  1:3064334/1‑150 (MQ=255)
                                                                                                                             gCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAAGAAACTAACAACCAGt                     <  1:602954/150‑1 (MQ=255)
                                                                                                                                cgggggTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAt                  <  1:2345551/147‑1 (MQ=255)
                                                                                                                                cgcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATTAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAt                  <  1:364934/150‑1 (MQ=255)
                                                                                                                                 gcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAtt                 <  1:3507541/150‑1 (MQ=255)
                                                                                                                                  cgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTc                <  1:501541/150‑1 (MQ=255)
                                                                                                                                     gTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTGTACGATCATGAACAGCCTAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAATAAac             >  2:4511754/1‑150 (MQ=255)
                                                                                                                                          ggAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacataca        >  1:4578784/1‑150 (MQ=255)
                                                                                                                                           gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatacat       <  1:1014875/150‑1 (MQ=255)
                                                                                                                                           gaagaaGTGACCTATTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatacat       <  1:2048058/150‑1 (MQ=255)
                                                                                                                                             agaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACatat     >  2:1105675/1‑150 (MQ=255)
                                                                                                                                             agaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACatat     >  2:1169303/1‑150 (MQ=255)
                                                                                                                                               aagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCg   <  2:1416646/150‑1 (MQ=255)
                                                                                                                                               aagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCg   <  2:870468/150‑1 (MQ=255)
                                                                                                                                               aagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCg   <  1:1045460/150‑1 (MQ=255)
                                                                                                                                                agTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCgg  >  2:621099/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGG  >  CP009273/1563129‑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: