Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,563,276 G→T 86.3% 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,2760GT86.3% 623.6 / 101.8 242H491N (CAC→AAC) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base T (111/98);  total (130/112)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 242 total aligned reads displayed.

CCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGGTA  >  CP009273/1563129‑1563424
                                                                                                                                                   |                                                                                                                                                    
cccGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTAc                                                                                                                                                    >  1:1813652/1‑150 (MQ=255)
               cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCt                                                                                                                                     <  2:1668471/150‑1 (MQ=255)
               cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCt                                                                                                                                     <  2:2699985/150‑1 (MQ=255)
               cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCt                                                                                                                                     <  1:3957603/150‑1 (MQ=255)
                gACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTg                                                                                                                                    >  1:1700853/1‑150 (MQ=255)
                   agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACt                                                                                                                                 <  1:2361226/150‑1 (MQ=255)
                      aagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGACCGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGtt                                                                                                                              >  2:4387866/1‑150 (MQ=255)
                      aagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGtt                                                                                                                              >  2:2047265/1‑150 (MQ=255)
                            gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTgg                                                                                                                        >  1:807455/1‑150 (MQ=255)
                              aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                      <  2:4138813/150‑1 (MQ=255)
                              aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                      <  2:3105293/150‑1 (MQ=255)
                              aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                      <  1:2814203/150‑1 (MQ=255)
                                aaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTCCCTTTTGGTGCGTTCTGACTGTTGATTGGCTcc                                                                                                                    >  2:1784910/1‑150 (MQ=255)
                                aaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTcc                                                                                                                    >  1:1577992/1‑150 (MQ=255)
                                 aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                   >  1:3287415/1‑150 (MQ=255)
                                 aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                   >  1:407370/1‑150 (MQ=255)
                                   ggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga                                                                                                                 >  2:1868451/1‑150 (MQ=255)
                                    gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag                                                                                                                <  2:617245/150‑1 (MQ=255)
                                      gCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTg                                                                                                              <  2:4108370/150‑1 (MQ=255)
                                       cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa                                                                                                             >  2:415946/1‑150 (MQ=255)
                                        gATGAATTATCGCTCCCTTGACTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc                                                                                                            <  1:2450344/150‑1 (MQ=255)
                                         aTGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACg                                                                                                           <  2:4463146/150‑1 (MQ=255)
                                           gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc                                                                                                         >  2:3592425/1‑150 (MQ=255)
                                               tATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTa                                                                                                     >  1:4294129/1‑150 (MQ=255)
                                                aTCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAt                                                                                                    >  1:3097994/1‑150 (MQ=255)
                                                 tCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtt                                                                                                   <  2:2915561/150‑1 (MQ=255)
                                                    cTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgc                                                                                                <  2:3714381/150‑1 (MQ=255)
                                                      cccTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgctt                                                                                              <  2:1838304/150‑1 (MQ=255)
                                                       ccTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgcttt                                                                                             <  2:686667/150‑1 (MQ=255)
                                                          tGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGcc                                                                                          >  1:552325/1‑150 (MQ=255)
                                                          tGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGcc                                                                                          <  2:4302976/150‑1 (MQ=255)
                                                          tGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGcc                                                                                          >  2:552887/1‑150 (MQ=255)
                                                             cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                       <  1:1468643/150‑1 (MQ=255)
                                                              ttATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                      >  1:1408968/1‑150 (MQ=255)
                                                               tataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTAc                                                                                     >  2:3959736/1‑150 (MQ=255)
                                                               tataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTAc                                                                                     >  1:3059533/1‑150 (MQ=255)
                                                                ataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACg                                                                                    >  1:1908884/1‑150 (MQ=255)
                                                                ataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACg                                                                                    >  2:1476963/1‑150 (MQ=255)
                                                                 taACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGa                                                                                   >  2:2433996/1‑150 (MQ=255)
                                                                  aaCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGAt                                                                                  >  2:3395456/1‑150 (MQ=255)
                                                                        tCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGaa                                                                            >  1:3670501/1‑150 (MQ=255)
                                                                        tCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGaa                                                                            >  1:1386300/1‑150 (MQ=255)
                                                                         cAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAAc                                                                           >  2:3704446/1‑150 (MQ=255)
                                                                         cAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAAc                                                                           >  2:2459171/1‑150 (MQ=255)
                                                                         cAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAAc                                                                           >  2:1397706/1‑150 (MQ=255)
                                                                         cAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAAc                                                                           >  2:1413439/1‑150 (MQ=255)
                                                                          aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                          >  2:4132612/1‑150 (MQ=255)
                                                                          aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                          <  2:2111195/150‑1 (MQ=255)
                                                                           gACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAg                                                                         <  2:4202381/150‑1 (MQ=255)
                                                                            aCATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGc                                                                        <  1:1013732/150‑1 (MQ=255)
                                                                            aCATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGc                                                                        <  1:415946/150‑1 (MQ=255)
                                                                              aTGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCAt                                                                      >  1:3518554/1‑150 (MQ=255)
                                                                               tGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATa                                                                     >  2:4171960/1‑150 (MQ=255)
                                                                                ggTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATag                                                                    <  1:954430/150‑1 (MQ=255)
                                                                                ggTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATag                                                                    <  2:2543494/150‑1 (MQ=255)
                                                                                 gTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATaga                                                                   <  2:1683142/150‑1 (MQ=255)
                                                                                  ttAGTGTTTCGTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTCCCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATGTGCTTTGCCTTTACGATCATGAACAGCATagag                                                                  >  2:2801768/1‑150 (MQ=255)
                                                                                    aGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGaa                                                                >  1:3420092/1‑150 (MQ=255)
                                                                                     gtgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAAt                                                               <  2:3988502/150‑1 (MQ=255)
                                                                                     gtgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAAt                                                               <  2:565564/150‑1 (MQ=255)
                                                                                        tttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaaa                                                            <  1:892234/150‑1 (MQ=255)
                                                                                                cattcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcaggg                                                    <  1:1278755/150‑1 (MQ=255)
                                                                                                cattcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcaggg                                                    <  1:4336469/150‑1 (MQ=255)
                                                                                                 attcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcagggc                                                   >  1:49015/1‑150 (MQ=255)
                                                                                                 attcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcagggc                                                   >  2:1531385/1‑150 (MQ=255)
                                                                                                  ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAgggcagggca                                                  <  2:2854705/150‑1 (MQ=255)
                                                                                                    catcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAg                                                >  1:2856344/1‑150 (MQ=255)
                                                                                                     atcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                               <  1:995126/150‑1 (MQ=255)
                                                                                                      tcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTa                                              <  2:191326/150‑1 (MQ=255)
                                                                                                      tcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTa                                              >  2:1726840/1‑150 (MQ=255)
                                                                                                       cacaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTCTTTGCTTTGCCTTTACGATCATGAACAGCATAGAGACTAAAGGGCAGGGCAAGTac                                             >  2:4283765/1‑150 (MQ=255)
                                                                                                        acaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTacc                                            >  2:3848601/1‑150 (MQ=255)
                                                                                                           atatAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTaccacc                                         >  2:2817854/1‑150 (MQ=255)
                                                                                                            tataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTAccacca                                        >  2:3157612/1‑150 (MQ=255)
                                                                                                            tataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTAccacca                                        <  1:3200844/150‑1 (MQ=255)
                                                                                                             ataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                       >  1:3946408/1‑150 (MQ=255)
                                                                                                                gtgtGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGaa                                    >  2:2333513/1‑150 (MQ=255)
                                                                                                                  gtggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAAc                                  >  1:3802512/1‑150 (MQ=255)
                                                                                                                   tggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACt                                 <  2:852762/150‑1 (MQ=255)
                                                                                                                     gtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaa                               >  2:1512248/1‑150 (MQ=255)
                                                                                                                     gtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATCGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaa                               <  2:288443/150‑1 (MQ=255)
                                                                                                                        aaCGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaacaa                            <  1:2938557/150‑1 (MQ=255)
                                                                                                                          cGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAAcc                          <  1:3255247/150‑1 (MQ=255)
                                                                                                                           gTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                         <  1:1208574/150‑1 (MQ=255)
                                                                                                                            tGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAg                        >  1:789397/1‑150 (MQ=255)
                                                                                                                              cACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTa                      >  1:1417498/1‑150 (MQ=255)
                                                                                                                               aCGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTaa                     <  1:2218637/150‑1 (MQ=255)
                                                                                                                                   gggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTACCAACCAGTACTTCa                 >  1:4408594/1‑150 (MQ=255)
                                                                                                                                   gggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa                 <  2:3171986/150‑1 (MQ=255)
                                                                                                                                     gTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAac               >  1:4060816/1‑150 (MQ=255)
                                                                                                                                      tGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAaca              >  1:1321276/1‑150 (MQ=255)
                                                                                                                                        cAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacata            >  1:1070057/1‑150 (MQ=255)
                                                                                                                                         aGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatac           <  2:3295457/150‑1 (MQ=255)
                                                                                                                                         aGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatac           <  2:1172671/150‑1 (MQ=255)
                                                                                                                                           gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatacat         >  1:675645/1‑150 (MQ=255)
                                                                                                                                            aagaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAAcatacata        >  1:266110/1‑150 (MQ=255)
                                                                                                                                             agaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACatat       <  1:2698720/150‑1 (MQ=255)
                                                                                                                                               aagTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCg     <  1:194054/150‑1 (MQ=255)
                                                                                                                                                agTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCgg    <  2:4019725/150‑1 (MQ=255)
                                                                                                                                                agTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCgg    >  2:2415864/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGGTA  >  CP009273/1563129‑1563424

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