Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,527,897 A→G T318A (ACC→GCC)  yfeH → putative inorganic ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,527,8970AG98.6% 254.4 / ‑3.8 71T318A (ACC→GCC) yfeHputative inorganic ion transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (34/36);  total (34/37)
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.62e-01

AAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCAT  >  NC_000913/2527753‑2528044
                                                                                                                                                |                                                                                                                                                   
aaaaaGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAAc                                                                                                                                               <  2:503385/151‑1 (MQ=255)
aaaaaGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAAc                                                                                                                                               <  2:65706/151‑1 (MQ=255)
aaaaaGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAAc                                                                                                                                               >  2:724313/1‑151 (MQ=255)
  aaaGAGTCTGGCAAATGGCATCCCGATGGCCCACATTCTGTTCCCCACATCGGTGATCGCTATGATGGTGCTGCCCCTGATGATTTTCCATTAGATCCAATTGATGGTCTGTGCGATGCTGGCCCGTCGATACCAACGCCCGGCCGCACAg                                                                                                                                             <  1:116607/151‑1 (MQ=255)
     gagTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTa                                                                                                                                          <  2:852094/151‑1 (MQ=255)
         cTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAgg                                                                                                                                      >  2:174394/1‑151 (MQ=255)
         cTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAgg                                                                                                                                      >  1:1146576/1‑151 (MQ=255)
          tGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGgc                                                                                                                                     >  2:1223840/1‑151 (MQ=255)
           ggCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAgg                                                                                                                                      <  2:86605/149‑1 (MQ=255)
           ggCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAgg                                                                                                                                      <  1:532312/149‑1 (MQ=255)
             cAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCgca                                                                                                                                  <  2:835648/151‑1 (MQ=255)
              aaaTGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCgcag                                                                                                                                 <  2:719343/151‑1 (MQ=255)
               aaTGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCgcagc                                                                                                                                >  1:764749/1‑151 (MQ=255)
                    cATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGaa                                                                                                                           <  1:1130641/151‑1 (MQ=255)
                      tCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGaa                                                                                                                           >  1:1263216/1‑149 (MQ=255)
                        ccGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAag                                                                                                                         >  1:506575/1‑149 (MQ=255)
                         cGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAagcag                                                                                                                      <  2:1172578/151‑1 (MQ=255)
                         cGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAagc                                                                                                                        >  2:1073836/1‑149 (MQ=255)
                                   cATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATaa                                                                                                             >  2:561843/1‑150 (MQ=255)
                                    aTTCTGTTCCACACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATa                                                                                                              >  1:1007043/1‑148 (MQ=255)
                                           tCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTaaa                                                                                                     <  2:580611/150‑1 (MQ=255)
                                            ccccACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAg                                                                                                    <  2:288479/150‑1 (MQ=255)
                                              ccACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCgg                                                                                                 >  1:562482/1‑151 (MQ=255)
                                                acaTCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGAc                                                                                               <  1:7650/151‑1 (MQ=255)
                                                  aTCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGc                                                                                             >  1:114352/1‑151 (MQ=255)
                                                      gTGAGCGGTATGATGGTGATGCCCCTGGTGATTTTCCATCAGATCCACTTGATGGTCTGTGCGGTGCGGGCGCGACGACACAAAAGCCAGGCCGGAGAGATTCAGGAGAAGCAGGAAAGCAGCGCCGATAAAACTTAAAGCGGTCg                                                                                              >  1:1066049/1‑146 (MQ=255)
                                                        gATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAgg                                                                                       <  2:859732/151‑1 (MQ=255)
                                                         aTCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAggg                                                                                      >  1:501471/1‑151 (MQ=255)
                                                            ggTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCt                                                                                   >  2:571444/1‑151 (MQ=255)
                                                               atgatgGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGa                                                                                <  1:152814/151‑1 (MQ=255)
                                                                   tgGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGAcc                                                                              >  1:712579/1‑149 (MQ=255)
                                                                    ggTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGt                                                                           >  1:326738/1‑151 (MQ=255)
                                                                    ggTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAg                                                                            <  2:1027970/150‑1 (MQ=255)
                                                                     gTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGtt                                                                          >  2:1032548/1‑151 (MQ=255)
                                                                       gctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGc                                                                        <  1:1184851/151‑1 (MQ=255)
                                                                         tgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGt                                                                      >  2:785779/1‑151 (MQ=255)
                                                                         tgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGt                                                                      >  2:1231309/1‑151 (MQ=255)
                                                                          gcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTc                                                                     >  2:429974/1‑151 (MQ=255)
                                                                            cccTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTc                                                                     >  2:69284/1‑149 (MQ=255)
                                                                                 atgatTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCt                                                              <  1:803305/151‑1 (MQ=255)
                                                                                 atgatTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCt                                                              <  2:541787/151‑1 (MQ=255)
                                                                                  tgatTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTc                                                             >  1:543158/1‑151 (MQ=255)
                                                                                   gatTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCg                                                            <  1:269084/151‑1 (MQ=255)
                                                                                    atTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCgg                                                           <  2:944681/151‑1 (MQ=255)
                                                                                        tCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGttt                                                        <  2:818663/150‑1 (MQ=255)
                                                                                         ccATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGtttt                                                       <  2:424945/150‑1 (MQ=255)
                                                                                            tCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATc                                                   >  2:302695/1‑151 (MQ=255)
                                                                                             cAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATca                                                  <  1:724279/151‑1 (MQ=255)
                                                                                                aTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATcagca                                               <  2:788443/151‑1 (MQ=255)
                                                                                                   cAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCg                                            >  2:987635/1‑151 (MQ=255)
                                                                                                   cAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCg                                            <  1:89362/151‑1 (MQ=255)
                                                                                                     aTTGATGGTCTGTGCGGTGCTGGCGCGTCGCTACAAACGCCAGGCCGAACAGTTACATGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGt                                           <  1:309214/150‑1 (MQ=255)
                                                                                                       tGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCTCCTCGGTTTTAATCAGCAGCGTGa                                         <  1:526056/150‑1 (MQ=255)
                                                                                                       tGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGAGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGAt                                        >  2:1141073/1‑151 (MQ=255)
                                                                                                         aTGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGAttt                                      >  1:298360/1‑151 (MQ=255)
                                                                                                         aTGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGAttt                                      >  2:214997/1‑151 (MQ=255)
                                                                                                             tCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAt                                  >  1:27015/1‑151 (MQ=255)
                                                                                                                gtgCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAtt                                 >  2:558317/1‑149 (MQ=255)
                                                                                                                gtgCGGTGCTGGCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAtt                                 >  2:326069/1‑149 (MQ=255)
                                                                                                                          ggCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                       <  2:983064/149‑1 (MQ=255)
                                                                                                                          ggCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                       <  2:61928/149‑1 (MQ=255)
                                                                                                                          ggCGCGTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGa                       <  2:630943/149‑1 (MQ=255)
                                                                                                                              cgTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTcc                 >  1:264629/1‑151 (MQ=255)
                                                                                                                              cgTCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTcc                 <  2:564010/151‑1 (MQ=255)
                                                                                                                                tCGATACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCt               >  2:693713/1‑151 (MQ=255)
                                                                                                                                   aTACAAACGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGcc            <  2:1122554/151‑1 (MQ=255)
                                                                                                                                       aaaCGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGa        <  1:221200/151‑1 (MQ=255)
                                                                                                                                        aaCGCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGaa       <  2:246565/151‑1 (MQ=255)
                                                                                                                                           gCCAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACt     <  1:1123827/150‑1 (MQ=255)
                                                                                                                                            ccAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCa   <  2:441276/151‑1 (MQ=255)
                                                                                                                                             cAGGCCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCAt  >  2:171661/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AAAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCACCGAGCTTCCCTGCCGGGAACTCAT  >  NC_000913/2527753‑2528044

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