Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 230,579 C→T 100% E231K (GAA→AAA)  gloB ← hydroxyacylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273230,5790CT93.9% 1024.8 / 41.5 294E231K (GAA→AAA) gloBhydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (11/7);  new base T (134/142);  total (145/149)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 299 total aligned reads displayed.

TTTGCCTTCATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATGGACAAATCG  >  CP009273/230432‑230727
                                                                                                                                                   |                                                                                                                                                    
tttGCCTTCATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTttctt                                                                                                                                                    <  2:2797873/150‑1 (MQ=255)
  tGCCTTCATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTca                                                                                                                                                  >  1:4189897/1‑150 (MQ=255)
    ccTTCATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATAAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCTGGTTGTTGCAATAATGTTTTTTcatt                                                                                                                                                >  2:2730042/1‑150 (MQ=255)
     cTTCATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcatta                                                                                                                                               <  1:1731970/150‑1 (MQ=255)
        cATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaatt                                                                                                                                            <  1:267336/150‑1 (MQ=255)
        cATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaatt                                                                                                                                            >  2:4790668/1‑150 (MQ=255)
             gtgtgtCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaattacatt                                                                                                                                       >  1:668037/1‑150 (MQ=255)
              tgtgtCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGATTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaattacatta                                                                                                                                      >  2:4890437/1‑150 (MQ=255)
               gtgtCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaattacattaa                                                                                                                                     >  1:4111599/1‑150 (MQ=255)
                tgtCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaattacattaat                                                                                                                                    >  2:1273346/1‑150 (MQ=255)
                  tCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTcattaattacattaatta                                                                                                                                  >  2:4967203/1‑150 (MQ=255)
                   cAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACAttaattaa                                                                                                                                 <  1:2856127/150‑1 (MQ=255)
                   cAATAGTAGCTTAAAAGACGACCGAACATAACGGCGACCGGAGCCGAAGACCAGAAAGTTTTAGCAGAACCTATCTTTCTTTGAACTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACAttaattaa                                                                                                                                 <  2:3678769/150‑1 (MQ=255)
                     aTAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAAt                                                                                                                               >  1:3764214/1‑150 (MQ=255)
                         ttGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAat                                                                                                                           >  2:4899643/1‑150 (MQ=255)
                         ttGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAat                                                                                                                           >  2:3085616/1‑150 (MQ=255)
                          tGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAata                                                                                                                          >  2:2086751/1‑150 (MQ=255)
                           gCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAatat                                                                                                                         <  2:430903/150‑1 (MQ=255)
                            cTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATc                                                                                                                        <  1:4922005/150‑1 (MQ=255)
                             ttAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCt                                                                                                                       >  2:841878/1‑150 (MQ=255)
                             ttAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCt                                                                                                                       >  2:3419260/1‑150 (MQ=255)
                              tAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCtt                                                                                                                      <  2:4033219/150‑1 (MQ=255)
                              tAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCtt                                                                                                                      >  1:3523493/1‑150 (MQ=255)
                                aaaGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTcc                                                                                                                    <  2:2488883/150‑1 (MQ=255)
                                aaaGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTcc                                                                                                                    <  2:4916347/150‑1 (MQ=255)
                                aaaGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTcc                                                                                                                    <  2:1773766/150‑1 (MQ=255)
                                 aaGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCg                                                                                                                   <  1:345947/150‑1 (MQ=255)
                                    acgacCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTACCCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTc                                                                                                                <  2:2868770/150‑1 (MQ=255)
                                      gacCGATCATAACGGCGAACGGAGCCTATGACAAGAAAGTTTTATCAGAACATATCTTTCTTTGACATTAACCATGAAAAACGCTCTTCAGGTTTTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTACGTTAtt                                                                                                              >  1:2186256/1‑150 (MQ=255)
                                      gacCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCtt                                                                                                              >  2:1346307/1‑150 (MQ=255)
                                       acCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTa                                                                                                             >  1:3255023/1‑150 (MQ=255)
                                          gATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTaaaa                                                                                                          <  1:4617300/150‑1 (MQ=255)
                                           aTCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTaaaaa                                                                                                         <  1:1633732/150‑1 (MQ=255)
                                                 aCGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACAtt                                                                                                   >  2:1034955/1‑150 (MQ=255)
                                                   ggCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTaa                                                                                                 <  2:4738097/150‑1 (MQ=255)
                                                   ggCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTaa                                                                                                 >  1:1303688/1‑150 (MQ=255)
                                                     cGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAAtt                                                                                               <  1:1142153/150‑1 (MQ=255)
                                                      gAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAAttt                                                                                              <  2:4248418/150‑1 (MQ=255)
                                                       aaCGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTg                                                                                             >  1:4253624/1‑150 (MQ=255)
                                                       aaCGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTg                                                                                             >  2:1881360/1‑150 (MQ=255)
                                                       aaCGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTg                                                                                             <  2:744274/150‑1 (MQ=255)
                                                        aCGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGc                                                                                            <  2:1333618/150‑1 (MQ=255)
                                                            aGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCt                                                                                        >  1:3343554/1‑150 (MQ=255)
                                                              ccGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCa                                                                                      <  1:1186641/150‑1 (MQ=255)
                                                              ccGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCa                                                                                      <  1:4419061/150‑1 (MQ=255)
                                                               cGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAt                                                                                     <  1:2319104/150‑1 (MQ=255)
                                                                gATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtt                                                                                    <  2:1551222/150‑1 (MQ=255)
                                                                gATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtt                                                                                    >  2:195558/1‑150 (MQ=255)
                                                                gATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtt                                                                                    <  1:1209339/150‑1 (MQ=255)
                                                                gATGAAAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACTCTATTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtt                                                                                    <  1:264258/150‑1 (MQ=255)
                                                                 aTGACAATAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAttt                                                                                   <  2:3322789/150‑1 (MQ=255)
                                                                  tGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCAtttt                                                                                  <  2:3652099/150‑1 (MQ=255)
                                                                    aCAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTc                                                                                >  2:1048174/1‑150 (MQ=255)
                                                                     cAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCa                                                                               >  2:621400/1‑150 (MQ=255)
                                                                         aaaGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAAt                                                                           >  1:2451577/1‑150 (MQ=255)
                                                                          aaGTTTTATCATAACATATCTTTATTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCTTTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAAtt                                                                          <  1:765750/150‑1 (MQ=255)
                                                                                tATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGt                                                                    >  1:3873872/1‑150 (MQ=255)
                                                                                 aTCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTa                                                                   >  1:729725/1‑150 (MQ=255)
                                                                                    aGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAgtg                                                                <  2:1439554/150‑1 (MQ=255)
                                                                                     gAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAgtgt                                                               >  1:704334/1‑150 (MQ=255)
                                                                                       aCCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTa                                                             <  2:950241/150‑1 (MQ=255)
                                                                                       aCCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTa                                                             >  1:5105541/1‑150 (MQ=255)
                                                                                        ccTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAt                                                            >  2:362953/1‑150 (MQ=255)
                                                                                        ccTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAt                                                            >  1:4772911/1‑150 (MQ=255)
                                                                                         cTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAtt                                                           <  2:4642561/150‑1 (MQ=255)
                                                                                          tATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTTATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTAttt                                                          <  1:5143838/150‑1 (MQ=255)
                                                                                           aTCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTg                                                         >  2:3731557/1‑150 (MQ=255)
                                                                                              tttctttGACCTTAACCATGCAAAAGGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGAtt                                                      <  1:3161610/150‑1 (MQ=255)
                                                                                                tctttGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGAtttt                                                    <  1:3600471/150‑1 (MQ=255)
                                                                                                   tttaCCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTg                                                 >  2:4410157/4‑150 (MQ=255)
                                                                                                      aCCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGccc                                              <  2:4548695/150‑1 (MQ=255)
                                                                                                        cTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGt                                            >  1:300854/1‑150 (MQ=255)
                                                                                                          tAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTaa                                          >  1:4074876/1‑150 (MQ=255)
                                                                                                            aCCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACt                                        <  2:4849011/150‑1 (MQ=255)
                                                                                                              cATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTcc                                      >  2:2534636/1‑150 (MQ=255)
                                                                                                               aTGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCt                                     <  2:2229188/150‑1 (MQ=255)
                                                                                                                tGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCtt                                    >  2:1904151/1‑150 (MQ=255)
                                                                                                                     aaCGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACtt                               >  2:2264698/1‑150 (MQ=255)
                                                                                                                     aaCGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACtt                               <  2:4394749/150‑1 (MQ=255)
                                                                                                                       cGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTa                             <  2:1024066/150‑1 (MQ=255)
                                                                                                                       cGCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTa                             >  2:3134043/1‑150 (MQ=255)
                                                                                                                        gCTCTTCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTAc                            <  2:16785/150‑1 (MQ=255)
                                                                                                                          tctTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGa                          <  2:3591954/150‑1 (MQ=255)
                                                                                                                            ttCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGata                        >  1:2102472/1‑150 (MQ=255)
                                                                                                                            ttCAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGata                        >  1:377814/1‑150 (MQ=255)
                                                                                                                              cAGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGataat                      >  2:1855642/1‑150 (MQ=255)
                                                                                                                               aGGTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGataata                     >  2:2539441/1‑150 (MQ=255)
                                                                                                                                ggTTGTTGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGAtaataa                    <  2:260041/150‑1 (MQ=255)
                                                                                                                                ggTTGTTGCAATAATGTTTTTTCATTAATTAAATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGAtaataa                    <  1:1939002/150‑1 (MQ=255)
                                                                                                                                 gttgttGCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATaataat                   >  1:2565422/1‑150 (MQ=255)
                                                                                                                                   tgttgCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCa                 >  1:1852503/1‑150 (MQ=255)
                                                                                                                                    gttgCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGGGGTCTTTGATTTTTTGCCCGTAACTCCTGAACTTTACGATAATAATCAt                >  2:4379715/1‑150 (MQ=255)
                                                                                                                                       gCAATAATGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTa             >  1:845997/1‑150 (MQ=255)
                                                                                                                                        cAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTAt            >  2:3854226/1‑150 (MQ=255)
                                                                                                                                         aataatGTTTTTTCATTAAGTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATg           <  2:1603547/150‑1 (MQ=255)
                                                                                                                                          ataatGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATgg          <  2:5082349/150‑1 (MQ=255)
                                                                                                                                            aatGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATGGAc        <  1:450059/150‑1 (MQ=255)
                                                                                                                                             atGTTTTTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATGGACa       >  1:915663/1‑150 (MQ=255)
                                                                                                                                                ttttttCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATGGACAAAt    <  1:1346307/150‑1 (MQ=255)
                                                                                                                                                  tcttcATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATGGACAAATCg  <  1:2622113/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTGCCTTCATGTGTGTGTCAATAGTTGCTTAAAAGACGACCGATCATAACGGCGAACGGAGCCGATGACAAGAAAGTTTTATCAGAACCTATCTTTCTTTGACCTTAACCATGCAAAACGCTCTTCAGGTTGTTGCAATAATGTTTCTTCATTAATTACATTAATTAAATCAATATCTTCCGTTCTTAAAAAAACATTAATTTGCCGCTCATTTTTCAGAATTACGGGTAGTGTTATTTGATTTTTTGCCCGTAACTCCTTAACTTTACGATAATAATCATTTATGGACAAATCG  >  CP009273/230432‑230727

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: