Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,464,329 A→T 5.9% Y127* (TAT→TAA treC ← trehalose‑6‑phosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,464,3290AT 5.9% 365.1 / 9.7 118Y127* (TAT→TAAtreCtrehalose‑6‑phosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (62/49);  new base T (4/3);  total (66/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 121 total aligned reads displayed.

GCTGGATTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCGTCCAGCGTG  >  NC_000913/4464180‑4464476
                                                                                                                                                     |                                                                                                                                                   
gCTGGATTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGata                                                                                                                                                     <  2:1057520/150‑1 (MQ=255)
gCTGGATTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGata                                                                                                                                                     >  2:1342841/1‑150 (MQ=255)
    gATTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTCGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATaaa                                                                                                                                                 <  2:78075/150‑1 (MQ=255)
     aTTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATaaaa                                                                                                                                                <  1:1127629/150‑1 (MQ=255)
     aTTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATaaaa                                                                                                                                                <  1:1179463/150‑1 (MQ=255)
         tcCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTgg                                                                                                                                            >  2:510169/1‑150 (MQ=255)
         tcCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTgg                                                                                                                                            <  2:1121166/150‑1 (MQ=255)
         tcCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTgg                                                                                                                                            >  2:1190139/1‑150 (MQ=255)
          cccAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGc                                                                                                                                           <  1:433230/150‑1 (MQ=255)
           ccAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCg                                                                                                                                          <  2:59500/150‑1 (MQ=255)
              gTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTa                                                                                                                                       >  1:574233/1‑150 (MQ=255)
                tGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAg                                                                                                                                     >  1:442043/1‑150 (MQ=255)
                 gagaTCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGTTTTATATCTGGCgcgatgg                                                                                                                                    <  2:880491/150‑7 (MQ=255)
                  agaTCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGTTTTATATCTGGCgcgatgga                                                                                                                                   <  1:969740/150‑8 (MQ=255)
                   gaTCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGc                                                                                                                                  >  2:649104/1‑150 (MQ=255)
                     tCCGCCTGTTCTGGTGCATAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGctt                                                                                                                                >  1:227802/1‑150 (MQ=255)
                      ccGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGcttt                                                                                                                               >  2:1254372/1‑150 (MQ=255)
                           tGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTg                                                                                                                          <  1:772923/150‑1 (MQ=255)
                           tGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTg                                                                                                                          <  1:845821/150‑1 (MQ=255)
                           tGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTg                                                                                                                          <  1:387035/150‑1 (MQ=255)
                           tGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTg                                                                                                                          <  1:241977/150‑1 (MQ=255)
                           tGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTg                                                                                                                          >  2:925075/1‑150 (MQ=255)
                              tCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTc                                                                                                                       <  2:410666/150‑1 (MQ=255)
                               cTGGTGCAAAGAGATGCAAATAGTACTGATCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAgg                                                                                                                                    >  2:492795/1‑136 (MQ=255)
                                tGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAg                                                                                                                     <  1:670649/150‑1 (MQ=255)
                                   tGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCTGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAgcgc                                                                                                                  >  1:1425740/1‑150 (MQ=255)
                                    gCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGcc                                                                                                                 >  2:74426/1‑150 (MQ=255)
                                     cAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCt                                                                                                                >  1:884356/1‑150 (MQ=255)
                                         gagaTGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTcgcg                                                                                                            <  2:764623/150‑1 (MQ=255)
                                               cAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCa                                                                                                      <  1:452433/150‑1 (MQ=255)
                                                aaaTAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAg                                                                                                     <  2:331549/150‑1 (MQ=255)
                                                aaaTAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAg                                                                                                     >  1:1183288/1‑150 (MQ=255)
                                                      tACTGTTGGCTTTCCGCATGCCAGCGCCACGCACTACTGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGt                                                                                               <  1:605292/150‑1 (MQ=255)
                                                       aCTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGtt                                                                                              >  1:224948/1‑150 (MQ=255)
                                                        cTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTg                                                                                             >  2:1029085/1‑150 (MQ=255)
                                                        cTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGGGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTg                                                                                             <  2:18529/150‑1 (MQ=255)
                                                          gTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTggg                                                                                           >  2:1192542/1‑150 (MQ=255)
                                                          gTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTggg                                                                                           >  1:1315216/1‑150 (MQ=255)
                                                           ttCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGt                                                                                          <  1:1190139/150‑1 (MQ=255)
                                                            tCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTa                                                                                         <  2:331490/150‑1 (MQ=255)
                                                            tCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTa                                                                                         >  1:1107840/1‑150 (MQ=255)
                                                             cGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTag                                                                                        >  2:758033/1‑150 (MQ=255)
                                                               cTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTagag                                                                                      >  2:1149042/1‑150 (MQ=255)
                                                               cTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTagag                                                                                      >  2:196436/1‑150 (MQ=255)
                                                                 ttCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAATCTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGt                                                                                    <  2:580194/150‑1 (MQ=255)
                                                                  tCCGCATGCCAGCGGCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTa                                                                                   >  1:391365/1‑150 (MQ=255)
                                                                  tCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGc                                                                                                                                           >  1:1041030/1‑94 (MQ=255)
                                                                     gCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATgg                                                                                <  2:574233/150‑1 (MQ=255)
                                                                       aTGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGtt                                                                              <  1:1378776/150‑1 (MQ=255)
                                                                        tGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTa                                                                             >  1:1218233/1‑150 (MQ=255)
                                                                          ccAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTaaa                                                                           >  1:1376474/1‑150 (MQ=255)
                                                                            aGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGTGGTATGGTTAAaca                                                                         <  1:602628/150‑1 (MQ=255)
                                                                                ccACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCat                                                                     <  2:1218233/150‑1 (MQ=255)
                                                                                 cACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCata                                                                    >  2:538314/1‑150 (MQ=255)
                                                                                    gCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCg                                                                 <  2:634123/150‑1 (MQ=255)
                                                                                      aCTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCgag                                                               <  2:442043/150‑1 (MQ=255)
                                                                                       cTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCgaga                                                              >  1:964231/1‑150 (MQ=255)
                                                                                           cGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGctttc                                                                                                                              >  1:1242272/1‑82 (MQ=255)
                                                                                           cGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGctttc                                                                                                                              <  2:1242272/82‑1 (MQ=255)
                                                                                            gCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGAt                                                         >  1:375687/1‑150 (MQ=255)
                                                                                            gCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGAt                                                         >  1:445044/1‑150 (MQ=255)
                                                                                              cAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATAc                                                       <  1:196436/150‑1 (MQ=255)
                                                                                                 aTTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGaa                                                    >  2:625671/1‑150 (MQ=255)
                                                                                                 aTTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGaa                                                    >  1:403811/1‑150 (MQ=255)
                                                                                                       aaCGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCcgc                                              >  1:1369284/1‑150 (MQ=255)
                                                                                                         cGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCcgcgc                                            <  2:803456/150‑1 (MQ=255)
                                                                                                         cGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCcgcgc                                            >  1:472173/1‑150 (MQ=255)
                                                                                                         cGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCcgcgc                                            >  2:639304/1‑150 (MQ=255)
                                                                                                                tgttCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTg                                     <  1:783031/150‑1 (MQ=255)
                                                                                                                tgttCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTg                                     >  1:865566/1‑150 (MQ=255)
                                                                                                                tgttCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTg                                     <  1:925075/150‑1 (MQ=255)
                                                                                                                  ttCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGcc                                   <  2:738412/150‑1 (MQ=255)
                                                                                                                     ggTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGAttt                                       >  1:1029399/1‑143 (MQ=255)
                                                                                                                     ggTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGAttt                                       <  2:1029399/143‑1 (MQ=255)
                                                                                                                     ggTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGc                                >  1:290318/1‑150 (MQ=255)
                                                                                                                     ggTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGc                                >  2:464107/1‑150 (MQ=255)
                                                                                                                         gCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCa                            >  1:509712/1‑150 (MQ=255)
                                                                                                                          cGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCAc                           >  1:1151961/1‑150 (MQ=255)
                                                                                                                           gTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCAcc                          <  1:583446/150‑1 (MQ=255)
                                                                                                                            tttCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCa                         <  1:17303/150‑1 (MQ=255)
                                                                                                                              tCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGt                       >  1:216352/1‑150 (MQ=255)
                                                                                                                               cTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGAttt                                       >  2:546561/1‑133 (MQ=255)
                                                                                                                               cTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGAttt                                       <  1:546561/133‑1 (MQ=255)
                                                                                                                                 ggTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACGATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCg                    >  1:615332/1‑150 (MQ=255)
                                                                                                                                    tctcCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCa                 <  2:280212/150‑1 (MQ=255)
                                                                                                                                     ctcCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCaa                >  2:883399/1‑150 (MQ=255)
                                                                                                                                     ctcCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCaa                <  1:74426/150‑1 (MQ=255)
                                                                                                                                      tcCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCaaa               >  1:851139/1‑150 (MQ=255)
                                                                                                                                      cccttaccgcCAGTTTTATATCTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCaaa               >  2:614502/8‑150 (MQ=255)
                                                                                                                                       ccATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCaaaa              >  1:795375/1‑150 (MQ=255)
                                                                                                                                        cttaccgcCAGTTTTATATCTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAAt             >  2:1009200/6‑150 (MQ=255)
                                                                                                                                        cATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAAt             <  2:509712/150‑1 (MQ=255)
                                                                                                                                        cATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAAt             >  2:1426911/1‑150 (MQ=255)
                                                                                                                                         aTCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATc            >  1:248371/1‑150 (MQ=255)
                                                                                                                                          taccgcCAGTTTTATATCTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCg           >  1:175478/4‑150 (MQ=255)
                                                                                                                                          tCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCg           <  2:1376474/150‑1 (MQ=255)
                                                                                                                                            ccgcCAGTTTTATATCTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCGTc         >  1:568218/2‑150 (MQ=255)
                                                                                                                                             ctccAGATATAAAACTGGCGGTAATGGCTTTCTTTTTTCAGCTCCTCGCGAAACCATGCATGTTGGTTAGAGTTATTGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCTTcc        <  1:1309355/148‑1 (MQ=255)
                                                                                                                                                cAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCGTCCAGc     <  1:639304/150‑1 (MQ=255)
                                                                                                                                                 aGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCGTCCAGCg    <  2:1107840/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTGGATTCTCCCAGTTGAGATCCGCCTGTTCTGGTGCAAAGAGATGCAAATAGTACTGTTCGCTTTCCGCATGCCAGCGCCACGCACTACCGCCAAATTTTGAACGCCAGTTGTTCGGTGGCGTTTCTGGTTCTCCATCGCGCCAGATATAAAACTGGCGGTAAGGGCTTTCTTTGTTCAGCGCCTCGCGAAACCAGGCATGTTGGGTAGAGGTATGGTTAAACACCATATCGAGAATGATACGAATCCCGCGCGATTTTGCCTGCGTCACCAGTTCGTCAAAATCGTCCAGCGTG  >  NC_000913/4464180‑4464476

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