Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,678,647 C→A 5.4% V74F (GTT→TTT)  glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,678,6470CA 5.4% 1183.0 / 26.1 295V74F (GTT→TTT) glyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (138/141);  new base A (3/13);  total (141/154)
Fisher's exact test for biased strand distribution p-value = 1.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 299 total aligned reads displayed.

ATGCGCCAGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTA  >  CP009273/2678498‑2678796
                                                                                                                                                     |                                                                                                                                                     
   cgcCAGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGat                                                                                                                                                    >  1:4518563/1‑150 (MQ=255)
     ccAGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGatat                                                                                                                                                  <  1:1030640/150‑1 (MQ=255)
      cAGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATc                                                                                                                                                 <  1:2468885/150‑1 (MQ=255)
       aGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCa                                                                                                                                                <  2:1897344/150‑1 (MQ=255)
       aGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAAAGATATCa                                                                                                                                                <  2:482063/150‑1 (MQ=255)
        ggTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCaa                                                                                                                                               >  2:4701310/1‑150 (MQ=255)
          ttCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAAAGATATCAACa                                                                                                                                             <  2:3187797/150‑1 (MQ=255)
          ttCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTGTCGCACGATCGAGCGCCAGTTGTTCAACGATATCACCa                                                                                                                                             >  2:2604027/1‑150 (MQ=255)
           tCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACAt                                                                                                                                            <  1:225971/150‑1 (MQ=255)
            cATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATa                                                                                                                                           >  1:5053448/1‑150 (MQ=255)
               aCCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTc                                                                                                                                        <  1:2097135/150‑1 (MQ=255)
               aCCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTc                                                                                                                                        >  2:1529828/1‑150 (MQ=255)
               aCCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTc                                                                                                                                        <  1:1373142/150‑1 (MQ=255)
                 ccAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGc                                                                                                                                      >  2:4395339/1‑150 (MQ=255)
                  cAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCa                                                                                                                                     >  1:609902/1‑150 (MQ=255)
                  cAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCa                                                                                                                                     <  2:875791/150‑1 (MQ=255)
                          tATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgt                                                                                                                             <  2:4259673/150‑1 (MQ=255)
                           aTCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgta                                                                                                                            <  1:3070642/150‑1 (MQ=255)
                           aTCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgta                                                                                                                            <  1:3825991/150‑1 (MQ=255)
                            tCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCCGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATCCTCGCAACCGCCgtag                                                                                                                           >  2:5025004/1‑150 (MQ=255)
                            tCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgtag                                                                                                                           <  2:4378101/150‑1 (MQ=255)
                            tCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgtag                                                                                                                           >  1:144037/1‑150 (MQ=255)
                             cACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgtagt                                                                                                                          <  1:236806/150‑1 (MQ=255)
                              aCCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgtagta                                                                                                                         <  1:4976959/150‑1 (MQ=255)
                               ccTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGTAGCCGGAGGGCGGCTGGACGTTAGCGTAGTCAGCGCCGGACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAAAATACTCGCAACAGCCGtagtag                                                                                                                        >  2:1948928/1‑150 (MQ=255)
                               ccTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGtagtag                                                                                                                        <  2:3272715/150‑1 (MQ=255)
                                cTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGc                                                                                                                       >  1:4285671/1‑150 (MQ=255)
                                cTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGc                                                                                                                       <  2:3877601/150‑1 (MQ=255)
                                  ggTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGt                                                                                                                     >  1:2689185/1‑150 (MQ=255)
                                       cagcagCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGcc                                                                                                                >  2:2445617/1‑150 (MQ=255)
                                        agcagcGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGccc                                                                                                               >  1:3460645/1‑150 (MQ=255)
                                        agcagcGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGccc                                                                                                               >  2:125631/1‑150 (MQ=255)
                                          cagcGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCgg                                                                                                             <  2:700515/150‑1 (MQ=255)
                                          cagcGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAAAGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCgg                                                                                                             >  1:2059835/1‑150 (MQ=255)
                                              gcgGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATaa                                                                                                         <  2:3521920/150‑1 (MQ=255)
                                                   gtAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTc                                                                                                    >  1:4815833/1‑150 (MQ=255)
                                                    tAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCa                                                                                                   >  2:4673833/1‑150 (MQ=255)
                                                    tAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCa                                                                                                   <  2:1775105/150‑1 (MQ=255)
                                                     aGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAg                                                                                                  <  1:1484768/150‑1 (MQ=255)
                                                       aCCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCa                                                                                                <  1:3314738/150‑1 (MQ=255)
                                                         cGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCata                                                                                              <  2:450213/150‑1 (MQ=255)
                                                          gCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCatat                                                                                             <  1:324155/150‑1 (MQ=255)
                                                            aaaGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAAAGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATAttt                                                                                           <  1:3495397/150‑1 (MQ=255)
                                                               gTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATttgtt                                                                                        >  2:586200/1‑150 (MQ=255)
                                                                ttAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATttgttg                                                                                       <  2:4543714/150‑1 (MQ=255)
                                                                   gCCTGGGAGACGGAGTGCTGCTGGACGTTAGCGTTGTCAGCTCCGAACAGTTCTTTCTCAAGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCATGCGCCTTAGTAGCTTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc                                                                                    <  2:523933/150‑1 (MQ=255)
                                                                    ccTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCa                                                                                   >  1:448873/1‑150 (MQ=255)
                                                                    ccTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCa                                                                                   <  1:3159748/150‑1 (MQ=255)
                                                                        ggAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTg                                                                               <  1:1529828/150‑1 (MQ=255)
                                                                           gCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTgaga                                                                            >  2:5126729/1‑150 (MQ=255)
                                                                             cGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAAc                                                                          >  2:3802378/1‑150 (MQ=255)
                                                                             cGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAAc                                                                          >  1:2804847/1‑150 (MQ=255)
                                                                               gAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAAccc                                                                        <  1:1661565/150‑1 (MQ=255)
                                                                                aGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCt                                                                       <  2:2105951/150‑1 (MQ=255)
                                                                                    cGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTgcgc                                                                   >  1:479759/1‑150 (MQ=255)
                                                                                        tGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGTGCCTGc                                                               <  1:1126009/150‑1 (MQ=255)
                                                                                           aCGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCAtt                                                            <  2:1142699/150‑1 (MQ=255)
                                                                                             gTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAc                                                          <  2:3085973/150‑1 (MQ=255)
                                                                                             gTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAc                                                          >  2:4464887/1‑150 (MQ=255)
                                                                                              ttAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcg                                                         <  2:2192236/150‑1 (MQ=255)
                                                                                               tAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgc                                                        <  2:781707/150‑1 (MQ=255)
                                                                                               tAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCAGATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgc                                                        >  1:1896598/1‑150 (MQ=255)
                                                                                               tAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTAAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgc                                                        >  1:4909422/1‑150 (MQ=255)
                                                                                                aGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgcg                                                       <  1:3371591/150‑1 (MQ=255)
                                                                                                aGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgcg                                                       <  1:3371451/150‑1 (MQ=255)
                                                                                                  cGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgcgcg                                                     >  2:2552148/1‑150 (MQ=255)
                                                                                                   gTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCgg                                                    <  1:2013523/150‑1 (MQ=255)
                                                                                                   gTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCgg                                                    >  1:1285803/1‑150 (MQ=255)
                                                                                                    tAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCggg                                                   >  2:1339694/1‑150 (MQ=255)
                                                                                                     aGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGc                                                  >  2:4145822/1‑150 (MQ=255)
                                                                                                      gTCAGCGCCGAACTGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCTTATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCt                                                 <  1:3667931/150‑1 (MQ=255)
                                                                                                       tCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTg                                                >  1:3985766/1‑150 (MQ=255)
                                                                                                          gcgcCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGgtg                                             <  1:2418023/150‑1 (MQ=255)
                                                                                                          gcgcCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAAAGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGgtg                                             <  1:4141632/150‑1 (MQ=255)
                                                                                                             ccGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                          <  1:4725163/150‑1 (MQ=255)
                                                                                                             ccGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                          <  1:4725234/150‑1 (MQ=255)
                                                                                                             ccGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                          >  1:1019011/1‑150 (MQ=255)
                                                                                                              cGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGt                                         >  1:516347/1‑150 (MQ=255)
                                                                                                              cGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGt                                         >  1:1621049/1‑150 (MQ=255)
                                                                                                                         ttCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCg                              <  1:2863647/150‑1 (MQ=255)
                                                                                                                           cGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGAt                            >  1:1246116/1‑150 (MQ=255)
                                                                                                                           cGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGAt                            >  1:916329/1‑150 (MQ=255)
                                                                                                                             cACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCa                          >  1:2996590/1‑150 (MQ=255)
                                                                                                                             cACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCa                          <  2:2293749/150‑1 (MQ=255)
                                                                                                                              aCGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAg                         >  2:1851248/1‑150 (MQ=255)
                                                                                                                               cgatcgatCGCCAGTTGTTCAAAGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGt                        <  2:3362729/150‑1 (MQ=255)
                                                                                                                                gatcgatcGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGtt                       >  1:2665897/1‑150 (MQ=255)
                                                                                                                                    gatcgCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAt                   <  1:4423143/150‑1 (MQ=255)
                                                                                                                                     atcgCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAtg                  >  1:2524948/1‑150 (MQ=255)
                                                                                                                                      tcgCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAtgt                 <  2:962820/150‑1 (MQ=255)
                                                                                                                                        gCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGc               >  2:1653668/1‑150 (MQ=255)
                                                                                                                                         ccAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGct              >  1:3635155/1‑150 (MQ=255)
                                                                                                                                           atttgttCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGctct            >  2:2930636/3‑150 (MQ=255)
                                                                                                                                           aGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGctct            >  2:2699724/1‑150 (MQ=255)
                                                                                                                                           aGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGctct            <  2:4173262/150‑1 (MQ=255)
                                                                                                                                           aGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGAGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGctct            >  2:4039513/1‑150 (MQ=255)
                                                                                                                                           aGTTGTTCAAAGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGctct            >  1:73546/1‑150 (MQ=255)
                                                                                                                                            gttgttCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCtt           <  2:4650584/150‑1 (MQ=255)
                                                                                                                                            gttgttCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTGTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCtt           <  1:3931534/150‑1 (MQ=255)
                                                                                                                                              tgttCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGAGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTcc         >  1:901970/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCGCCAGGTTCATACCCAGAACGGTATCACCTGGTTCCAGCAGCGCGGTGTAGACCGCAAAGTTAGCCTGGGAGCCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTA  >  CP009273/2678498‑2678796

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: