Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,610,720 T→A 6.1% L436Q (CTG→CAG)  prfC → peptide chain release factor RF‑3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,610,7200TA 6.1% 283.0 / 7.4 82L436Q (CTG→CAG) prfCpeptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (29/48);  new base A (2/3);  total (31/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

ACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAGTTC  >  NC_000913/4610575‑4610868
                                                                                                                                                 |                                                                                                                                                    
aCCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGAt                                                                                                                                                   >  4:157437/1‑149 (MQ=255)
   ggTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAACGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACACCTTTCCGTAGAGGGCGCGGTGCAGGTGTTCCGTCCACTCTCCAACAACGATCTGATCGt                                                                                                                                                <  4:637738/149‑1 (MQ=255)
     tATTCCGAACTTCTCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCAAAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTg                                                                                                                                              <  2:34237/149‑1 (MQ=255)
       ttCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGTCGCGGTGCAGGTGTTCCGTCCAATCTCCACCAACGATCAGATCGTTgtt                                                                                                                                            >  1:82166/1‑147 (MQ=255)
        tCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGttgttg                                                                                                                                           <  2:78027/149‑1 (MQ=255)
         ccGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGc                                                                                                                                          >  3:680244/1‑149 (MQ=255)
            aaCTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCACCAACGATCAGATCGTTGTTGGAga                                                                                                                                       <  2:82166/149‑2 (MQ=255)
             aCTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGtt                                                                                                                                      <  2:692543/149‑1 (MQ=255)
               ttCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCt                                                                                                                                                      >  3:302350/1‑131 (MQ=255)
               ttCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCt                                                                                                                                                      <  4:302350/131‑1 (MQ=255)
               ttCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTgg                                                                                                                                    >  2:93125/1‑149 (MQ=255)
                tCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGgt                                                                                                                                   <  3:517594/149‑1 (MQ=255)
                 cGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGgtg                                                                                                                                  <  3:755090/149‑1 (MQ=255)
                  gCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGgtgt                                                                                                                                 <  2:617583/149‑1 (MQ=255)
                       aGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                               <  3:345876/130‑1 (MQ=255)
                       aGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGtt                                                                                                                                               >  4:345876/1‑130 (MQ=255)
                        gAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCa                                                                                                                           >  1:92292/1‑149 (MQ=255)
                            tGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGttt                                                                                                                       >  1:222/1‑149 (MQ=255)
                               tCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGAt                                                                                                                    <  2:576707/149‑1 (MQ=255)
                                    cgtATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATgtggt                                                                                                               <  4:393478/149‑1 (MQ=255)
                                              tGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAgcgcgc                                                                                                     >  2:285705/1‑149 (MQ=255)
                                                aaaGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCt                                                                                                   >  4:405924/1‑149 (MQ=255)
                                                   gATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGaa                                                                                                <  4:44528/149‑1 (MQ=255)
                                                   gATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGaa                                                                                                <  4:748581/149‑1 (MQ=255)
                                                    aTCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAag                                                                                               <  1:93125/149‑1 (MQ=255)
                                                    aTCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAag                                                                                               >  3:461014/1‑149 (MQ=255)
                                                    aTCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAag                                                                                               >  3:317303/1‑149 (MQ=255)
                                                     tCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAaga                                                                                              <  2:255661/149‑1 (MQ=255)
                                                          tGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGaa                                                                                         <  2:550187/149‑1 (MQ=255)
                                                          tGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGaa                                                                                         <  4:461014/149‑1 (MQ=255)
                                                          tGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGaa                                                                                         <  2:92292/149‑1 (MQ=255)
                                                              gCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACa                                                                                     >  1:412153/1‑149 (MQ=255)
                                                                 aaaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAGTCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGc                                                                                            <  4:349496/139‑1 (MQ=255)
                                                                 aaaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGc                                                                                            >  3:349496/1‑139 (MQ=255)
                                                                  aaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGt                                                                                 <  3:748691/149‑1 (MQ=255)
                                                                  aaaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGt                                                                                 <  4:67763/149‑1 (MQ=255)
                                                                   aaCAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGtt                                                                                <  4:184799/149‑1 (MQ=255)
                                                                      aGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGaa                                                                             <  2:198020/149‑1 (MQ=255)
                                                                      aGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGaa                                                                             <  3:399602/149‑1 (MQ=255)
                                                                      aGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGaa                                                                             <  2:437856/149‑1 (MQ=255)
                                                                         tgctCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCa                                                                          >  4:724644/1‑149 (MQ=255)
                                                                          gctCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAg                                                                         >  1:466353/1‑149 (MQ=255)
                                                                           ctCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGTAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAgt                                                                        >  3:586826/1‑149 (MQ=255)
                                                                            tCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAgtg                                                                       >  3:88731/1‑149 (MQ=255)
                                                                              aaaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGttgttg                                                                                                                                           >  2:735432/1‑79 (MQ=255)
                                                                              aaaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCAGATCGttgttg                                                                                                                                           <  1:735432/79‑1 (MQ=255)
                                                                               aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTAt                                                                    <  4:297425/149‑1 (MQ=255)
                                                                               aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTAt                                                                    <  4:618102/149‑1 (MQ=255)
                                                                               aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTAt                                                                    <  3:405924/149‑1 (MQ=255)
                                                                               aaGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTAt                                                                    <  3:157437/149‑1 (MQ=255)
                                                                                 gggCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGa                                                                  <  4:317303/149‑1 (MQ=255)
                                                                                           aGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAAc                                                        >  4:692115/1‑149 (MQ=255)
                                                                                           aGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAAc                                                        <  1:285705/149‑1 (MQ=255)
                                                                                              tttCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGtt                                                     <  4:555682/149‑1 (MQ=255)
                                                                                              tttCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGtt                                                     >  2:763305/1‑149 (MQ=255)
                                                                                               ttCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTg                                                    <  4:316553/149‑1 (MQ=255)
                                                                                                tCCGAAGAGGGCGCGGTGCCGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGc                                                   >  3:626943/1‑149 (MQ=255)
                                                                                                  cGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCa                                                 <  1:692403/149‑1 (MQ=255)
                                                                                                  cGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCa                                                 <  4:26179/149‑1 (MQ=255)
                                                                                                  cGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCa                                                 >  3:252286/1‑149 (MQ=255)
                                                                                                   gAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCAc                                                >  2:454948/1‑149 (MQ=255)
                                                                                                       agGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGcc                                            >  4:713732/1‑149 (MQ=255)
                                                                                                            gcgGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTg                                       <  4:446345/149‑1 (MQ=255)
                                                                                                              ggTGCAGGTGTTCCGTCCAATCTCCAACAACGGTCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAATCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTggg                                     <  4:579416/149‑1 (MQ=255)
                                                                                                              ggTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTggg                                     <  3:68256/149‑1 (MQ=255)
                                                                                                                  cAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGa                                 >  1:493168/1‑149 (MQ=255)
                                                                                                                   aGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGaa                                <  2:489261/149‑1 (MQ=255)
                                                                                                                     gtgtTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAAtg                              >  2:664077/1‑149 (MQ=255)
                                                                                                                      tgtTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAAtgt                             <  1:454948/149‑1 (MQ=255)
                                                                                                                        ttCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTGGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGc                           >  2:369208/1‑149 (MQ=255)
                                                                                                                         tCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCa                          <  4:252286/149‑1 (MQ=255)
                                                                                                                         tCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCa                          <  4:654693/149‑1 (MQ=255)
                                                                                                                         tCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCa                          <  2:366219/149‑1 (MQ=255)
                                                                                                                          ccGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAg                         >  3:378593/1‑149 (MQ=255)
                                                                                                                            gTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGAc                       <  4:75838/149‑1 (MQ=255)
                                                                                                                              ccAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGAcgc                     <  1:266544/149‑1 (MQ=255)
                                                                                                                               cAATCTCCAACAACGATCTGATCGTTGGTGTAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGAcgcg                    >  3:571870/1‑149 (MQ=255)
                                                                                                                                 aTCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCgaa                  <  3:713732/149‑1 (MQ=255)
                                                                                                                                  tctcCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCgaag                 >  1:725718/1‑149 (MQ=255)
                                                                                                                                       aacaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAAtt            <  4:118512/149‑1 (MQ=255)
                                                                                                                                       aacaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATAAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTc           <  2:365000/149‑1 (MQ=255)
                                                                                                                                        acaacGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTc           <  2:493168/149‑1 (MQ=255)
                                                                                                                                              aTCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTc           <  4:597166/143‑1 (MQ=255)
                                                                                                                                              aTCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTc           >  3:597166/1‑143 (MQ=255)
                                                                                                                                                 tGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGtt                                                     >  2:216063/1‑98 (MQ=255)
                                                                                                                                                 tGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGtt                                                     <  1:216063/98‑1 (MQ=255)
                                                                                                                                                 tGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGGCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAGTTc  <  2:466353/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACCGGTATTCCGAACTTCGCACCAGAACTGTTCCGTCGTATCCGCCTGAAAGATCCGCTGAAGCAAAAACAGCTGCTCAAAGGGCTGGTACAGCTTTCCGAAGAGGGCGCGGTGCAGGTGTTCCGTCCAATCTCCAACAACGATCTGATCGTTGGTGCAGTTGGTGTGCTGCAGTTTGATGTGGTGGTAGCGCGCCTGAAGAGCGAATACAACGTTGAAGCAGTGTATGAGTCAGTCAACGTTGCCACTGCCCGCTGGGTAGAATGTGCAGACGCGAAGAAATTCGAAGAGTTC  >  NC_000913/4610575‑4610868

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