Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,272,016 T→A 6.0% F431I (TTT→ATT)  rtn → resistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,272,0160TA 6.0% 246.1 / 4.9 92F431I (TTT→ATT) rtnresistance protein for phages lambda and N4, putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (36/48);  new base A (3/3);  total (41/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.33e-01

TCCAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAA  >  NC_000913/2271870‑2272162
                                                                                                                                                  |                                                                                                                                                  
tCCAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGtttt                                                                                                                                                  <  2:1363925/149‑1 (MQ=255)
 ccAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGATAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCCGTCGGCGAAGAAATTGATATTGATGACTTCTGCACCGGGCACAGCGCGCTTATCTATATTGAGCCTATTa                                                                                                                                                 >  2:1243384/1‑149 (MQ=255)
  cAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGACGCCACACAACTCGTCGCCGTGCTGCACTCGGTCGGCGGAGAAATTGCTATTGATGACTTCGGCACCGGGCCCAGCGCGCTTATCGATCTTGAGCGTGTTAc                                                                                                                                                >  1:321684/1‑149 (MQ=255)
         tgCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTAAAACGCTcaag                                                                                                                                         <  1:829570/149‑4 (MQ=255)
           cTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGAATTAGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTa                                                                                                                                       >  1:464586/1‑149 (MQ=255)
            tGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAATCCCCCCAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTTTCTTTCTTGAGCGGTTTTCGCTCGAGTAt                                                                                                                                      >  2:876081/1‑149 (MQ=255)
            tGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTAAAACGCTcaagata                                                                                                                                      <  1:713958/149‑7 (MQ=255)
             ggAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATc                                                                                                                                     <  1:1108790/149‑1 (MQ=255)
             ggAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGGAGAAATTGCGATTGCTGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGGTTTTCGCTCGATTAta                                                                                                                                     >  2:1095881/1‑148 (MQ=255)
                 aatACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGaa                                                                                                                                 <  2:1500971/147‑1 (MQ=255)
                 aTTACCGAGCGCGATATGTTGATAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGTGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGaa                                                                                                                                 <  1:1194631/149‑1 (MQ=255)
                 aTTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGaa                                                                                                                                 <  1:1513369/149‑1 (MQ=255)
                   tACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGGTTTACGCTCGATTATCTGaaaa                                                                                                                               <  2:1255106/149‑1 (MQ=255)
                     ccGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAAtt                                                                                                                             <  1:1163618/149‑1 (MQ=255)
                      cGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTg                                                                                                                            >  2:250064/1‑149 (MQ=255)
                      cGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTg                                                                                                                            >  2:10044/1‑149 (MQ=255)
                             gATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTgg                                                                                                                     >  1:56042/1‑149 (MQ=255)
                               tatGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGTCGTAGAAATTGCTATTGATGACTTCGTCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTGGATTATCTGAAAATTGACCGTGGAt                                                                                                                   <  1:1043770/149‑1 (MQ=255)
                               tatGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTAAAACGCTcaa                                                                                                                                          <  2:451008/126‑3 (MQ=255)
                               tatGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTAAAACGCTcaa                                                                                                                                          >  1:451008/1‑124 (MQ=255)
                               tatGTTGAAAGAGCAAGAAGCCACACAACACTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGAt                                                                                                                   <  1:942659/149‑1 (MQ=255)
                                 tGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGAttt                                                                                                                 <  1:1142706/149‑1 (MQ=255)
                                  gTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGCCCGTGGATTTa                                                                                                                >  2:670559/1‑149 (MQ=255)
                                  gTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTa                                                                                                                >  2:940556/1‑149 (MQ=255)
                                   ttGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGCCCGTGGATTTag                                                                                                               >  2:1012755/1‑148 (MQ=255)
                                    tGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATc                                                                                                              <  2:125062/149‑1 (MQ=255)
                                      aaaGCGCAAGAAGCCACACAACTCTTCGCCTTGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCTCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCaa                                                                                                            <  2:1023596/149‑1 (MQ=255)
                                       aaGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAAc                                                                                                           <  2:540060/149‑1 (MQ=255)
                                           gCAAGAAGCCACACAACTCTTCGCCTGGCAGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCa                                                                                                       >  2:709184/1‑149 (MQ=255)
                                             aagaagACACACAACTCTTCGACTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGCCACCGGCCAAAGCGCGCTTATATATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGAAAGTGGATTTATCAACGCCATc                                                                                                     >  1:122191/1‑149 (MQ=255)
                                              agaagCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCg                                                                                                    <  1:1530462/149‑1 (MQ=255)
                                               gaagCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCgg                                                                                                   <  2:993478/149‑1 (MQ=255)
                                                  gCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTAc                                                                                                <  1:453733/149‑1 (MQ=255)
                                                   ccACACAACTCTTCGCCTGGCTGCACTCGGTCGGAGTAGAAATTGCTCTTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATATTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGCGGATTTATCAACGCCATCGGTACg                                                                                               <  2:1428949/149‑1 (MQ=255)
                                                    cacacaACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGCCTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACgg                                                                                              >  1:1174678/1‑149 (MQ=255)
                                                     acacaACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCCCCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGa                                                                                             >  1:463149/1‑149 (MQ=255)
                                                     acacaACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGa                                                                                             >  2:1498840/1‑149 (MQ=255)
                                                     acacaACTCCTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACACTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGa                                                                                             >  2:1287500/1‑149 (MQ=255)
                                                      cacaACTCTTCGCCTGGCTGCACTAGGTCGGCGTACAAATTGCTATTGATGACTTCGGCACCGGGCACACCGCGCTTATCTAACTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCCTGGATTTATCAACGCCATCGCTAAGGaa                                                                                            >  2:616880/1‑149 (MQ=255)
                                                        caACTCTTCGCCGGGCTGCACTCGGTCGGCGTAGAAATTGCTCTTTATGACTTCGGCACGGGGAACATCGCTCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGAGTTATCAACGCCATCGTTACGGAAAc                                                                                          <  2:952481/149‑1 (MQ=255)
                                                         aaCTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACg                                                                                         >  1:557843/1‑149 (MQ=255)
                                                           ctctTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGAt                                                                                       <  1:188237/149‑1 (MQ=255)
                                                           ctctTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGAt                                                                                       <  2:158912/149‑1 (MQ=255)
                                                            tctTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATc                                                                                      <  2:799857/149‑1 (MQ=255)
                                                             gttCGCCTGGCTGGACTCTGTCGGCGTAGAAATTGCTATTTATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCa                                                                                     <  2:96847/148‑1 (MQ=255)
                                                                cGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACtt                                                                                  <  2:571811/149‑1 (MQ=255)
                                                                   cTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc                                                                               <  2:247346/149‑1 (MQ=255)
                                                                   cTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc                                                                               <  1:1312064/149‑1 (MQ=255)
                                                                     ggCTGCACTCGCTCCGCGTAGAAATTGCTATTGATGACTTCGGCACCTGGGATTGCGCGCTTATATATCTTGATCGTTTTACGCTCGATTATCTGAATTTTGACCGTGGATTTATCAACGCCATCGGTACTGAAACTATCTATTccccc                                                                             <  2:122191/149‑1 (MQ=255)
                                                                        tGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCAGAAAATTGACCTTGGATTTATCAACACCATCGGTACGGATACGATCACTGCCCCCGTa                                                                          >  2:839531/1‑149 (MQ=255)
                                                                        tGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTAATCTATCTTGAGCGTATTACGCACGATAATCTGAAAATTGCCCGTCGATTTATCAACGCCATCGGTACGGAACCGATCACTTCCCCCGTa                                                                          >  2:260088/1‑149 (MQ=255)
                                                                          cACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACt                                                                        <  2:256476/149‑1 (MQ=255)
                                                                            cTCGGTCGGCGTAGAAATTGCTATTTATGACTTCGGCACCGGGCACAGCGCGCTTATATATCTTGAGCGTTTTACGCTCGATTATCTGAATATTGTACGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTg                                                                      <  2:464586/149‑1 (MQ=255)
                                                                            cTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTg                                                                      >  1:740628/1‑149 (MQ=255)
                                                                                gtcggCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGAcgc                                                                  >  2:946377/1‑149 (MQ=255)
                                                                                    gCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTg                                                              <  1:334374/149‑1 (MQ=255)
                                                                                      gTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACACTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCt                                                            <  1:1287500/149‑1 (MQ=255)
                                                                                          aaaTTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTTAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCTACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGATGACg                                                        >  2:1248474/1‑149 (MQ=255)
                                                                                           aaTTGCGATTGATGACTTGGGCACCGGGGACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGc                                                       <  1:1191995/149‑1 (MQ=255)
                                                                                              tGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCAACGGTACGGAAACGATCGCTTCCCCCGTACTTGACGCGGTGCTGACGCTgg                                                    >  1:1547725/1‑149 (MQ=255)
                                                                                               cctctagatgaCTACCGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTATACGCTCGATTATCTGAAAATTGACCGTCGATTTATCAACGCCATCGGTACTGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGc                                                   <  2:1500720/143‑1 (MQ=255)
                                                                                                 tATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTGTTACGCACGATTATCTGAAAATTGACCGTGGATTTATCAACCCCATCGGTACGGAAACGATCACATCACCCGTACTGGACGCGGTGCTGAAGCTGGCGa                                                 >  2:1003367/1‑149 (MQ=255)
                                                                                                  aTTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGaa                                                >  2:410195/1‑149 (MQ=255)
                                                                                                   ttGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGaaa                                               >  2:1312027/1‑149 (MQ=255)
                                                                                                   ttGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGaaa                                               <  2:477196/149‑1 (MQ=255)
                                                                                                    tgatgaCTTCGGCACCGGGCACAGCGCGGTTATCTATCTTGAGCGTTTTAATCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAAc                                              <  1:634273/149‑1 (MQ=255)
                                                                                                    tgatgaCTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAAc                                              <  1:250064/149‑1 (MQ=255)
                                                                                                         aCTTCGGCACCGGGCACAGCGCGCTTAGCTATCTTGAGCGTTTTACGCTCGTTAATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTc                                         <  2:1172346/149‑1 (MQ=255)
                                                                                                           ttCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCaa                                       >  2:1434943/1‑149 (MQ=255)
                                                                                                           ttCGGCACCGGGCACAGCGCGCTTATCTATCTCGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCaa                                       >  2:492613/1‑149 (MQ=255)
                                                                                                            tCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAat                                      <  2:942109/149‑1 (MQ=255)
                                                                                                               gCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATg                                   <  2:647057/149‑1 (MQ=255)
                                                                                                                cACCGGGCACAGCGCTCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTTGCGAAACGCCTCAATATGc                                  <  2:1361218/149‑1 (MQ=255)
                                                                                                                cACCGGGCACAGCGCGCTTAGCTATCTTGAGCGTTTTACTCTCGATTATCTGAACATTGACCGTGGTTTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGc                                  <  1:1462767/149‑1 (MQ=255)
                                                                                                                    gggCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAACCGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGAc                              >  2:232224/1‑149 (MQ=255)
                                                                                                                    gggCACAGCGCGATTATCTTTCTTGAGCGTTTTACGCTCGATTATCTGAAAAATGACCGTGGATTTATCAACGCCATCGGTACGGACACGATCACTTCCCCCGTACTTGCCGCGGTGCTGACGCGGGCGAAACGCCTCAATATGCTCAc                              <  1:735578/149‑1 (MQ=255)
                                                                                                                       cacaTCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACCGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGAGAAACGCCTCAATATGCTGACGGt                           <  1:1243384/149‑1 (MQ=255)
                                                                                                                            cgcgcTTATCTATCTTGAGCGATTTACGCTCGATTATCTGAAAATTCACCGTGGAATTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTTCTGACGCTGGCGAAACGCCTCAATATGATGCCGGTAGCtc                      >  2:417578/1‑148 (MQ=255)
                                                                                                                              cgcTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGaa                    <  2:988928/149‑1 (MQ=255)
                                                                                                                               gcTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAg                   <  1:1263396/149‑1 (MQ=255)
                                                                                                                               gcTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCGGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAg                   <  1:670559/149‑1 (MQ=255)
                                                                                                                                 ttATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGCCGCTGGCGAACCGCCTCAATATGCTGACGGTTGCTGAAggg                 >  2:1387147/1‑149 (MQ=255)
                                                                                                                                  tagctatctTGTGCGTTTTACGCTAGATTATCAGAAAATTGACCGTGGATTTCTCAACGCCATCGGTACGGAAACGATCACAGCCCCCGTACTTGACCCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCGGAAgggg                <  2:706479/146‑1 (MQ=255)
                                                                                                                                   atctatctTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGt               <  1:10044/149‑1 (MQ=255)
                                                                                                                                   atctatctTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCGGAAGGGGt               >  1:366977/1‑149 (MQ=255)
                                                                                                                                     ctatctTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCg             <  2:557843/149‑1 (MQ=255)
                                                                                                                                     ctatctTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACATCCCCCGAACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCg             >  1:360023/1‑149 (MQ=255)
                                                                                                                                       atctTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATAAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTCACGCTGGCGAAACGCCTCGATATGCTGACGGTTGCTGAAGGGGTCGaa           >  2:969302/1‑149 (MQ=255)
                                                                                                                                          ttGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACg        >  2:821896/1‑149 (MQ=255)
                                                                                                                                           tGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGc       >  1:1193297/1‑149 (MQ=255)
                                                                                                                                            gAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGcc      <  2:56042/149‑1 (MQ=255)
                                                                                                                                               cGTTTTACGCGCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGa   >  2:90190/1‑149 (MQ=255)
                                                                                                                                                gTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGCAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGaa  >  2:210180/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCCAGATTGTGCTGGAAATTACCGAGCGCGATATGTTGAAAGAGCAAGAAGCCACACAACTCTTCGCCTGGCTGCACTCGGTCGGCGTAGAAATTGCTATTGATGACTTCGGCACCGGGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCCCCGTACTTGACGCGGTGCTGACGCTGGCGAAACGCCTCAATATGCTGACGGTTGCTGAAGGGGTCGAAACGCCGGAA  >  NC_000913/2271870‑2272162

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: