Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,937,613 G→T 12.2% T306T (ACG→ACT csdA → cysteine sulfinate desulfinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,937,6130GT12.2% 951.1 / 71.1 270T306T (ACG→ACTcsdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base G (113/124);  new base T (9/24);  total (122/148)
Fisher's exact test for biased strand distribution p-value = 3.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 273 total aligned reads displayed.

TTTGACGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGCGGGCCGGGCA  >  CP009273/2937464‑2937762
                                                                                                                                                     |                                                                                                                                                     
    aCGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTGAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTgg                                                                                                                                                   >  1:3313266/1‑150 (MQ=255)
     cGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGc                                                                                                                                                  <  1:1928949/150‑1 (MQ=255)
      ggCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCg                                                                                                                                                 <  1:494206/150‑1 (MQ=255)
       gCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCgg                                                                                                                                                <  1:4006118/150‑1 (MQ=255)
       gCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCgg                                                                                                                                                >  2:2481378/1‑150 (MQ=255)
        cTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGa                                                                                                                                               >  2:3218413/1‑150 (MQ=255)
          tCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAg                                                                                                                                             >  2:1058003/1‑150 (MQ=255)
          tCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAg                                                                                                                                             >  2:1614910/1‑150 (MQ=255)
            aCGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGAt                                                                                                                                           >  2:3425581/1‑150 (MQ=255)
             cGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATg                                                                                                                                          <  2:2242698/150‑1 (MQ=255)
              gACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATgc                                                                                                                                         >  1:1134642/1‑150 (MQ=255)
                 tCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCt                                                                                                                                      <  2:3607711/150‑1 (MQ=255)
                     tCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCg                                                                                                                                  >  2:2494578/1‑150 (MQ=255)
                       tGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGaa                                                                                                                                >  1:3524195/1‑150 (MQ=255)
                        gcgcCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGaaa                                                                                                                               <  2:3023634/150‑1 (MQ=255)
                        gcgcCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGaaa                                                                                                                               >  2:4172061/1‑150 (MQ=255)
                              tGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTccc                                                                                                                         <  1:3374949/150‑1 (MQ=255)
                                gAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCgg                                                                                                                       <  1:303339/150‑1 (MQ=255)
                                  aaCTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCt                                                                                                                     <  2:2152428/150‑1 (MQ=255)
                                  aaCTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCt                                                                                                                     <  1:3800142/150‑1 (MQ=255)
                                  aaCTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCt                                                                                                                     <  1:3147587/150‑1 (MQ=255)
                                   aCTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCtt                                                                                                                    >  2:2818718/1‑150 (MQ=255)
                                    cTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCttt                                                                                                                   <  1:1479036/150‑1 (MQ=255)
                                    cTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCttt                                                                                                                   <  2:3971985/150‑1 (MQ=255)
                                     tGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTc                                                                                                                  >  2:2725599/1‑150 (MQ=255)
                                     tGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTc                                                                                                                  <  1:4160682/150‑1 (MQ=255)
                                      ggAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCg                                                                                                                 <  1:2170711/150‑1 (MQ=255)
                                       gAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGt                                                                                                                <  1:535478/150‑1 (MQ=255)
                                       gAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGt                                                                                                                >  1:1763498/1‑150 (MQ=255)
                                         aGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTc                                                                                                              >  2:1008747/1‑150 (MQ=255)
                                         aGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTc                                                                                                              >  1:1063285/1‑150 (MQ=255)
                                           cTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCAt                                                                                                            <  1:960158/150‑1 (MQ=255)
                                            tGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCAtt                                                                                                           >  1:3407777/1‑150 (MQ=255)
                                              gAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTcc                                                                                                         <  2:3715762/150‑1 (MQ=255)
                                               aaCGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCg                                                                                                        <  2:2601335/150‑1 (MQ=255)
                                                   ccAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGc                                                                                                    >  1:2027752/1‑150 (MQ=255)
                                                   ccAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGc                                                                                                    >  2:1728312/1‑150 (MQ=255)
                                                            gCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTcc                                                                                           <  2:4141908/150‑1 (MQ=255)
                                                             cTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCa                                                                                          >  2:4340968/1‑150 (MQ=255)
                                                             cTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCa                                                                                          >  1:1106769/1‑150 (MQ=255)
                                                             cTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCa                                                                                          <  2:3894776/150‑1 (MQ=255)
                                                               ggTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGc                                                                                        >  1:2771659/1‑150 (MQ=255)
                                                                 tgtCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCct                                                                                      <  1:4371305/150‑1 (MQ=255)
                                                                 tgtCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCct                                                                                      <  1:1605992/150‑1 (MQ=255)
                                                                   tCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCctgc                                                                                    >  2:4064453/1‑150 (MQ=255)
                                                                    cATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCctgct                                                                                   <  2:181614/150‑1 (MQ=255)
                                                                    cATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCGGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCctgct                                                                                   <  1:4319196/150‑1 (MQ=255)
                                                                     aTAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCctgctg                                                                                  <  2:615772/150‑1 (MQ=255)
                                                                     aTAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCctgctg                                                                                  >  1:697801/1‑150 (MQ=255)
                                                                       aGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGc                                                                                >  2:2524617/1‑150 (MQ=255)
                                                                        ggATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGcc                                                                               >  1:289308/1‑150 (MQ=255)
                                                                        ggATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGcc                                                                               <  1:1677954/150‑1 (MQ=255)
                                                                         gATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCt                                                                              <  2:2982395/150‑1 (MQ=255)
                                                                          aTTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCtt                                                                             >  1:882599/1‑150 (MQ=255)
                                                                          aTTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCtt                                                                             <  1:2584508/150‑1 (MQ=255)
                                                                          aTTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCtt                                                                             >  1:2641831/1‑150 (MQ=255)
                                                                           ttAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCttt                                                                            <  2:3647960/150‑1 (MQ=255)
                                                                            tAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTg                                                                           >  2:1919887/1‑150 (MQ=255)
                                                                            tAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTg                                                                           >  2:4246880/1‑150 (MQ=255)
                                                                              aGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGAt                                                                         >  1:2754514/1‑150 (MQ=255)
                                                                                cgcgGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGAttt                                                                       <  1:2448977/150‑1 (MQ=255)
                                                                                 gcggcgCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGAtttt                                                                      <  1:3492878/150‑1 (MQ=255)
                                                                                 gcggcgCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGAtttt                                                                      <  2:3189611/150‑1 (MQ=255)
                                                                                    gcgcTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCt                                                                   <  2:648737/150‑1 (MQ=255)
                                                                                    gcgcTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCt                                                                   <  1:1507762/150‑1 (MQ=255)
                                                                                      gcTGGAATGGCAGGCAGATTACGATATAAACCAGGCCGAAAGCTGGAGCCGTAGCATAGCAACGCTGGCGGAAGATGCGCTGGCGAAACTTCCCGGCATTAGATCAATCCGCTGACAGGATTCCACCCTGCTGGCCTTTGATTTTGCTgg                                                                 >  2:3482126/1‑150 (MQ=255)
                                                                                        tGGAATGGCTGGCAGATTACGATATCAACCGGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCg                                                               >  1:3368477/1‑150 (MQ=255)
                                                                                          gAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGtt                                                             <  2:2220399/150‑1 (MQ=255)
                                                                                                    cAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCg                                                   >  2:945184/1‑150 (MQ=255)
                                                                                                     aGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGa                                                  <  1:1298505/150‑1 (MQ=255)
                                                                                                      gATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGat                                                 <  2:3756213/150‑1 (MQ=255)
                                                                                                         tACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATg                                              <  1:4295819/150‑1 (MQ=255)
                                                                                                         tACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATg                                              <  1:1198834/150‑1 (MQ=255)
                                                                                                           cGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGt                                            <  2:2744746/150‑1 (MQ=255)
                                                                                                            gATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTg                                           <  1:4565826/150‑1 (MQ=255)
                                                                                                              tatCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGAc                                         >  1:915326/1‑150 (MQ=255)
                                                                                                               atCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACg                                        >  1:31174/1‑150 (MQ=255)
                                                                                                               atCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACg                                        >  2:530968/1‑150 (MQ=255)
                                                                                                                tCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACgc                                       >  2:1391709/1‑150 (MQ=255)
                                                                                                                 cAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACgct                                      >  1:2734544/1‑150 (MQ=255)
                                                                                                                 cAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACgct                                      >  1:3628234/1‑150 (MQ=255)
                                                                                                                       ggCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGc                                <  2:3824456/150‑1 (MQ=255)
                                                                                                                       ggCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGACTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGc                                >  2:686576/1‑150 (MQ=255)
                                                                                                                         ccGAAAGCTGGAGCCGTAGCTTAGCAACTCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCgg                              <  1:791139/150‑1 (MQ=255)
                                                                                                                            aaaGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGt                           <  2:198057/150‑1 (MQ=255)
                                                                                                                              aGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTAc                         <  1:2647691/150‑1 (MQ=255)
                                                                                                                                  ggAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTa                     >  1:397679/1‑150 (MQ=255)
                                                                                                                                       cGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGcc                >  1:1169830/1‑150 (MQ=255)
                                                                                                                                          aGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTg             >  2:3782345/1‑150 (MQ=255)
                                                                                                                                          aGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTg             >  2:1320801/1‑150 (MQ=255)
                                                                                                                                            gttAGCATCGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTGGTTCATTCCGCTGCAAGGATTCCAGCCTTCTGGCCTTTGATTTTGCTGGCGTTCTTCTTTGCGATTTGGTGACGCTGCTGGCGGAGTTGGGTTTTGCCCTGcg           <  2:2069456/149‑1 (MQ=255)
                                                                                                                                            cTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcg           <  1:3469593/150‑1 (MQ=255)
                                                                                                                                             ttAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcgg          >  2:889380/1‑150 (MQ=255)
                                                                                                                                             ttAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcgg          <  2:2679332/150‑1 (MQ=255)
                                                                                                                                              tAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcggg         >  1:2984704/1‑150 (MQ=255)
                                                                                                                                               aGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcgggc        >  2:148300/1‑150 (MQ=255)
                                                                                                                                                  aaCGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcgggccgg     >  1:1004063/1‑150 (MQ=255)
                                                                                                                                                    cGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcgggccgggc   <  2:3560183/150‑1 (MQ=255)
                                                                                                                                                    cGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGcgggccgggc   >  2:4532696/1‑150 (MQ=255)
                                                                                                                                                     gCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGCGGGCCGGgca  <  1:1728312/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTGACGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTGGAACGCCAAATGTCGCTGGTGTCATAGGATTAAGCGCGGCGCTGGAATGGCTGGCAGATTACGATATCAACCAGGCCGAAAGCTGGAGCCGTAGCTTAGCAACGCTGGCGGAAGATGCGCTGGCGAAACGTCCCGGCTTTCGTTCATTCCGCTGCCAGGATTCCAGCCTGCTGGCCTTTGATTTTGCTGGCGTTCATCATAGCGATATGGTGACGCTGCTGGCGGAGTACGGTATTGCCCTGCGGGCCGGGCA  >  CP009273/2937464‑2937762

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: