Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,510,923 G→T 24.9% T363T (ACG→ACT patD → gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,510,9230GT24.9% 778.1 / 194.0 337T363T (ACG→ACTpatDgamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (138/115);  new base T (45/39);  total (183/154)
Fisher's exact test for biased strand distribution p-value = 9.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 352 total aligned reads displayed.

GACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAA  >  CP009273/1510774‑1511072
                                                                                                                                                     |                                                                                                                                                     
   gAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgctg                                                                                                                                                    >  1:153864/1‑150 (MQ=255)
      tCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCACCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGctgctg                                                                                                                                                 <  2:2272461/150‑1 (MQ=255)
        tACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggc                                                                                                                                               <  2:3574454/150‑1 (MQ=255)
         aCGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTgctggct                                                                                                                                              <  2:5038906/150‑1 (MQ=255)
          cGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggctg                                                                                                                                             <  2:790915/150‑1 (MQ=255)
           ggAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGctggctgg                                                                                                                                            <  2:3674980/150‑1 (MQ=255)
                 tGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCAtt                                                                                                                                      >  2:2677183/1‑150 (MQ=255)
                  ggACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTa                                                                                                                                     <  1:5195189/150‑1 (MQ=255)
                   gACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc                                                                                                                                    >  2:4895259/1‑150 (MQ=255)
                    aCCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACa                                                                                                                                   <  1:1785530/150‑1 (MQ=255)
                       tttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCTGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGa                                                                                                                                <  1:5136040/150‑1 (MQ=255)
                        ttAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGAc                                                                                                                               <  1:2003590/150‑1 (MQ=255)
                         tAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACg                                                                                                                              >  1:5052327/1‑150 (MQ=255)
                          aaGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGa                                                                                                                             >  2:4045164/1‑150 (MQ=255)
                           tgCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGAt                                                                                                                            <  2:4494891/149‑1 (MQ=255)
                            gCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATg                                                                                                                           <  1:1902381/150‑1 (MQ=255)
                             cTCGCTGGCGCATCTCGAACGCGTCGGCAAGGATGTATAAGTGGCGAAAGCGACAGGGCACTTCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGc                                                                                                                          <  1:643200/150‑1 (MQ=255)
                                      gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                                                                                                 <  2:4655634/150‑1 (MQ=255)
                                      gcATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCa                                                                                                                 >  1:2171963/1‑150 (MQ=255)
                                        aTCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaa                                                                                                               >  1:4293927/1‑150 (MQ=255)
                                          ctcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCaaaaa                                                                                                             <  1:4405486/150‑1 (MQ=255)
                                           tcGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAag                                                                                                            <  2:2606924/150‑1 (MQ=255)
                                            cGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAaga                                                                                                           <  2:236370/150‑1 (MQ=255)
                                             gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                                                                                          >  2:179717/1‑150 (MQ=255)
                                             gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                                                                                          >  1:4815746/1‑150 (MQ=255)
                                             gAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAagag                                                                                                          >  1:2899607/1‑150 (MQ=255)
                                               aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGTCGATGCCATCGTGCAAAAAGAGGt                                                                                                        <  2:2944957/150‑1 (MQ=255)
                                               aCGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGt                                                                                                        <  2:5185525/150‑1 (MQ=255)
                                                  cgTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAtt                                                                                                     >  2:3499034/1‑150 (MQ=255)
                                                   gTCGGCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                                                                                                    <  1:4591275/150‑1 (MQ=255)
                                                   gTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTAttt                                                                                                    <  2:3400988/150‑1 (MQ=255)
                                                      ggCAAGGCAGTAGAAGAGGCGTAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                                                                                                 >  1:5022856/1‑150 (MQ=255)
                                                      ggCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGt                                                                                                 >  1:1844929/1‑150 (MQ=255)
                                                       gCAAGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTc                                                                                                <  2:649937/150‑1 (MQ=255)
                                                         aaGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                                                                                              >  2:175036/1‑150 (MQ=255)
                                                         aaGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                                                                                              >  1:1179962/1‑150 (MQ=255)
                                                         aaGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                                                                                              >  2:621693/1‑150 (MQ=255)
                                                         aaGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCa                                                                                              >  1:2433071/1‑150 (MQ=255)
                                                          aGGCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCag                                                                                             >  2:1573033/1‑150 (MQ=255)
                                                           ggCAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagt                                                                                            <  1:588239/150‑1 (MQ=255)
                                                             cAGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAAAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtag                                                                                          <  1:3779425/150‑1 (MQ=255)
                                                             cAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtag                                                                                          <  2:4369881/150‑1 (MQ=255)
                                                              aGTAGAAGAGGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCagtagt                                                                                         <  1:1618910/150‑1 (MQ=255)
                                                                   aagaGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagtg                                                                                    >  2:1380081/1‑150 (MQ=255)
                                                                   aagaGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTagtgagtg                                                                                    <  2:4273520/150‑1 (MQ=255)
                                                                      agGCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTa                                                                                 <  1:3960814/150‑1 (MQ=255)
                                                                       ggCGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTAc                                                                                <  2:3713362/150‑1 (MQ=255)
                                                                         cGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGATGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGc                                                                              <  2:868194/150‑1 (MQ=255)
                                                                         cGAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGc                                                                              >  2:479086/1‑150 (MQ=255)
                                                                          gAAAGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCAAGTAGTGAGTGTTACGcc                                                                             >  2:1291835/1‑150 (MQ=255)
                                                                           aaaGCTACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGccc                                                                            >  2:1460946/1‑150 (MQ=255)
                                                                           aaaGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGccc                                                                            >  2:2830688/1‑150 (MQ=255)
                                                                              gctaCAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTc                                                                         >  2:4716308/4‑150 (MQ=255)
                                                                                gACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGa                                                                       >  2:3111155/1‑150 (MQ=255)
                                                                                 aCAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGAc                                                                      <  1:1000326/150‑1 (MQ=255)
                                                                                    gggCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAAc                                                                   >  1:1422483/1‑150 (MQ=255)
                                                                                     ggCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATGACAGGACGATGCCATCGAGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACg                                                                  >  1:134311/1‑150 (MQ=255)
                                                                                       cacaTCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACgaa                                                                <  1:1287932/150‑1 (MQ=255)
                                                                                        acaTCAAAGTCAACACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCAGGCTGGCGCATTACAGGACGATGCCAGCGTGCACAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTGCGACAACgaag                                                               >  2:2032413/1‑150 (MQ=255)
                                                                                           tcaccgTGATCACTGGCGGTGAACAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAAc                                                            <  1:4079586/145‑1 (MQ=255)
                                                                                             aaaGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAg                                                          >  2:3644330/1‑150 (MQ=255)
                                                                                             aaaGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGGCTTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGACCAg                                                          >  2:1374634/1‑150 (MQ=255)
                                                                                               aGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAggt                                                        <  1:84267/150‑1 (MQ=255)
                                                                                               aGTGATCACTGGAGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACTCCCTACGACAACGAAGAACAggt                                                        >  2:4613507/1‑150 (MQ=255)
                                                                                                 tGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAggtgg                                                      <  2:1626159/150‑1 (MQ=255)
                                                                                                  gATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAggtggt                                                     >  1:3163814/1‑150 (MQ=255)
                                                                                                  gATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAggtggt                                                     <  2:3259127/150‑1 (MQ=255)
                                                                                                   aTCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGGATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGgtggtg                                                    >  1:3448610/1‑150 (MQ=255)
                                                                                                      aCTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAAc                                                 <  2:4963983/150‑1 (MQ=255)
                                                                                                      aCTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAAc                                                 <  1:1060381/150‑1 (MQ=255)
                                                                                                       cTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACt                                                <  1:2625669/150‑1 (MQ=255)
                                                                                                         ggCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTgg                                              >  2:1123322/1‑150 (MQ=255)
                                                                                                            ggTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCg                                           <  1:2023540/150‑1 (MQ=255)
                                                                                                              tGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGaa                                         >  1:230448/1‑150 (MQ=255)
                                                                                                               gAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGact                                        >  2:1233043/1‑148 (MQ=255)
                                                                                                                aaaaGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATg                                       <  1:1578244/150‑1 (MQ=255)
                                                                                                                 aaaGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGa                                      <  2:3921943/150‑1 (MQ=255)
                                                                                                                  aaGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGAc                                     >  2:938122/1‑150 (MQ=255)
                                                                                                                     cgcAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGc                                  >  1:3311565/1‑150 (MQ=255)
                                                                                                                      gcAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGcc                                 <  2:2316926/150‑1 (MQ=255)
                                                                                                                      gcAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGcc                                 <  1:858730/150‑1 (MQ=255)
                                                                                                                      gcAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGcc                                 <  1:4955035/150‑1 (MQ=255)
                                                                                                                       cAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCa                                >  2:1725788/1‑150 (MQ=255)
                                                                                                                        aaGGGTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAg                               >  1:3398120/1‑150 (MQ=255)
                                                                                                                         aGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGt                              >  2:2180395/1‑150 (MQ=255)
                                                                                                                          gggTAATGGCTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTa                             >  1:3543177/1‑150 (MQ=255)
                                                                                                                            gTAATGGCTTTTTCTTTGCGCAGACGCTGCTGGCTGGAGCATTACAGGACGATGGCATCGTGCATTAAGAGGTTTTTGGTCCAGTAGTGTGTGTTTCGGCCTTCGTCAACGTAGAACTGGTGGTGAACTGGGCGAATGACAGCCAGTACg                           <  1:2995184/150‑1 (MQ=255)
                                                                                                                            gTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACg                           >  1:2098678/1‑150 (MQ=255)
                                                                                                                             tAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACgg                          >  1:2228372/1‑150 (MQ=255)
                                                                                                                               aTGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGAc                        >  1:2686242/1‑150 (MQ=255)
                                                                                                                                tGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACt                       >  2:1515719/1‑150 (MQ=255)
                                                                                                                                  gCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTg                     >  1:44747/1‑150 (MQ=255)
                                                                                                                                   cTATTACTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGc                    >  1:1876693/1‑150 (MQ=255)
                                                                                                                                        aCTATGCGCCGACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCtt               <  2:1054733/150‑1 (MQ=255)
                                                                                                                                          tATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTCGTGAGTGTTACGCCCTTCGACAACGAAGACCAGGTGGTGAACTGGGCGACTGACAGCCAGTACGGACTTGCATCTGCg             >  2:4175942/1‑150 (MQ=255)
                                                                                                                                          tATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCg             >  1:4750264/1‑150 (MQ=255)
                                                                                                                                                 cGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGGATGGa      >  2:513569/1‑150 (MQ=255)
                                                                                                                                                  gACTCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGAc     >  1:3270220/1‑150 (MQ=255)
                                                                                                                                                  aaCGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAATAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGAc     >  1:3668053/2‑150 (MQ=255)
                                                                                                                                                     gctgctGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGaa  >  1:733653/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTACAGGACGATGCCATCGTGCAAAAAGAGGTATTTGGTCCAGTAGTGAGTGTTACGCCCTTCGACAACGAAGAACAGGTGGTGAACTGGGCGAATGACAGCCAGTACGGACTTGCATCTTCGGTATGGACGAA  >  CP009273/1510774‑1511072

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: