Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,946,256 A→T 10.4% Q138L (CAG→CTG)  ilvE → branched‑chain amino‑acid aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,946,2560AT10.4% 675.2 / 55.4 230Q138L (CAG→CTG) ilvEbranched‑chain amino‑acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (102/104);  new base T (9/15);  total (111/119)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 235 total aligned reads displayed.

AAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGA  >  CP009273/3946109‑3946403
                                                                                                                                                   |                                                                                                                                                   
aaaaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCa                                                                                                                                                     <  2:2567509/148‑1 (MQ=255)
  aaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAgg                                                                                                                                                   <  2:2212356/148‑1 (MQ=255)
  aaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAgg                                                                                                                                                   <  2:1738215/148‑1 (MQ=255)
  aaCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAgg                                                                                                                                                   <  2:1738102/148‑1 (MQ=255)
   aCAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAggg                                                                                                                                                  <  1:90457/148‑1 (MQ=255)
     aaTCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGa                                                                                                                                                <  1:1899646/148‑1 (MQ=255)
      aTCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGAt                                                                                                                                               <  1:1788936/148‑1 (MQ=255)
        ctcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcg                                                                                                                                             <  1:1043405/148‑1 (MQ=255)
         tcACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcga                                                                                                                                            <  2:1907674/148‑1 (MQ=255)
          cACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgat                                                                                                                                           >  2:558920/1‑148 (MQ=255)
            ccAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatgc                                                                                                                                         >  2:503713/1‑148 (MQ=255)
             cAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatgcg                                                                                                                                        >  2:3108884/1‑148 (MQ=255)
                cgcCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatgcgatg                                                                                                                                     <  2:2850185/148‑1 (MQ=255)
                cgcCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATcgatgcgatg                                                                                                                                     <  2:2154203/148‑1 (MQ=255)
                 gcCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATgg                                                                                                                                    <  2:910150/148‑1 (MQ=255)
                   cTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGtt                                                                                                                                  >  2:584734/1‑148 (MQ=255)
                     atatCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtc                                                                                                                                >  1:3555139/1‑148 (MQ=255)
                       atCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcct                                                                                                                              <  1:820119/148‑1 (MQ=255)
                       atCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcct                                                                                                                              >  2:2934408/1‑148 (MQ=255)
                       atCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcct                                                                                                                              >  2:393203/1‑148 (MQ=255)
                       atCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcct                                                                                                                              >  2:1156754/1‑148 (MQ=255)
                        tccgtccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcctc                                                                                                                             <  1:709220/148‑1 (MQ=255)
                        tccgtccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcctc                                                                                                                             >  1:2875046/1‑148 (MQ=255)
                         ccgtccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTtcctcc                                                                                                                            <  1:2210707/148‑1 (MQ=255)
                          cgtccgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTcctcct                                                                                                                           >  1:818547/1‑148 (MQ=255)
                              cgCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGaa                                                                                                                       <  1:2119998/148‑1 (MQ=255)
                                 tGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACcg                                                                                                                    >  2:2074125/1‑148 (MQ=255)
                                 tGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACcg                                                                                                                    <  2:2598237/148‑1 (MQ=255)
                                 tGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACcg                                                                                                                    <  2:1229567/148‑1 (MQ=255)
                                 tGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACcg                                                                                                                    <  1:2029701/148‑1 (MQ=255)
                                 tGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCGTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACcg                                                                                                                    >  1:2967225/1‑148 (MQ=255)
                                 tGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGATTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACcg                                                                                                                    >  2:1171505/1‑148 (MQ=255)
                                   aTCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACcgcg                                                                                                                  <  2:2401351/148‑1 (MQ=255)
                                              tGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAac                                                                                                       >  2:2022270/1‑148 (MQ=255)
                                                 tGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACAcc                                                                                                    <  1:3047008/148‑1 (MQ=255)
                                                  gTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCa                                                                                                   >  2:2753972/1‑148 (MQ=255)
                                                   ttGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCAt                                                                                                  <  2:3128468/148‑1 (MQ=255)
                                                       cATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCg                                                                                              >  1:584098/1‑148 (MQ=255)
                                                        aTGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCTCCAAACACCATCCCGa                                                                                             >  1:348625/1‑148 (MQ=255)
                                                            gAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggc                                                                                         >  2:2602253/1‑148 (MQ=255)
                                                             aGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                                                        >  1:1032944/1‑148 (MQ=255)
                                                             aGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcg                                                                                        >  2:398044/1‑148 (MQ=255)
                                                              gTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGAcggcgg                                                                                       >  1:2601866/1‑148 (MQ=255)
                                                               tAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACggcggc                                                                                      >  2:2167445/1‑148 (MQ=255)
                                                                  aCCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaa                                                                                   <  1:1456587/148‑1 (MQ=255)
                                                                   cccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                                                  >  1:1210544/1‑148 (MQ=255)
                                                                   cccGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCaaaa                                                                                  <  2:1345547/148‑1 (MQ=255)
                                                                     cGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGc                                                                                >  2:1958653/1‑148 (MQ=255)
                                                                      gCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGcc                                                                               >  2:3303542/1‑148 (MQ=255)
                                                                       ccAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCg                                                                              <  1:2365083/148‑1 (MQ=255)
                                                                        cAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCgg                                                                             >  2:1426413/1‑148 (MQ=255)
                                                                          gCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCggtg                                                                           >  1:2857549/1‑148 (MQ=255)
                                                                             ggATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCACAAGCCGGTGGTa                                                                        >  2:61739/1‑148 (MQ=255)
                                                                              gATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGGTGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                                       <  1:1426247/148‑1 (MQ=255)
                                                                              gATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                                       >  2:1397244/1‑148 (MQ=255)
                                                                              gATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                                       <  1:439752/148‑1 (MQ=255)
                                                                              gATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                                       >  1:1816483/1‑148 (MQ=255)
                                                                              gATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTaa                                                                       >  1:1212439/1‑148 (MQ=255)
                                                                                 aCTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTa                                                                    >  2:3206439/1‑148 (MQ=255)
                                                                                   tCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTAcc                                                                  >  2:1873282/1‑148 (MQ=255)
                                                                                    cAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                                                 <  1:3487319/148‑1 (MQ=255)
                                                                                    cAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACct                                                                 <  2:275616/148‑1 (MQ=255)
                                                                                      aCCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctct                                                               >  1:389882/1‑148 (MQ=255)
                                                                                       ccGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctc                                                              <  2:2270040/148‑1 (MQ=255)
                                                                                        cGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctct                                                             <  1:531143/148‑1 (MQ=255)
                                                                                        cGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACctctct                                                             >  2:832019/1‑148 (MQ=255)
                                                                                          aCGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTc                                                           >  2:1495462/1‑148 (MQ=255)
                                                                                            gTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                                         <  1:884027/148‑1 (MQ=255)
                                                                                            gTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTccc                                                         <  1:2668904/148‑1 (MQ=255)
                                                                                             tGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCct                                                        >  1:2533891/1‑148 (MQ=255)
                                                                                              gATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctg                                                       >  2:2026086/1‑148 (MQ=255)
                                                                                                ttATCGCTGCGTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgct                                                     <  1:1330516/148‑1 (MQ=255)
                                                                                                 tATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCctgctg                                                    <  2:3513865/148‑1 (MQ=255)
                                                                                                    cGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCtggtg                                                 >  1:1948885/1‑148 (MQ=255)
                                                                                                       tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                                              <  2:3076824/148‑1 (MQ=255)
                                                                                                       tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                                              <  2:3082373/148‑1 (MQ=255)
                                                                                                       tgctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGt                                              >  1:3326526/1‑148 (MQ=255)
                                                                                                         ctTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAg                                            <  2:452087/148‑1 (MQ=255)
                                                                                                          tttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGc                                           <  2:2100817/148‑1 (MQ=255)
                                                                                                           ttCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCg                                          >  2:687939/1‑148 (MQ=255)
                                                                                                              ccGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAg                                       <  2:3294541/148‑1 (MQ=255)
                                                                                                                  ggggAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAgcgcg                                   <  1:119760/148‑1 (MQ=255)
                                                                                                                                ggCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTtctc                     >  2:368249/1‑148 (MQ=255)
                                                                                                                                   gcAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAgg                  <  2:2460181/148‑1 (MQ=255)
                                                                                                                                    cAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGa                 <  2:2985033/148‑1 (MQ=255)
                                                                                                                                       aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAgg              >  2:2917271/1‑148 (MQ=255)
                                                                                                                                       aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAgg              <  1:2816297/148‑1 (MQ=255)
                                                                                                                                       aaGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAgg              >  2:1138834/1‑148 (MQ=255)
                                                                                                                                         gcgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTa            <  1:1089713/148‑1 (MQ=255)
                                                                                                                                         gcgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTa            <  2:1039727/148‑1 (MQ=255)
                                                                                                                                          cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAATCACCATCCCGACGGCGGCAAACGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTAt           <  2:227389/148‑1 (MQ=255)
                                                                                                                                          cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTAt           >  1:792654/1‑148 (MQ=255)
                                                                                                                                          cgcTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTAt           <  1:2084635/148‑1 (MQ=255)
                                                                                                                                             tGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATcgc        <  2:1676776/148‑1 (MQ=255)
                                                                                                                                              ggAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATcgcg       <  1:2385855/148‑1 (MQ=255)
                                                                                                                                               gAGCTGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCgcgc      >  2:249496/1‑148 (MQ=255)
                                                                                                                                                  cAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTgg   >  1:3357901/1‑148 (MQ=255)
                                                                                                                                                  cAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTgg   <  1:3243759/148‑1 (MQ=255)
                                                                                                                                                  cAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTgg   >  1:493108/1‑148 (MQ=255)
                                                                                                                                                   aGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGa  <  1:2167445/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAACAATCTCACCAGCGCCTATATCCGTCCGCTGATCTTCGTCGGTGATGTTGGCATGGGAGTAAACCCGCCAGCGGGATACTCAACCGACGTGATTATCGCTGCTTTCCCGTGGGGAGCGTATCTGGGCGCAGAAGCGCTGGAGCAGGGGATCGATGCGATGGTTTCCTCCTGGAACCGCGCAGCACCAAACACCATCCCGACGGCGGCAAAAGCCGGTGGTAACTACCTCTCTTCCCTGCTGGTGGGTAGCGAAGCGCGCCGCCACGGTTATCAGGAAGGTATCGCGCTGGA  >  CP009273/3946109‑3946403

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: