Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,040,663 A→C 5.6% intergenic (+27/+243) rrfA → / ← mobB 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,040,6630AC 5.6% 284.6 / 8.2 89intergenic (+27/+243)rrfA/mobB5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base C (2/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATC  >  NC_000913/4040515‑4040805
                                                                                                                                                    |                                                                                                                                              
tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCCTGCGAGAGTAGGGAACTGCCAGGCATCAAATAAATCGAAAGCCTCATTCCGaa                                                                                                                                                <  2:637123/149‑1 (MQ=11)
tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGaa                                                                                                                                                <  2:421032/149‑1 (MQ=34)
     cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTg                                                                                                                                           <  1:350420/149‑1 (MQ=35)
      tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTgg                                                                                                                                          <  4:165822/149‑1 (MQ=35)
        gcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGc                                                                                                                                        <  3:198409/149‑1 (MQ=35)
              gTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCg                                                                                                                                  >  4:808655/1‑149 (MQ=34)
                aGCGCGGTGGTCCCTCCTGACACCATGCCGAACTCAGAAGTGAAACGCTGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGtt                                                                                                                                <  1:243157/149‑1 (MQ=255)
                        ggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTCGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCTTTCGTTTTATttga                                                                                                                        <  1:667890‑M1/149‑5 (MQ=255)
                           cccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgttgt                                                                                                                     <  2:132796/149‑1 (MQ=35)
                            ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTgttgtt                                                                                                                    <  3:860849/149‑1 (MQ=35)
                             cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgttt                                                                                                                   <  2:677963/149‑1 (MQ=35)
                             cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgttt                                                                                                                   <  4:631904/149‑1 (MQ=35)
                              aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT                                                                                                                            <  3:724176‑M1/139‑1 (MQ=255)
                              aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT                                                                                                                            >  4:724176‑M1/1‑139 (MQ=255)
                                cTGAACCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTc                                                                                                                <  3:775798/149‑1 (MQ=21)
                                      ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGaa                                                                                                          <  1:61544/149‑1 (MQ=34)
                                      ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGaa                                                                                                          >  4:559461/1‑149 (MQ=35)
                                            cGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGctct                                                                                                    <  4:109696/149‑1 (MQ=35)
                                             gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTc                                                                                                                                   >  4:683673/1‑117 (MQ=38)
                                             gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTc                                                                                                                                   <  3:683673/117‑1 (MQ=38)
                                               aCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCt                                                                                                 <  4:681717/149‑1 (MQ=35)
                                                 taagaAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATttgatg                                                                                                                      <  3:110049‑M1/124‑7 (MQ=255)
                                                 tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                               >  2:906720/1‑149 (MQ=34)
                                                 tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGa                                                                                               <  3:407040/149‑1 (MQ=35)
                                                 tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATttgatg                                                                                                                      >  4:110049‑M1/1‑120 (MQ=255)
                                                  cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                              <  4:277377/149‑1 (MQ=35)
                                                   aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGc                                                                                                       >  4:640937/1‑139 (MQ=35)
                                                   aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGc                                                                                                       <  3:640937/139‑1 (MQ=37)
                                                    gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                                                            >  3:907890/1‑149 (MQ=34)
                                                     aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCTAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                                                           <  4:665600/149‑1 (MQ=25)
                                                     aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                                                           <  4:907890/149‑1 (MQ=35)
                                                       gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAGATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCTTTCGTTTTATt                                                                                                                           <  2:817728‑M1/115‑2 (MQ=255)
                                                       gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAGATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCTTTCGTTTTATt                                                                                                                           >  1:817728‑M1/1‑114 (MQ=255)
                                                        tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                                                           <  1:873786/146‑1 (MQ=34)
                                                        tGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                                                           >  2:873786/1‑146 (MQ=35)
                                                          aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaa                                                                                      >  2:302215/1‑149 (MQ=35)
                                                           aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaaa                                                                                     <  3:308154/149‑1 (MQ=35)
                                                                 gTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                                                               <  3:808655/149‑1 (MQ=35)
                                                                  tAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgctc                                                                              >  1:264053/1‑147 (MQ=25)
                                                                         gATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCgg                                                                       >  1:812568/1‑149 (MQ=35)
                                                                            ggTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAtt                                                                    <  2:812568/149‑1 (MQ=34)
                                                                             gTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAttt                                                                   >  2:783468/1‑149 (MQ=34)
                                                                                    ggggTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGt                                                                                                                                 >  4:442224/1‑80 (MQ=255)
                                                                                    ggggTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGt                                                                                                                                 <  3:442224/80‑1 (MQ=40)
                                                                                     gggTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                               >  3:670683/1‑145 (MQ=34)
                                                                                     gggTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                               <  4:670683/145‑1 (MQ=34)
                                                                                     gggTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTg                                                           <  3:559461/149‑1 (MQ=35)
                                                                                         ctcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGaa                                                       <  4:830836/149‑1 (MQ=35)
                                                                                         ctcctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGaa                                                       <  2:175926/149‑1 (MQ=35)
                                                                                            ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGc                                                          <  1:62703/143‑1 (MQ=34)
                                                                                            ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGc                                                          >  2:62703/1‑143 (MQ=34)
                                                                                            ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                    >  4:546944/1‑149 (MQ=33)
                                                                                            ctcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                    >  1:800703/1‑149 (MQ=34)
                                                                                             tcATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCTGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGATCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCaa                                                   <  4:582063/149‑1 (MQ=18)
                                                                                                  cGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCAAAGCAACGGcc                                              <  1:693782/149‑1 (MQ=17)
                                                                                                   gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                             >  2:893605/1‑149 (MQ=32)
                                                                                                      agTAGGGAACTGCCAGGCATCAAATAAAACGAAGGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGa                                          <  1:893605/149‑1 (MQ=17)
                                                                                                         aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGcctggcagt                                       >  1:332440/1‑142 (MQ=17)
                                                                                                         aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGCCGGAAGACTGGGCCTTTCGTTTTATCTATTGTTTGTCGGTGAACGCTCTCCTGAGTAGGTCAAATACGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAggg                                       >  3:776339/1‑149 (MQ=255)
                                                                                                              aCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg                                  <  1:896545/149‑1 (MQ=32)
                                                                                                                  ccAGGCATCAAATAAAACGAAAGGCACAGGCGGAAGGCTGGGCTTTTCGTTTTATCTGTGGTTTGTCGGCGAACGCTCTCTTGAGTAGGACAACTCCGCCGGGAGCGGATTTGAACGTTGCGAAACACCGGCCCGGAGGGTGGCGGCCa                              <  4:776339/149‑1 (MQ=11)
                                                                                                                         tCAAATAAATCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCGGGAGGGTGGCGGGCAGGAcgcc                       >  1:802715/1‑149 (MQ=17)
                                                                                                                           aaaTAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccac                     >  1:172811‑M1/1‑46 (MQ=255)
                                                                                                                            aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                    >  4:394332/1‑149 (MQ=32)
                                                                                                                                 aacaaaaggctCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggg                                    <  4:286644‑M1/124‑89 (MQ=255)
                                                                                                                                 aaCAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggg                                    >  3:286644‑M1/1‑40 (MQ=255)
                                                                                                                                  aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgcta              >  3:115762‑M1/1‑39 (MQ=255)
                                                                                                                                  aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgcta              >  3:138307‑M1/1‑39 (MQ=255)
                                                                                                                                   caaaaggCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctac             >  1:915517‑M1/3‑38 (MQ=255)
                                                                                                                                    gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt            <  3:546944/149‑1 (MQ=32)
                                                                                                                                    gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt            >  1:553614/1‑149 (MQ=32)
                                                                                                                                    gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt            >  1:173346/1‑149 (MQ=32)
                                                                                                                                    aaaagGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctact            >  1:245092‑M1/2‑37 (MQ=255)
                                                                                                                                      aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgc          >  3:401982‑M1/1‑35 (MQ=255)
                                                                                                                                       aGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc         >  3:437560‑M1/1‑34 (MQ=255)
                                                                                                                                       aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc         >  4:557897‑M1/1‑34 (MQ=255)
                                                                                                                                       aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgcc         >  4:868984‑M1/1‑34 (MQ=255)
                                                                                                                                       aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgttt                                                         >  1:559288‑M1/1‑34 (MQ=255)
                                                                                                                                       aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgttt                                                         <  2:559288‑M1/101‑68 (MQ=255)
                                                                                                                                        ggCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgttt                                                         >  1:581489‑M1/1‑33 (MQ=255)
                                                                                                                                        ggCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgttt                                                         <  2:581489‑M1/100‑68 (MQ=255)
                                                                                                                                        ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCa        <  2:553614/149‑1 (MQ=32)
                                                                                                                                        ggCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccg        >  1:57058‑M1/1‑33 (MQ=255)
                                                                                                                                        ggCTCAGTCGAAAGACTGGGCCTTTCGGTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccg        >  4:444095‑M1/1‑33 (MQ=255)
                                                                                                                                         gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAg       >  4:558706/1‑149 (MQ=32)
                                                                                                                                           tCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgcca     >  3:136958‑M1/1‑30 (MQ=255)
                                                                                                                                           tCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgcca     <  3:557897‑M1/149‑120 (MQ=255)
                                                                                                                                           tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc     >  4:615937/1‑149 (MQ=32)
                                                                                                                                            cAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgccag    <  4:136958‑M1/149‑121 (MQ=255)
                                                                                                                                              gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc  >  1:937352/1‑149 (MQ=32)
                                                                                                                                              gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc  >  2:191939/1‑149 (MQ=32)
                                                                                                                                                    |                                                                                                                                              
TTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATC  >  NC_000913/4040515‑4040805

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