Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,100 C→A 12.9% F89L (TTC→TTA crr → glucose‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,536,1000CA12.9% 395.4 / 35.7 123F89L (TTC→TTAcrrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (49/58);  new base A (8/8);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 127 total aligned reads displayed.

TCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGG  >  NC_000913/2535952‑2536248
                                                                                                                                                    |                                                                                                                                                    
tCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc                                                                                                                                                      <  4:709678/149‑1 (MQ=255)
          ggAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCg                                                                                                                                            >  3:334133/1‑149 (MQ=255)
            aaaaaaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGt                                                                                                                                          >  4:54678/1‑149 (MQ=255)
             aaaaaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTa                                                                                                                                         >  4:901571/1‑149 (MQ=255)
             aaaaaTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCGTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGCTa                                                                                                                                         >  2:7936/1‑149 (MQ=255)
                 aTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGa                                                                                                                                     <  2:507762/149‑1 (MQ=255)
                 aTCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGa                                                                                                                                     <  2:5889/149‑1 (MQ=255)
                    gTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGacac                                                                                                                                  <  3:479399/149‑1 (MQ=255)
                      tGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCg                                                                                                                                <  1:479030/149‑1 (MQ=255)
                        gTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGtt                                                                                                                              <  1:892349/149‑1 (MQ=255)
                         tGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTg                                                                                                                             <  2:397533/149‑1 (MQ=255)
                         tGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTg                                                                                                                             >  2:589698/1‑149 (MQ=255)
                           aTGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGaa                                                                                                                           <  4:700018/149‑1 (MQ=255)
                                aTTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaa                                                                                                                      <  3:54678/149‑1 (MQ=255)
                                 ttGCTATCAAACCAATGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGaaa                                                                                                                     >  4:905026/1‑149 (MQ=255)
                                   gCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAgg                                                                                                                   <  1:867182/149‑1 (MQ=255)
                                     tATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCg                                                                                                                 >  3:689954/1‑149 (MQ=255)
                                     tATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCg                                                                                                                 >  4:326322/1‑149 (MQ=255)
                                      aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGa                                                                                                                >  1:638804/1‑149 (MQ=255)
                                      aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGa                                                                                                                >  3:109747/1‑149 (MQ=255)
                                      aTCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGa                                                                                                                >  4:261065/1‑149 (MQ=255)
                                       tCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGaa                                                                                                               >  4:547241/1‑149 (MQ=255)
                                        cAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAg                                                                                                              <  4:109747/149‑1 (MQ=255)
                                        cAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAg                                                                                                              <  4:391712/149‑1 (MQ=255)
                                         aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGacac                                                                                                                                  <  3:421986/128‑1 (MQ=255)
                                         aaaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGacac                                                                                                                                  >  4:421986/1‑128 (MQ=255)
                                          aaCCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGc                                                                                                            >  4:877225/1‑149 (MQ=255)
                                             cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                                                                                         >  1:161337/1‑149 (MQ=255)
                                             cAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                                                                                         >  4:675051/1‑149 (MQ=255)
                                                  ggTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCg                                                                                                    >  4:2614/1‑149 (MQ=255)
                                                      aCAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTAtt                                                                                                <  4:870134/149‑1 (MQ=255)
                                                        aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCaa                                                                                                       >  3:721210/1‑140 (MQ=255)
                                                        aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCaa                                                                                                       <  4:721210/140‑1 (MQ=255)
                                                        aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGc                                                                                              >  2:691208/1‑149 (MQ=255)
                                                        aaaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGc                                                                                              >  1:486826/1‑149 (MQ=255)
                                                         aaaTGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCt                                                                                             >  2:240231/1‑149 (MQ=255)
                                                             ggTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTgaag                                                                                         <  2:161337/149‑1 (MQ=255)
                                                              gTCGCGCCAGTAGGCGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTgaaga                                                                                        >  3:857346/1‑149 (MQ=255)
                                                                cgcgCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                                                                                      <  2:232726/149‑1 (MQ=255)
                                                                cgcgCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                                                                                      <  3:2614/149‑1 (MQ=255)
                                                                cgcgCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCAAGCACTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                                                                                      <  1:862285/149‑1 (MQ=255)
                                                                              gCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGaaa                                                                        <  2:638804/149‑1 (MQ=255)
                                                                                 ccATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGtt                                                                     <  4:595304/149‑1 (MQ=255)
                                                                                  cATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTCCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTg                                                                    <  4:524490/149‑1 (MQ=255)
                                                                                   aTTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTgg                                                                   <  2:410304/149‑1 (MQ=255)
                                                                                    ttGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGc                                                                  >  2:563556/1‑149 (MQ=255)
                                                                                      ggTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGa                                                                <  3:473413/149‑1 (MQ=255)
                                                                                       gTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGAt                                                               >  1:500485/1‑149 (MQ=255)
                                                                                       gTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGAt                                                               >  1:890084/1‑149 (MQ=255)
                                                                                         aaaaTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATAc                                                             <  4:334133/149‑1 (MQ=255)
                                                                                         aaaaTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATAc                                                             <  3:696917/149‑1 (MQ=255)
                                                                                          aaaTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACt                                                            >  4:210126/1‑149 (MQ=255)
                                                                                           aaTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTg                                                           >  2:139655/1‑149 (MQ=255)
                                                                                            aTCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGt                                                          >  3:152953/1‑149 (MQ=255)
                                                                                             tCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTc                                                         >  2:348862/1‑149 (MQ=255)
                                                                                               tttGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCAt                                                       >  1:635201/1‑149 (MQ=255)
                                                                                               tttGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCAt                                                       <  2:711912/149‑1 (MQ=255)
                                                                                               tttGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCAt                                                       >  1:635862/1‑149 (MQ=255)
                                                                                               tttGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCAt                                                       >  4:867944/1‑149 (MQ=255)
                                                                                                    aaccaaccACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAt                                                  <  1:589698/149‑1 (MQ=255)
                                                                                                    aaccaaccACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCaa                                                                                                       <  4:432065/96‑1 (MQ=255)
                                                                                                     accaaccaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGaagaag                                                                                      <  3:182854/112‑1 (MQ=255)
                                                                                                     accaaccaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAtt                                                 <  3:680505/149‑1 (MQ=255)
                                                                                                        aaccaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGa                                              <  1:691208/149‑1 (MQ=255)
                                                                                                        aaccaCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGa                                              >  1:96827/1‑149 (MQ=255)
                                                                                                           caCGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCt                                           <  3:874796/149‑1 (MQ=255)
                                                                                                             cGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGc                                         >  4:477930/1‑149 (MQ=255)
                                                                                                             cGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGc                                         <  4:677949/149‑1 (MQ=255)
                                                                                                              gCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGcc                                        >  4:410536/1‑149 (MQ=255)
                                                                                                              gCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGcc                                        <  4:448826/149‑1 (MQ=255)
                                                                                                                aTTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                                                                                         >  2:419075/1‑82 (MQ=255)
                                                                                                                aTTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTc                                                                                                         <  1:419075/82‑1 (MQ=255)
                                                                                                                   ctctATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCgctgc                                   <  4:689954/149‑1 (MQ=255)
                                                                                                                    tctATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGc                                                                  <  3:172497/117‑1 (MQ=255)
                                                                                                                    tctATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGc                                                                  >  4:172497/1‑117 (MQ=255)
                                                                                                                      tATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTgg                                <  1:370720/149‑1 (MQ=255)
                                                                                                                      tATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTAATTGAATTTGATCTGCCGCTGCTgg                                <  2:373300/149‑1 (MQ=255)
                                                                                                                      tATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTAATTGAATTTGATCTGCCGCTGCTgg                                <  4:273026/149‑1 (MQ=255)
                                                                                                                        tCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGaa                              >  2:838478/1‑149 (MQ=255)
                                                                                                                        tCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGaa                              >  1:929984/1‑149 (MQ=255)
                                                                                                                        tCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGaa                              <  2:328888/149‑1 (MQ=255)
                                                                                                                          gAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAaga                            >  1:274972/1‑149 (MQ=255)
                                                                                                                            aTCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGaaa                                                                        >  4:248586/1‑103 (MQ=255)
                                                                                                                            aTCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGaaa                                                                        <  3:248586/103‑1 (MQ=255)
                                                                                                                              cTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGaaa                        <  3:210126/149‑1 (MQ=255)
                                                                                                                                gATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGc                      <  3:905026/149‑1 (MQ=255)
                                                                                                                                gATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGc                      <  3:413135/149‑1 (MQ=255)
                                                                                                                                 aTAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACGCCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGcc                     <  4:767298/149‑1 (MQ=255)
                                                                                                                                 aTAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAtt                                                 >  3:469454/1‑121 (MQ=255)
                                                                                                                                 aTAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAAtt                                                 <  4:469454/121‑1 (MQ=255)
                                                                                                                                    gcgtcgTTGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAg                  <  2:484278/149‑1 (MQ=255)
                                                                                                                                          ttGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTAcc            <  4:588804/149‑1 (MQ=255)
                                                                                                                                          ttGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTAcc            >  2:303776/1‑149 (MQ=255)
                                                                                                                                           tGAACTGTTAGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTAccc           <  1:300094/149‑1 (MQ=255)
                                                                                                                                             aaCTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTg         >  3:697350/1‑149 (MQ=255)
                                                                                                                                              aCTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGa        <  4:857346/149‑1 (MQ=255)
                                                                                                                                              aCTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGa        <  2:570359/149‑1 (MQ=255)
                                                                                                                                               cTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGAc       <  4:792604/149‑1 (MQ=255)
                                                                                                                                                 gTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTc     <  3:410536/149‑1 (MQ=255)
                                                                                                                                                    cGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCgg  <  3:92956/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGG  >  NC_000913/2535952‑2536248

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: