Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,171 T→A 9.0% noncoding (170/184 nt) glmY ← sRNA activator of glmS mRNA, glmZ processing antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1710TA 9.0% 423.1 / 24.8 122noncoding (170/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base T (58/53);  new base A (5/6);  total (63/59)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 122 total aligned reads displayed.

AGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCT  >  NC_000913/2691031‑2691315
                                                                                                                                             |                                                                                                                                                
agAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGc                                                                                                                                          <  1:334676/149‑1 (MQ=255)
agAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGAGGCTAACCGAGGGGAAGTTCAGATACAACAAAGACGGG‑AATTACCCGGc                                                                                                                                          <  2:807335/149‑1 (MQ=255)
 gAATCAGCAAAAATGCCAGCATTACCAGTTTTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                         <  1:1452251/149‑1 (MQ=255)
   cccatctaaaaTGCCAGCATTACCAGTTATCGTAATGAGCGGGGAAAAACGAGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCTGGCttt                                                                                                                                       <  2:822719/142‑1 (MQ=255)
   aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCttt                                                                                                                                       <  1:912928/149‑1 (MQ=255)
       gCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTa                                                                                                                                   >  2:568938/1‑149 (MQ=255)
         aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCttt                                                                                                                                       >  2:1362883/1‑143 (MQ=255)
                 cAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                         >  4:646359/1‑133 (MQ=255)
                 cAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                         <  3:646359/133‑1 (MQ=255)
                            gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                 <  1:1149138/130‑1 (MQ=255)
                                       gAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCg                                                                                                   >  3:146409/1‑149 (MQ=255)
                                        aGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                    >  1:522445/1‑147 (MQ=255)
                                        aGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                    <  2:522445/147‑1 (MQ=255)
                                          cGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCCGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACt                                                                                                             <  2:1531399/136‑1 (MQ=255)
                                          cGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATAGAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGc                                                                                                                   >  1:544193/1‑130 (MQ=255)
                                          cGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGc                                                                                                                   <  2:544193/130‑1 (MQ=255)
                                          cGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGccc                                                                                                >  2:57539/1‑149 (MQ=255)
                                          cGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATCCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACt                                                                                                             >  1:1531399/1‑136 (MQ=255)
                                           ggggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                 >  3:1078451/1‑131 (MQ=255)
                                            gggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTg                                                                                                                                      >  3:182585/1‑109 (MQ=255)
                                            gggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                                  >  4:1887254/1‑113 (MQ=255)
                                            gggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGa                                                                                              >  1:1621117/1‑149 (MQ=255)
                                            gggAAAAACGGGCCAGCGGTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTg                                                                                                                                      <  4:182585/109‑1 (MQ=39)
                                                 aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCtt                                                                                                                                        >  3:656611/1‑102 (MQ=255)
                                                 aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                         >  3:330904/1‑101 (MQ=255)
                                                 aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                         <  4:330904/101‑1 (MQ=255)
                                                 aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                         <  2:1621117/149‑1 (MQ=255)
                                                   aCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGa                                                                                       <  1:57539/149‑1 (MQ=255)
                                                     gggCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                         >  2:320785/1‑145 (MQ=255)
                                                      ggCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGatat                                                                                    <  1:568938/149‑1 (MQ=255)
                                                       gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCt                                                                                                                                         <  4:450961/95‑1 (MQ=255)
                                                           gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                         >  4:752206/1‑139 (MQ=255)
                                                              tttCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGc                                                                            >  4:14819/1‑149 (MQ=255)
                                                                 cAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                         <  3:14819/149‑1 (MQ=255)
                                                                  aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTg                                                                                                                                      >  1:1634009/1‑87 (MQ=255)
                                                                  aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTg                                                                                                                                      <  2:1634009/87‑1 (MQ=255)
                                                                  aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                          <  3:1004005/131‑1 (MQ=255)
                                                                  aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                          <  3:514630/131‑1 (MQ=255)
                                                                  aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                          >  4:1004005/1‑131 (MQ=255)
                                                                  aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGt                                                                                          >  4:514630/1‑131 (MQ=255)
                                                                        ttACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGa                                                                  >  2:339881/1‑149 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTg                                                                                                                                      <  4:1604301/77‑1 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCt                                                                                                                 <  1:959883/98‑1 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                    >  3:186063/1‑111 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                    <  4:186063/111‑1 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGa                                                                                       <  2:110134/124‑1 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGa                                                                                       >  1:110134/1‑124 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGata                                                                       >  3:39147/1‑140 (MQ=255)
                                                                            tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGata                                                                       <  4:39147/140‑1 (MQ=255)
                                                                             ctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACt                                                                                                             <  1:1002483/101‑1 (MQ=255)
                                                                             ctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACt                                                                                                             >  2:1002483/1‑101 (MQ=255)
                                                                             ctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTg                                                             <  2:1560940/149‑1 (MQ=255)
                                                                              tcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTgg                                                            >  1:191548/1‑149 (MQ=255)
                                                                                gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGc                                                                                                                                          >  4:1663732/1‑69 (MQ=255)
                                                                                gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGc                                                                                                                                          <  3:1663732/69‑1 (MQ=255)
                                                                                gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCt                                                                           <  2:866010/132‑1 (MQ=255)
                                                                                gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCt                                                                           >  1:866010/1‑132 (MQ=255)
                                                                                gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                          >  1:1229207/1‑149 (MQ=255)
                                                                                  cAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACgg                                                        >  4:678810/1‑149 (MQ=255)
                                                                                   aGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggc                                                       >  3:195660/1‑149 (MQ=255)
                                                                                    gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                      <  1:339881/149‑1 (MQ=255)
                                                                                         gAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCAc                                                 <  4:195660/149‑1 (MQ=255)
                                                                                           aTAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAg                                                                                  >  3:1767503/1‑114 (MQ=255)
                                                                                               ccTGATGCTAACCGAGGGGAAGTTTAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                           <  2:1229207/149‑1 (MQ=255)
                                                                                                cTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGat                                                                        >  4:1470380/1‑119 (MQ=255)
                                                                                                atGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGat                                                                        <  3:1470380/118‑1 (MQ=255)
                                                                                                    tGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                 <  3:1253852/58‑1 (MQ=255)
                                                                                                    tGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                                 >  4:1253852/1‑58 (MQ=255)
                                                                                                     gCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGt                                     <  2:191548/149‑1 (MQ=255)
                                                                                                         aCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                                  <  4:764029/52‑1 (MQ=255)
                                                                                                         aCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                                  >  3:764029/1‑52 (MQ=255)
                                                                                                              ggggAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCtt                                                                                                                                        >  1:399152/1‑41 (MQ=255)
                                                                                                                  aaGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGa                                                           >  3:30817/1‑114 (MQ=255)
                                                                                                                    gTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCAc                                                 >  2:567796/1‑122 (MQ=255)
                                                                                                                     ttCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                         >  3:1446493/1‑129 (MQ=255)
                                                                                                                     ttCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                         <  4:1446493/129‑1 (MQ=255)
                                                                                                                       cAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAg                   <  3:678810/149‑1 (MQ=255)
                                                                                                                         gATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                               <  1:1806349/55‑1 (MQ=255)
                                                                                                                         gATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                               >  2:1806349/1‑55 (MQ=255)
                                                                                                                          aTACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCa                                   >  2:1445028/1‑130 (MQ=255)
                                                                                                                          aTACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCa                                   <  1:1445028/130‑1 (MQ=255)
                                                                                                                           tACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                         >  4:1035966/1‑91 (MQ=255)
                                                                                                                           tACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                         <  3:1035966/91‑1 (MQ=255)
                                                                                                                           tACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGa                                                                         <  3:895987/91‑1 (MQ=255)
                                                                                                                              aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                          >  4:1331364/1‑103 (MQ=255)
                                                                                                                              aacaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                          <  1:1491288/103‑1 (MQ=255)
                                                                                                                               acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc           >  1:1410887/1‑149 (MQ=255)
                                                                                                                               acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc           >  3:222161/1‑149 (MQ=255)
                                                                                                                                  aaGCCGGGTAATT‑CCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa        >  4:1148074/1‑149 (MQ=255)
                                                                                                                                  aaGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGt              <  1:790220/143‑1 (MQ=255)
                                                                                                                                  aaGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGt              >  2:790220/1‑143 (MQ=255)
                                                                                                                                  aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa        >  2:1388333/1‑149 (MQ=255)
                                                                                                                                  aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa        >  2:1555497/1‑149 (MQ=255)
                                                                                                                                   aGCCGGGTAATT‑CCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                               <  3:1004995/93‑1 (MQ=39)
                                                                                                                                   aGCCGGGTAATT‑CCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                               >  4:1004995/1‑93 (MQ=39)
                                                                                                                                   aGCCGGGTAATT‑CCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                           >  3:1487131/1‑113 (MQ=255)
                                                                                                                                   aGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                           <  4:1487131/113‑1 (MQ=255)
                                                                                                                                   aGCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg       >  1:1050458/1‑149 (MQ=255)
                                                                                                                                        gg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                       >  1:717294/1‑112 (MQ=255)
                                                                                                                                        gg‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCt  >  4:799998/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                
AGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCT  >  NC_000913/2691031‑2691315

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