Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,158,139 A→T 8.2% K91* (AAA→TAA) ‡ ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,1390AT 8.2% 334.3 / 17.4 97K91* (AAA→TAA) ‡ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base T (4/4);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

Only 100 of 101 total aligned reads displayed.

TGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAA  >  NC_000913/1157997‑1158284
                                                                                                                                              |                                                                                                                                                 
tGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAAcc                                                                                                                                              >  1:919766/1‑148 (MQ=255)
 gTATCGCAAGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCa                                                                                                                                             <  2:919766/148‑1 (MQ=255)
   aTCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTTTAACCATg                                                                                                                                           >  1:452582/1‑148 (MQ=255)
    tCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTTTAACCATga                                                                                                                                          <  2:452582/148‑2 (MQ=255)
                     aGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTaaa                                                                                                                                                 <  1:518535/124‑1 (MQ=255)
                     aGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTaaa                                                                                                                                                 >  2:518535/1‑124 (MQ=255)
                       cAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCAt                                                                                                                                            <  3:1136890/127‑1 (MQ=255)
                       cAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCAt                                                                                                                                            >  4:1136890/1‑127 (MQ=255)
                               tCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACt                                                                                                                          >  1:489646/1‑137 (MQ=255)
                               tCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACt                                                                                                                          <  2:489646/137‑1 (MQ=255)
                                ccGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTTTAACCATga                                                                                                                                          >  2:301760/1‑119 (MQ=255)
                                ccGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTTTAACCATga                                                                                                                                          <  1:301760/120‑2 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTaaa                                                                                                                                                 <  4:1105063/110‑1 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTaaa                                                                                                                                                 >  3:1105063/1‑110 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTgcg                                                                                                                               >  3:564654/1‑128 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTgcg                                                                                                                               <  4:564654/128‑1 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTa                                                                                                                      <  4:111022/137‑1 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTa                                                                                                                      >  3:111022/1‑137 (MQ=255)
                                   tCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTAcc                                                                                                           >  2:93454/1‑148 (MQ=255)
                                           aaCATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATg                                                                                                                                           >  4:459213/1‑108 (MQ=255)
                                           aaCATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATg                                                                                                                                           <  3:459213/108‑1 (MQ=255)
                                               tGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAAt                                                                                               <  1:575335/148‑1 (MQ=255)
                                                 ccACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCg                                                                                             >  1:616809/1‑148 (MQ=255)
                                                   aCTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCa                                                                                                                            >  1:721546/1‑115 (MQ=255)
                                                   aCTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCa                                                                                                                            >  1:412597/1‑115 (MQ=255)
                                                   aCTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCa                                                                                                                            <  2:412597/115‑1 (MQ=255)
                                                   aCTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTAACAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCa                                                                                                                            <  2:721546/115‑1 (MQ=255)
                                                    cTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCct                                                                                          >  3:695590/1‑148 (MQ=255)
                                                     tGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTgaaga                                                                                                  >  2:1084174/1‑139 (MQ=255)
                                                     tGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTgaaga                                                                                                  <  1:1084174/139‑1 (MQ=255)
                                                      gATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCAt                                                                                                                >  2:30875/1‑124 (MQ=255)
                                                      gATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCAt                                                                                                                <  1:30875/124‑1 (MQ=255)
                                                         ttttGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCATCTaaa                                                                                     <  2:616809/148‑1 (MQ=255)
                                                                aTCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCTCCAGTTGTTGCCTATGGCATCATGGTTAAAGCCATGGCCGTGGTTGCGCCACTGGTACTGc                                                                                                                  <  1:51495/112‑1 (MQ=255)
                                                                aTCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCa                                                                                                                                             <  2:956649/85‑1 (MQ=255)
                                                                aTCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGc                                                                                                                  >  2:51495/1‑112 (MQ=255)
                                                                aTCAGTGTCGCCCTCGGCTTTACTAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCa                                                                                                                                             >  1:956649/1‑85 (MQ=255)
                                                                          cccTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCgg                                                                           <  4:962996/141‑1 (MQ=255)
                                                                          cccTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCgg                                                                           >  3:962996/1‑141 (MQ=255)
                                                                          cccTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATAGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTgg                                                                    >  3:1114651/1‑148 (MQ=255)
                                                                            cTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCa                                                                                                                            >  2:595246/1‑90 (MQ=255)
                                                                            cTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCa                                                                                                                            <  1:595246/90‑1 (MQ=255)
                                                                             tCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCg                                                                                             <  2:136579/120‑1 (MQ=255)
                                                                             tCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCg                                                                                             >  1:136579/1‑120 (MQ=255)
                                                                                   ttACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTgaagaa                                                                                                 <  3:1094530/110‑1 (MQ=255)
                                                                                   ttACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTgaagaa                                                                                                 >  4:1094530/1‑110 (MQ=255)
                                                                                   ttACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCgg                                                           <  2:113918/148‑1 (MQ=255)
                                                                                   ttACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCgg                                                           >  1:113918/1‑148 (MQ=255)
                                                                                    tACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGa                                                          <  1:93454/148‑1 (MQ=255)
                                                                                       cAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACt                                                                      >  1:708817/1‑133 (MQ=255)
                                                                                       cAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACt                                                                      <  2:708817/133‑1 (MQ=255)
                                                                                        aaTAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGa                                                      <  3:964517/148‑1 (MQ=255)
                                                                                         aTAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGc                                                                                                                  >  2:1042319/1‑87 (MQ=255)
                                                                                         aTAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGc                                                                                                                  <  1:1042319/87‑1 (MQ=255)
                                                                                           aaCGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGAt                                                     <  1:842867/146‑1 (MQ=255)
                                                                                           aaCGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGAt                                                     >  2:842867/1‑146 (MQ=255)
                                                                                                  gCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCgg                                            <  4:695590/148‑1 (MQ=255)
                                                                                                   cGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTAtctc                                               <  1:238252/144‑1 (MQ=255)
                                                                                                   cGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTAtctc                                               >  2:238252/1‑144 (MQ=255)
                                                                                                         cgcgCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATc                                     >  4:1136157/1‑148 (MQ=255)
                                                                                                                 ccGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATc                                                                                              <  2:655294/83‑1 (MQ=255)
                                                                                                                 ccGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATc                                                                                              >  1:655294/1‑83 (MQ=255)
                                                                                                                  cGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTAc                            <  3:1136157/148‑1 (MQ=255)
                                                                                                                   gCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGc                                                                                            >  2:91059/1‑83 (MQ=255)
                                                                                                                   gCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGc                                                                                            <  1:91059/83‑1 (MQ=255)
                                                                                                                          ttgCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTgg                                                                              >  1:651316/1‑90 (MQ=255)
                                                                                                                          ttgCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTgg                                                                              <  2:651316/90‑1 (MQ=255)
                                                                                                                               tATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAAt                                                                                               <  4:508037/68‑1 (MQ=255)
                                                                                                                               tATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAAt                                                                                               >  3:508037/1‑68 (MQ=255)
                                                                                                                               tATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAAt                                                                                               >  2:349120/1‑68 (MQ=255)
                                                                                                                               tATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAAt                                                                                               <  1:349120/68‑1 (MQ=255)
                                                                                                                               tATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCg                                                                                             >  4:449620/1‑70 (MQ=255)
                                                                                                                               tATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCg                                                                                             <  3:449620/70‑1 (MQ=255)
                                                                                                                                  ggCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCt            <  4:1114651/148‑1 (MQ=255)
                                                                                                                                    catcatGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGAt                                      >  1:171473/1‑120 (MQ=255)
                                                                                                                                    catcatGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGAt                                      <  2:171473/120‑1 (MQ=255)
                                                                                                                                    catcatGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACa                           <  4:419932/131‑1 (MQ=255)
                                                                                                                                    catcatGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACGCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACa                           >  3:419932/1‑131 (MQ=255)
                                                                                                                                     atcatGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGc                                         <  1:239237/116‑1 (MQ=255)
                                                                                                                                     atcatGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGc                                         >  2:239237/1‑116 (MQ=255)
                                                                                                                                       catGGTTTTAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCt                                                                                >  2:850549/1‑75 (MQ=255)
                                                                                                                                       catGGTTTTAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCt                                                                                <  1:850549/75‑1 (MQ=255)
                                                                                                                                       catGGTTTTAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTa                     <  4:641226/134‑1 (MQ=255)
                                                                                                                                       catGGTTTTAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTa                     >  3:641226/1‑134 (MQ=255)
                                                                                                                                       catGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGa                                                          <  2:718291/97‑1 (MQ=255)
                                                                                                                                       catGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGa                                                          >  1:718291/1‑97 (MQ=255)
                                                                                                                                       catGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGAt                                                     >  3:766262/1‑102 (MQ=255)
                                                                                                                                       catGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGAt                                                     <  4:766262/102‑1 (MQ=255)
                                                                                                                                         tGGTTAAAACCATGGCCGTGGTTGCGCCAATGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGt                <  3:487953/137‑1 (MQ=255)
                                                                                                                                         tGGTTAAAACCATGGCCGTGGTTGCGCCAATGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGt                >  4:487953/1‑137 (MQ=255)
                                                                                                                                           gTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTa   >  1:373815/1‑148 (MQ=255)
                                                                                                                                            ttAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTga                                                                                                     <  2:93526/49‑1 (MQ=255)
                                                                                                                                            ttAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTga                                                                                                     >  1:93526/1‑49 (MQ=255)
                                                                                                                                            ttAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTaa  >  3:775388/1‑148 (MQ=255)
                                                                                                                                             tAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTAt                                                  >  3:721712/1‑99 (MQ=255)
                                                                                                                                             tAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTAc                            >  2:1073619/1‑121 (MQ=255)
                                                                                                                                             tAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTAc                            <  1:1073619/121‑1 (MQ=255)
                                                                                                                                             tAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTaa                    >  1:151622/1‑129 (MQ=255)
                                                                                                                                             tAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTaa                    <  2:151622/129‑1 (MQ=255)
                                                                                                                                             tAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATCCACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATCAt                                                  <  4:721712/99‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGCAAACATGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACCAATAACGATGGCGTATCCGCGCTGGCCGCAGTTGTTGCCTATGGCATCATGGTTAAAACCATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTACCGTATTAA  >  NC_000913/1157997‑1158284

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