Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,535,976 G→T 100% G48V (GGT→GTT)  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,535,9760GT96.9% 224.8 / ‑0.1 64G48V (GGT→GTT) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (23/39);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACC  >  NC_000913/2535831‑2536121
                                                                                                                                                 |                                                                                                                                                 
aTCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGa                                                                                                                                                <  1:483854/149‑1 (MQ=255)
  cATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATg                                                                                                                                              <  3:402196/149‑1 (MQ=255)
  cATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATg                                                                                                                                              <  1:35521/149‑1 (MQ=255)
   aTGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATgg                                                                                                                                             >  1:168642/1‑149 (MQ=255)
   aTGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATgg                                                                                                                                             >  1:507518/1‑149 (MQ=255)
    tGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATAGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGt                                                                                                                                            <  1:394801/149‑1 (MQ=255)
       gTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGAt                                                                                                                                               <  3:229188/143‑1 (MQ=255)
       gTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGAt                                                                                                                                               >  4:229188/1‑143 (MQ=255)
         gtGTTCGATAAACTGAAAACTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGTCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATAGTTGTTGATGGTATTGc                                                                                                                                       <  2:153063/148‑1 (MQ=255)
           gTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTa                                                                                                                                     <  1:409558/149‑1 (MQ=255)
            ttCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATTCCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTAt                                                                                                                                    >  4:489073/1‑149 (MQ=255)
               gATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCaa                                                                                                                                 >  3:30026/1‑149 (MQ=255)
                aTAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCaaa                                                                                                                                <  1:562620/149‑1 (MQ=255)
                   aaCTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCa                                                                                                                             >  4:512073/1‑149 (MQ=255)
                     cTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAAc                                                                                                                           <  3:386470/149‑1 (MQ=255)
                       gAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACgg                                                                                                                         <  4:130413/149‑1 (MQ=255)
                       gAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAAATATTGAGATCATTGCTCCGCCCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACgg                                                                                                                         <  3:376973/149‑1 (MQ=255)
                        aaaTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACggg                                                                                                                        <  3:124329/149‑1 (MQ=255)
                          aTCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTa                                                                                                                      >  4:247810/1‑149 (MQ=255)
                           tctcTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTaa                                                                                                                     >  4:334251/1‑149 (MQ=255)
                                ggTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACaaaa                                                                                                                <  1:484922/149‑1 (MQ=255)
                                  tttCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGc                                                                                                                                       >  1:143036/1‑124 (MQ=255)
                                  tttCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGc                                                                                                                                       <  2:143036/124‑1 (MQ=255)
                                    tCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGt                                                                                                            <  4:385674/149‑1 (MQ=255)
                                       gacgacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTcgc                                                                                                         >  2:430986/1‑149 (MQ=255)
                                         cgacAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCgcgc                                                                                                       <  1:430986/149‑1 (MQ=255)
                                                 agGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGAc                                                                                               <  4:583500/149‑1 (MQ=255)
                                                 agGATACCGGAAATATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGAc                                                                                               <  3:489073/149‑1 (MQ=255)
                                                  ggATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACg                                                                                              >  4:127714/1‑149 (MQ=255)
                                                  ggATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACg                                                                                              >  1:289750/1‑149 (MQ=255)
                                                         ggAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCAt                                                                                       <  1:66380/149‑1 (MQ=255)
                                                           aaCTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTg                                                                                     <  4:65444/149‑1 (MQ=255)
                                                              tATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTa                                                                                  <  1:125269/149‑1 (MQ=255)
                                                                   agaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTgg                                                                                    <  3:66549/142‑1 (MQ=255)
                                                                   agaTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTgg                                                                                    >  4:66549/1‑142 (MQ=255)
                                                                     aTCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                           <  4:181649/149‑1 (MQ=255)
                                                                      tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                          >  4:460348/1‑149 (MQ=255)
                                                                      tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                          <  1:412874/149‑1 (MQ=255)
                                                                      tCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTACTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCttt                                                                          <  1:309251/149‑1 (MQ=255)
                                                                        aTTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGa                                                                        >  2:506960/1‑149 (MQ=255)
                                                                                gCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTTACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAaaccaacc                                                                >  1:103851/1‑149 (MQ=255)
                                                                                        gCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtc                                                        <  3:127714/149‑1 (MQ=255)
                                                                                        gCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtc                                                        <  3:460348/149‑1 (MQ=255)
                                                                                         cGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtct                                                       <  3:562898/149‑1 (MQ=255)
                                                                                          gagaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATtctc                                                      <  2:289750/149‑1 (MQ=255)
                                                                                            gaTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTa                                                    <  2:103517/149‑1 (MQ=255)
                                                                                             aTCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTAt                                                   <  4:30026/149‑1 (MQ=255)
                                                                                                gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                           <  2:30387/122‑1 (MQ=255)
                                                                                                gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCtt                                                                           >  1:30387/1‑122 (MQ=255)
                                                                                                gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                                >  4:346876/1‑149 (MQ=255)
                                                                                                gTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGa                                                >  2:565072/1‑149 (MQ=255)
                                                                                                       tCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTc                                                                                                           <  2:41596/83‑1 (MQ=255)
                                                                                                       tCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTc                                                                                                           >  1:41596/1‑83 (MQ=255)
                                                                                                       tCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGAt                                         <  2:534869/149‑1 (MQ=255)
                                                                                                        cGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATa                                        <  2:507518/149‑1 (MQ=255)
                                                                                                         gAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAg                                       >  1:152191/1‑149 (MQ=255)
                                                                                                          aaGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgc                                      >  1:359461/1‑149 (MQ=255)
                                                                                                             aCGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAgcggc                                   <  1:228578/149‑1 (MQ=255)
                                                                                                                tGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGtt                                <  1:108209/149‑1 (MQ=255)
                                                                                                                         tCGTTTTTGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTc                       <  3:334251/149‑1 (MQ=255)
                                                                                                                                tGCGGAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACt                <  4:259129/149‑1 (MQ=255)
                                                                                                                                    gAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg            >  1:454841/1‑149 (MQ=255)
                                                                                                                                    gAAAAAATCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCgg            >  1:395043/1‑149 (MQ=255)
                                                                                                                                        aaaTCGTTGTTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATc        <  2:519348/149‑1 (MQ=255)
                                                                                                                                              ttgttgATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACAcc  >  3:364895/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATCATGGGTTTGTTCGATAAACTGAAATCTCTGGTTTCCGACGACAAGAAGGATACCGGAACTATTGAGATCATTGCTCCGCTCTCTGGCGAGATCGTCAATATCGAAGACGTGCCGGATGTCGTTTTTGCGGAAAAAATCGTTGGTGATGGTATTGCTATCAAACCAACGGGTAACAAAATGGTCGCGCCAGTAGACGGCACCATTGGTAAAATCTTTGAAACCAACCACGCATTCTCTATCGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACC  >  NC_000913/2535831‑2536121

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