Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,882,392 C→A 15.2% W356L (TGG→TTG)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,882,3920CA15.2% 646.5 / 65.7 214W356L (TGG→TTG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base C (56/125);  new base A (14/19);  total (70/144)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

TCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTG  >  CP009273/2882245‑2882538
                                                                                                                                                   |                                                                                                                                                  
   ggTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGt                                                                                                                                                 <  1:559074/148‑1 (MQ=255)
    gTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGtt                                                                                                                                                <  2:2932996/148‑1 (MQ=255)
    gTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGtt                                                                                                                                                >  1:1228730/1‑148 (MQ=255)
     tACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTa                                                                                                                                               <  1:1661484/148‑1 (MQ=255)
      acacCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTat                                                                                                                                              <  2:3375618/148‑1 (MQ=255)
        acCGTCGTTTATCAGTTTCTTGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCATGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGTTGCGACCATTTTCGTTAAACAGCGTCAGGTGCCAGTTatca                                                                                                                                            <  1:3661612/148‑1 (MQ=255)
        acCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTatca                                                                                                                                            <  1:1233575/148‑1 (MQ=255)
           ggCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTAtcatca                                                                                                                                         <  1:1347138/148‑1 (MQ=255)
           ggCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTAtcatca                                                                                                                                         <  1:1734452/148‑1 (MQ=255)
            gCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCaa                                                                                                                                        <  1:1478599/148‑1 (MQ=255)
            gCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCaa                                                                                                                                        <  1:3377361/148‑1 (MQ=255)
             cGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAAt                                                                                                                                       <  2:140638/148‑1 (MQ=255)
              gatgatCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATg                                                                                                                                      <  1:83312/148‑1 (MQ=255)
                   tCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCtt                                                                                                                                 >  2:3797136/1‑148 (MQ=255)
                    cAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTa                                                                                                                                >  1:2712400/1‑148 (MQ=255)
                    cAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTa                                                                                                                                >  2:2059256/1‑148 (MQ=255)
                    cAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTa                                                                                                                                >  1:2846781/1‑148 (MQ=255)
                      gATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAAc                                                                                                                              >  2:3464623/1‑148 (MQ=255)
                      gATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAAc                                                                                                                              <  2:2958855/148‑1 (MQ=255)
                        ttCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTa                                                                                                                                               >  1:832624/1‑129 (MQ=255)
                        ttCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTa                                                                                                                                               <  2:832686/129‑1 (MQ=255)
                           tGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGc                                                                                                                         <  2:1186117/148‑1 (MQ=255)
                           tGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGc                                                                                                                         <  2:3755843/148‑1 (MQ=255)
                           tGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGc                                                                                                                         <  2:532284/148‑1 (MQ=255)
                               tCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAAt                                                                                                                     <  1:364699/148‑1 (MQ=255)
                               tCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAAt                                                                                                                     <  1:1835180/148‑1 (MQ=255)
                                   cGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGa                                                                                                                 <  2:3067710/148‑1 (MQ=255)
                                   cGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGa                                                                                                                 <  1:791963/148‑1 (MQ=255)
                                   cGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGa                                                                                                                 <  1:3803150/148‑1 (MQ=255)
                                   cGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGa                                                                                                                 <  2:1514833/148‑1 (MQ=255)
                                    gTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGAt                                                                                                                >  1:2519508/1‑148 (MQ=255)
                                    gTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGAt                                                                                                                <  1:710498/148‑1 (MQ=255)
                                        tGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGcc                                                                                                            >  1:2325445/1‑148 (MQ=255)
                                        tGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGcc                                                                                                            <  1:1920475/148‑1 (MQ=255)
                                        tGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGcc                                                                                                            >  2:757741/1‑148 (MQ=255)
                                         gCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCt                                                                                                           <  2:2701140/148‑1 (MQ=255)
                                             aaaTCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTc                                                                                                       >  2:3527656/1‑148 (MQ=255)
                                               aTCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCg                                                                                                     <  1:1102469/148‑1 (MQ=255)
                                                tCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCgg                                                                                                    <  1:2388911/148‑1 (MQ=255)
                                                tCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCgg                                                                                                    <  1:3486669/148‑1 (MQ=255)
                                                tCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCgg                                                                                                    <  1:3207188/148‑1 (MQ=255)
                                                tCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCgg                                                                                                    <  2:475476/148‑1 (MQ=255)
                                                      ttGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACt                                                                                              >  1:2200374/1‑148 (MQ=255)
                                                        gtgGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCa                                                                                            <  2:667977/148‑1 (MQ=255)
                                                         tgGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCat                                                                                           >  1:2257516/1‑148 (MQ=255)
                                                         tgGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCat                                                                                           >  1:328476/1‑148 (MQ=255)
                                                          ggATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                                                          <  2:820910/148‑1 (MQ=255)
                                                          ggATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                                                          <  1:3341871/148‑1 (MQ=255)
                                                           gagctacgcgATCTCCCGCAGGCCGCTAGACAGCGGACGCGCCCGAGAATCAAGGATGCGACCATGTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAAGGCCCGTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                                                         <  2:2615210/142‑1 (MQ=255)
                                                           gATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                                                         <  2:2615166/148‑1 (MQ=255)
                                                           gATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                                                         >  1:2783758/1‑148 (MQ=255)
                                                           gATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                                                         >  1:1437520/1‑148 (MQ=255)
                                                             tCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATgg                                                                                       <  1:56312/148‑1 (MQ=255)
                                                              cTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGa                                                                                      >  2:2286098/1‑148 (MQ=255)
                                                                tCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACg                                                                                    >  1:217115/1‑148 (MQ=255)
                                                                tCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACg                                                                                    <  1:2920563/148‑1 (MQ=255)
                                                                tCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACg                                                                                    <  1:1011282/148‑1 (MQ=255)
                                                                 cgcgATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACgg                                                                                   <  2:2783758/148‑1 (MQ=255)
                                                                 cgcgATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACgg                                                                                   <  1:25480/148‑1 (MQ=255)
                                                                  gcgATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGa                                                                                  <  2:1651668/148‑1 (MQ=255)
                                                                  gcgATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGa                                                                                  <  1:3688296/148‑1 (MQ=255)
                                                                   cgATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGAt                                                                                 <  2:1516984/148‑1 (MQ=255)
                                                                     aTCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCg                                                                               <  1:2217793/148‑1 (MQ=255)
                                                                      tctcCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCgg                                                                              <  2:3643003/148‑1 (MQ=255)
                                                                       ctcCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGt                                                                             <  1:760682/148‑1 (MQ=255)
                                                                        tcCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGtt                                                                            <  1:1125111/148‑1 (MQ=255)
                                                                                gCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTg                                                                    <  1:3565779/148‑1 (MQ=255)
                                                                                gCCGGTTTTCAGAGGACGCGCCGTATAAACAAGGATTCGACAATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACACATATGGACGGATCGGTTCAAATTTg                                                                    >  2:840818/1‑148 (MQ=255)
                                                                                 ccgtatttAAGCTGACGCGCAGGATAATCAAGGATGCGACAAATTTCGATAAACAACGTCAGGTGCCAGTTATAATCAATCCCCTAAACCCAGCCAATACGATCGCAACGTACGGTGAACTCATATGGACGGATCGGTTCAAATTAGa                                                                   >  2:2238425/6‑148 (MQ=255)
                                                                                 ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGAGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGa                                                                   <  2:2164771/148‑1 (MQ=255)
                                                                                    tttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATcc                                                                <  1:1397005/147‑1 (MQ=255)
                                                                                     ttttCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATccc                                                               <  2:2499880/148‑1 (MQ=255)
                                                                                            cGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGa                                                        <  2:2257472/148‑1 (MQ=255)
                                                                                              gACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCACTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACg                                                      <  2:3117921/148‑1 (MQ=255)
                                                                                                cgcgCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCt                                                    <  2:2458512/148‑1 (MQ=255)
                                                                                                      ggATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTAATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACtt                                              <  2:2431532/148‑1 (MQ=255)
                                                                                                       gATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTc                                             >  1:1694941/1‑148 (MQ=255)
                                                                                                        aTAATCAAGGATGCGACCATTTCCGATAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTcc                                            >  2:1533999/1‑148 (MQ=255)
                                                                                                             ctaGGATGCGACCATTTTCGAGAAACAGCGTCATGTGCCAGTTTTCATCAATGCCCTTCACCCAGCATATCCGATCTCCTAGTCCGGTGAACTCATATGGACGGATCGGTTCATATTTTAGCCCCGCGAGACGATCCACTTCCGCttt                                       <  2:688819/146‑1 (MQ=255)
                                                                                                             cAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCttt                                       >  1:3224536/1‑148 (MQ=255)
                                                                                                             cAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCttt                                       <  2:3486565/148‑1 (MQ=255)
                                                                                                                ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaaa                                    <  1:862627/148‑1 (MQ=255)
                                                                                                                ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaaa                                    <  1:606793/148‑1 (MQ=255)
                                                                                                                   tGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGt                                 <  2:1234419/148‑1 (MQ=255)
                                                                                                                   tGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGt                                 <  2:3599675/148‑1 (MQ=255)
                                                                                                                       aCCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCa                             >  2:158303/1‑148 (MQ=255)
                                                                                                                            tttCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAcccc                        >  1:2164897/1‑148 (MQ=255)
                                                                                                                             ttCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAccccc                       <  2:3444494/148‑1 (MQ=255)
                                                                                                                                 aTAAACAGCGTCAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcgc                   >  1:495443/1‑148 (MQ=255)
                                                                                                                                    aaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCt                <  2:3431584/148‑1 (MQ=255)
                                                                                                                                     aCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTc               <  2:2765576/148‑1 (MQ=255)
                                                                                                                                         cGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGc           >  2:1518033/1‑148 (MQ=255)
                                                                                                                                         cGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGc           >  2:1974260/1‑148 (MQ=255)
                                                                                                                                          gTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCg          <  1:264437/148‑1 (MQ=255)
                                                                                                                                            cAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtg        >  2:125944/1‑148 (MQ=255)
                                                                                                                                            cAGGTGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtg        >  1:2907355/1‑148 (MQ=255)
                                                                                                                                             aGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtgt       <  1:3630579/148‑1 (MQ=255)
                                                                                                                                                tGCAAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAtt    <  2:2523243/148‑1 (MQ=255)
                                                                                                                                                 gCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAttt   <  1:2836640/148‑1 (MQ=255)
                                                                                                                                                 gCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAttt   <  1:3744156/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTG  >  CP009273/2882245‑2882538

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: