Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,707,179 A→G 5.3% F51F (TTT→TTC)  USA300HOU_2559 ← phytoene dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,707,1790AG 5.3% 233.2 / 2.4 77F51F (TTT→TTC) USA300HOU_2559phytoene dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (32/40);  new base G (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

TAGTCGTTCCATCTGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGCCACCACCAATAACGATG  >  CP000730/2707040‑2707317
                                                                                                                                           |                                                                                                                                          
tAGTCGTTCCATCTGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACATAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCACAAATATAAGGCATCGTTAAAATAGATGGACATAAATcaa                                                                                                                                           >  2:339701/1‑141 (MQ=255)
  gTCGTTCCATCTGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATaaaag                                                                                                                                         >  2:92997/1‑141 (MQ=255)
    cGTTCCATCTGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATcaaagcc                                                                                                                                       >  2:686637/1‑141 (MQ=255)
    cGTTCCATCTGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGAAATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCCTCGTTAAAATAGATGGACCTAAACcaaagcc                                                                                                                                       >  2:300772/1‑141 (MQ=255)
             tGGAAAAACGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATc                                                                                                                              <  1:742383/141‑1 (MQ=255)
             tGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTCAATCTAAGCCAAAGCCATc                                                                                                                              <  1:446485/141‑1 (MQ=255)
             tGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATc                                                                                                                              <  1:514837/141‑1 (MQ=255)
                    aaGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATACGGCCTCGTTAAAATAGATGGGCCTAAATCAAAGCCAAAGCCATCTGATTCa                                                                                                                       >  1:150074/1‑141 (MQ=255)
                      gCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATg                                                                                                                     >  2:25866/1‑141 (MQ=255)
                         cgcccATTGTGGGGGCATTCGTTTGATGGAAACGTATTCGGACTTTGGTTCCTGGCATTTTTCGATTATTTTTCAAAAAATATAAGGCATAGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGAGTCATGACg                                                                                                                  <  2:597122/138‑1 (MQ=255)
                          agCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCGGACATTTGATTCTTGCTATATTAGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGAGGGACCTAACGCAAAGCCAAAGCCATCTGAGTCAtgatga                                                                                                                 >  2:153444/2‑141 (MQ=255)
                           gCCATTGTTGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGAt                                                                                                                <  1:474069/141‑1 (MQ=255)
                             cATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTc                                                                                                              <  2:285385/141‑1 (MQ=255)
                             cATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTc                                                                                                              <  2:286332/141‑1 (MQ=255)
                              cttGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCa                                                                                                             <  2:602042/140‑1 (MQ=255)
                              actGATGCGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCa                                                                                                             <  1:110283/139‑1 (MQ=255)
                               ttGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCAc                                                                                                            <  2:154057/141‑1 (MQ=255)
                               ttGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTCGTTCTTGCTAAATTCGAATTATTTTTCAAAAATATAAGGCTTCGTTAAAATAGATGGACCTAAATCAAATCCATGGCTATCTTATTCATTACGATTCAc                                                                                                            <  2:223371/141‑1 (MQ=255)
                               ttGATGTGGCAATCGCTTGACTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTCAAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCAc                                                                                                            <  2:606355/141‑1 (MQ=255)
                                tGATGTGGCAATCGCTTGATTATAACGTAGTCTGACATTTGCTTCGTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACt                                                                                                           <  1:66467/141‑1 (MQ=255)
                                    gtgGCATTCGTTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGTCAAAGCCATCTGATTCATGACGATTCACTTTGc                                                                                                       <  2:50698/141‑1 (MQ=255)
                                    gtgGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGc                                                                                                       <  2:218253/141‑1 (MQ=255)
                                    gtgGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGc                                                                                                       >  2:192979/1‑141 (MQ=255)
                                    gtgGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGc                                                                                                       <  2:472695/141‑1 (MQ=255)
                                    gtgGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAACGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGc                                                                                                       <  2:691932/141‑1 (MQ=255)
                                     tgGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATACTCGAATAATTTTTCAAAAATATAAGGCCGCGGTAACTTAGATGGACCTAAATCAAAGCCAAAGCCATCGGATTCATTACGATTCACTTTGcc                                                                                                      >  2:271397/1‑141 (MQ=255)
                                       gcaatagttaaaTTGTCAAGGAGTCTGACATTTGCTGCTTGCTATATTCGAACATTTTTTCAAAAATATAAGACATCGGTATAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTAACTTTGcccc                                                                                                    >  2:385921/12‑139 (MQ=255)
                                        cAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGcctcc                                                                                                   <  2:185307/141‑1 (MQ=255)
                                          aTCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCta                                                                                                 >  2:242086/1‑141 (MQ=255)
                                           tCGCTTGTTTGTAACGTAGTATGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCGAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCtat                                                                                                >  1:507869/1‑141 (MQ=255)
                                           tCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCGAAGCCAAAGCCATCTGATTCATGACGATACACTTTGCCTCCtaa                                                                                                >  1:16785/1‑140 (MQ=255)
                                           tCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCtat                                                                                                >  2:335494/1‑141 (MQ=255)
                                           tCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCtat                                                                                                >  2:267811/1‑141 (MQ=255)
                                                tgcttGTCACGTATTTAGACACTTGTCTCCTGATATAATCGCAGACGTTTTCAAAAATATAAGGTATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGat                                                                                           <  1:167766/138‑1 (MQ=255)
                                                tGATTGTAACGTAGTCTGATATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGat                                                                                           <  2:346721/141‑1 (MQ=255)
                                                tGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATATTTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGat                                                                                           >  1:554797/1‑141 (MQ=255)
                                                tGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGTTCCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGat                                                                                           <  2:289203/141‑1 (MQ=255)
                                                tGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGat                                                                                           <  2:557000/141‑1 (MQ=255)
                                                tGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTCTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGat                                                                                           <  2:98951/141‑1 (MQ=255)
                                                   ttGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGattat                                                                                        <  1:16913/141‑1 (MQ=255)
                                                   ttGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGattat                                                                                        <  2:461005/141‑1 (MQ=255)
                                                   ttGTAACGGAGTCCGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGattat                                                                                        <  1:214020/141‑1 (MQ=255)
                                                         cGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGtt                                                                                  <  1:25866/141‑1 (MQ=255)
                                                           tAGTCTGACATTTGCTTCTTGCTATATTTGAATAATTTTTCAAAAATATAAGGCATCCTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCa                                                                                <  2:701245/141‑1 (MQ=255)
                                                           tAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCa                                                                                <  2:436285/141‑1 (MQ=255)
                                                           tAGTCAGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTTTATGATTCTTTTGtgct                                                                                >  2:460069/1‑138 (MQ=255)
                                                            aGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTAATAACGATTCACTTTGCCCCCGATATGATAATTTTGCTCat                                                                               >  2:213762/1‑141 (MQ=255)
                                                             gTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAACTCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCata                                                                              >  2:733876/1‑141 (MQ=255)
                                                               cTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCAtata                                                                            >  1:804934/1‑141 (MQ=255)
                                                                tGACATTTGCTTTTTTCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATaa                                                                           <  1:178566/141‑1 (MQ=255)
                                                                 gACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAAt                                                                          >  2:283181/1‑141 (MQ=255)
                                                                    aTTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATgag                                                                       <  1:91760/141‑1 (MQ=255)
                                                                    aTTTGCTTCTTGCTAGATTCGAATAATTTTGCAAAAAGATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATgag                                                                       <  1:395679/141‑1 (MQ=255)
                                                                       tGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGAcc                                                                    >  2:753842/1‑141 (MQ=255)
                                                                               gCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGATGAGcc                                                            >  1:619616/1‑141 (MQ=255)
                                                                                       cGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTtgtg                                                    >  1:460479/1‑141 (MQ=255)
                                                                                         aataatTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGcc                                                  >  2:160627/1‑141 (MQ=255)
                                                                                         aataatTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCAAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGcc                                                  >  2:432498/1‑141 (MQ=255)
                                                                                          ataatTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGcct                                                 >  2:70582/1‑140 (MQ=255)
                                                                                          ataatTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCa                                                 >  1:244788/1‑141 (MQ=255)
                                                                                          ataatTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCa                                                 >  1:320496/1‑141 (MQ=255)
                                                                                           taatTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCAt                                                >  1:133178/1‑141 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCTTTAAAATAGATTGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTTCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             <  2:221562/141‑1 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCGAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             <  2:16785/141‑1 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCGAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTATTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             <  2:507869/141‑1 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCAttgg                                             >  1:222728/1‑139 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             >  1:334913/1‑141 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             <  1:335494/141‑1 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             <  1:267811/141‑1 (MQ=255)
                                                                                              tttttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCg                                             <  2:626206/141‑1 (MQ=255)
                                                                                                tttCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGaa                                           >  2:696073/1‑141 (MQ=255)
                                                                                                    aaaaaGATACTGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGCTTCATTACCATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATtgc                                       <  1:385921/141‑1 (MQ=255)
                                                                                                       aaTATAATGCATCGTTAAAATAGATGTACCTACATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATtgctgc                                    <  2:173092/141‑1 (MQ=255)
                                                                                                              ggCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATc                             >  1:671065/1‑141 (MQ=255)
                                                                                                              ggCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATc                             >  2:265050/1‑141 (MQ=255)
                                                                                                                cATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATccc                           <  1:79824/141‑1 (MQ=255)
                                                                                                                            aGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATACGAATTGCTGCAGAAATCCCACCTAAGccacc               <  1:192979/141‑1 (MQ=255)
                                                                                                                                ggACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTCTATGATTATTTTGTTCATATAATGAGACCGAATAACCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGccaccacca           >  2:626867/1‑141 (MQ=255)
                                                                                                                                        aTCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGCCACCACCAATAACgat   >  1:261875/1‑141 (MQ=255)
                                                                                                                                        aTCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGCCACCACCAATAACgat   >  1:678072/1‑141 (MQ=255)
                                                                                                                                         tCAAAGCCAAAGCCATCTGATTTATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGCCACCACCAATAACgatg  <  2:452359/141‑1 (MQ=255)
                                                                                                                                         ccaaagccaaagcCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGCCACCACCAATAACgatg  <  2:602937/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TAGTCGTTCCATCTGGAAAAAAGCTACGCCATTGATGTGGCAATCGCTTGATTGTAACGTAGTCTGACATTTGCTTCTTGCTATATTCGAATAATTTTTCAAAAATATAAGGCATCGTTAAAATAGATGGACCTAAATCAAAGCCAAAGCCATCTGATTCATGACGATTCACTTTGCCTCCTATATGATTATTTTGTTCATATAATGAGACCGAATAGCCACTTTGTGCCATTCGAATTGCTGCAGAAATCCCACCTAAGCCACCACCAATAACGATG  >  CP000730/2707040‑2707317

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑