Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,554 G→A 5.1% L171L (CTG→CTA cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,5540GA 5.1% 398.1 / 6.5 117L171L (CTG→CTAcydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base G (56/55);  new base A (3/3);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 124 total aligned reads displayed.

CATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTAT  >  NC_000913/773407‑773702
                                                                                                                                                   |                                                                                                                                                    
cATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTgg                                                                                                                                                     >  2:139209/1‑149 (MQ=255)
 aTCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGc                                                                                                                                                    >  2:72493/1‑149 (MQ=255)
   cTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAg                                                                                                                                                  >  1:506952/1‑149 (MQ=255)
     tCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAgg                                                                                                                                                 >  3:851385/1‑148 (MQ=255)
        ttGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                                                                                                                                              <  4:726278/148‑1 (MQ=255)
        ttGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgt                                                                                                                                              <  4:851385/148‑1 (MQ=255)
         tGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtg                                                                                                                                             >  3:758744/1‑148 (MQ=255)
         tGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCCAGCAGGccga                                                                                                                                             <  3:255777/148‑4 (MQ=255)
          ggTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggt                                                                                                                                           <  1:427542/149‑1 (MQ=255)
            tAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAGTATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCgtggtg                                                                                                                                          <  3:62767/148‑1 (MQ=255)
            tAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGa                                                                                                                                         <  1:174566/149‑1 (MQ=255)
            tAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGa                                                                                                                                         <  1:902326/149‑1 (MQ=255)
               cTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGc                                                                                                                                       >  4:158805/1‑148 (MQ=255)
               cTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCg                                                                                                                                      >  1:68790/1‑149 (MQ=255)
                ttCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                                                                                                                                     >  1:345227/1‑149 (MQ=255)
                 tCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                                                                                                                                     <  3:608890/148‑1 (MQ=255)
                 tCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGt                                                                                                                                     >  4:679105/1‑148 (MQ=255)
                  cGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAg                                                                                                                                   <  2:506952/149‑1 (MQ=255)
                  cGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAg                                                                                                                                   <  2:578731/149‑1 (MQ=255)
                  cGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAg                                                                                                                                   <  1:489488/149‑1 (MQ=255)
                      ccgccgCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgat                                                                                                                               <  1:445205/149‑1 (MQ=255)
                         ccgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatgat                                                                                                                            <  2:898369/149‑1 (MQ=255)
                          cgcTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGgatgat                                                                                                                            <  4:758744/148‑1 (MQ=255)
                               gTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTAc                                                                                                                      <  1:2137/149‑1 (MQ=255)
                                  aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                   <  1:1012946/149‑1 (MQ=255)
                                  aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                   >  2:922050/1‑149 (MQ=255)
                                  aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                   <  2:761245/149‑1 (MQ=255)
                                  aTTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCa                                                                                                                   >  2:618537/1‑149 (MQ=255)
                                    tGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGCACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAg                                                                                                                  <  4:254068/148‑1 (MQ=255)
                                      gtgtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAggg                                                                                                                <  4:284364/148‑1 (MQ=255)
                                      gtgtAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGgc                                                                                                               >  2:278283/1‑149 (MQ=255)
                                      gtgtAGCGTTAGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTCAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAggg                                                                                                                >  3:697489/1‑148 (MQ=255)
                                        gtAGCGTTAGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGTATATTTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACACGTTCGGCCTGCTGGCAGGCGTGATGAGCGTAGGGATGATCATTACTCAGGgcg                                                                                                              >  4:781620/1‑148 (MQ=255)
                                                 ggTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATc                                                                                                     <  4:823627/148‑1 (MQ=255)
                                                     aCCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTCACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCACTACTCAGGGCGCAACCTATCTGCa                                                                                                 >  3:611407/1‑148 (MQ=255)
                                                      ccTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCaaa                                                                                               <  1:187881/149‑1 (MQ=255)
                                                      ccTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCaa                                                                                                >  3:349896/1‑148 (MQ=255)
                                                        tGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAAt                                                                                              >  4:947955/1‑148 (MQ=255)
                                                         gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGc                                                                                            <  1:960743/149‑1 (MQ=255)
                                                         gTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGc                                                                                            >  2:25669/1‑149 (MQ=255)
                                                            gtaGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTa                                                                                         >  2:122161/3‑149 (MQ=255)
                                                               gggCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTAt                                                                                                      <  2:352015/133‑1 (MQ=255)
                                                               gggCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTAt                                                                                                      >  1:352015/1‑133 (MQ=255)
                                                               gggCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCg                                                                                      >  2:962350/1‑149 (MQ=255)
                                                                  cGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGAAGGCGTGGTGAGCGTAGGGATGATTACTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTg                                                                                    <  4:611407/148‑1 (MQ=255)
                                                                      ccGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCg                                                                                <  3:1029805/148‑1 (MQ=255)
                                                                      ccGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGa                                                                               >  1:267253/1‑149 (MQ=255)
                                                                      ccGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGa                                                                               <  2:853000/149‑1 (MQ=255)
                                                                              cGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGaa                                                                              <  4:160784/142‑1 (MQ=255)
                                                                              cGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGaa                                                                              >  3:160784/1‑142 (MQ=255)
                                                                               gTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCt                                                                      <  2:505429/149‑1 (MQ=255)
                                                                                   atgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATc                                                                                                     <  1:945624/114‑1 (MQ=255)
                                                                                   atgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTCCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATc                                                                                                     >  2:945624/1‑114 (MQ=255)
                                                                                    tgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGt                                                                  <  4:3773/148‑1 (MQ=255)
                                                                                    tgaATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTa                                                                 <  2:825472/149‑1 (MQ=255)
                                                                                      aaTATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTg                                                                                    <  4:740333/128‑1 (MQ=255)
                                                                                      aaTATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTg                                                                                    >  3:740333/1‑128 (MQ=255)
                                                                                      aaTATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTAc                                                                >  3:859484/1‑148 (MQ=255)
                                                                                       atatCTTCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTAccc                                                              >  2:264119/1‑149 (MQ=255)
                                                                                       atatCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTAccc                                                              <  1:278283/149‑1 (MQ=255)
                                                                                              cGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCaa                                                        >  4:893776/1‑148 (MQ=255)
                                                                                                     actacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCa                                                 >  3:943860/1‑148 (MQ=255)
                                                                                                      ctacACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAg                                                <  3:593306/148‑1 (MQ=255)
                                                                                                        acacCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTg                                             >  2:290196/1‑149 (MQ=255)
                                                                                                         cacCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTg                                             >  4:636486/1‑148 (MQ=255)
                                                                                                            cGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCt                                          >  3:305952/1‑148 (MQ=255)
                                                                                                               tAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCAc                                                                        <  4:656098/115‑1 (MQ=255)
                                                                                                                 aCTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTg                                    <  1:1062723/149‑1 (MQ=255)
                                                                                                                   ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTgg                                   >  4:1000145/1‑148 (MQ=255)
                                                                                                                   ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                  <  1:31644/149‑1 (MQ=255)
                                                                                                                   ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                  >  2:218030/1‑149 (MQ=255)
                                                                                                                   ttcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                  >  1:446587/1‑149 (MQ=255)
                                                                                                                    tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCt                                     <  3:195167/145‑1 (MQ=255)
                                                                                                                    tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGt                                  >  3:511281/1‑148 (MQ=255)
                                                                                                                    tcttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTg                                 <  1:27711/149‑1 (MQ=255)
                                                                                                                      ttcCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGac                               >  1:985719/1‑149 (MQ=255)
                                                                                                                          aGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACt                            <  3:430384/148‑1 (MQ=255)
                                                                                                                          aGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTg                           <  1:310861/149‑1 (MQ=255)
                                                                                                                           gTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                          >  1:543135/1‑149 (MQ=255)
                                                                                                                            ttGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTgg                          <  3:198268/148‑1 (MQ=255)
                                                                                                                              gCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTc                        <  4:943860/148‑1 (MQ=255)
                                                                                                                                ttAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCt                                                                                                    >  2:622279/1‑70 (MQ=255)
                                                                                                                                ttAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCt                                                                                                    <  1:622279/70‑1 (MQ=255)
                                                                                                                                  aaCCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGtt                    >  4:157013/1‑148 (MQ=255)
                                                                                                                                   aCCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc                  >  1:974454/1‑149 (MQ=255)
                                                                                                                                    cccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGc                                           >  4:801308/1‑123 (MQ=255)
                                                                                                                                    cccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc                  >  3:811632/1‑148 (MQ=255)
                                                                                                                                    cccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCg                 >  1:1030380/1‑149 (MQ=255)
                                                                                                                                     ccGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCg                 >  4:166841/1‑148 (MQ=255)
                                                                                                                                       gTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAc              >  1:724836/1‑149 (MQ=255)
                                                                                                                                       accacgCCTGCGAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAc              >  1:561769/6‑149 (MQ=255)
                                                                                                                                       accacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc                  >  1:43039/6‑145 (MQ=255)
                                                                                                                                       accacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTc                  <  2:43039/140‑1 (MQ=255)
                                                                                                                                        ttCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCAct             >  2:151939/1‑149 (MQ=255)
                                                                                                                                         cacgCCTGCCAGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCCGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActg            >  1:522681/4‑149 (MQ=255)
                                                                                                                                             ccTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActggct         >  4:352692/1‑148 (MQ=255)
                                                                                                                                               tgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTgg                                            >  4:56713/1‑111 (MQ=255)
                                                                                                                                               tgctgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCActggctgg       >  4:590591/1‑148 (MQ=255)
                                                                                                                                                  tgGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGt    <  3:721412/148‑1 (MQ=255)
                                                                                                                                                   ggCAGGCGTGGTGAGCGTAGGGGTGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTAt  >  2:592920/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATCTTCATTGGTAGCTTCGTTCCGCCGCTGGTAATTGGTGTAGCGTTCGGTAACCTGTTGCAGGGCGTACCGTTCAACGTTGATGAATATCTGCGTCTGTACTACACCGGTAACTTCTTCCAGTTGCTTAACCCGTTCGGCCTGCTGGCAGGCGTGGTGAGCGTAGGGATGATCATTACTCAGGGCGCAACCTATCTGCAAATGCGTACCGTGGGCGAACTGCACCTGCGTACCCGTGCAACGGCTCAGGTGGCTGCGCTGGTGACACTGGTCTGTTTCGCACTGGCTGGCGTAT  >  NC_000913/773407‑773702

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