Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,251,140 G→C L239V (CTT→GTT)  yhaJ ← putative transcriptional regulator LYSR‑type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,251,1400GC94.4% 269.8 / 7.3 90L239V (CTT→GTT) yhaJputative transcriptional regulator LYSR‑type
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base C (56/28);  total (61/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

TCTCGCTGGTCGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGG  >  NC_000913/3251050‑3251228
                                                                                          |                                                                                        
tctcGCTGGTCGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAAccc                                                                                        >  2:4194348/1‑93 (MQ=255)
tctcGCTGGTCGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAAccc                                                                                        >  2:169564/1‑93 (MQ=255)
tctcGCTGGGCGATTCCGGGCGGACGACACGCAACCGCCCTTCCGCAATAGCTTTTTCGACCATAGGATACGGCATCGTCGCCACTCCAAccc                                                                                        >  2:4194347/1‑93 (MQ=255)
  tcGCTGGTCGATTCCGGGCTGACGACAAGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACcccc                                                                                      >  2:1110857/1‑93 (MQ=255)
      tGGTCGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCa                                                                                  >  1:1476105/1‑93 (MQ=255)
      tGGTCGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCa                                                                                  >  2:1327282/1‑93 (MQ=255)
         gcGATTCCGGGCTGACGACCCGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTa                                                                               <  1:1631324/92‑1 (MQ=255)
          cGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCTACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTaa                                                                              <  2:2650526/93‑1 (MQ=255)
            aTTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATg                                                                            <  2:1591969/93‑1 (MQ=255)
              tCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGcc                                                                          >  1:1366009/1‑93 (MQ=255)
               ccGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCt                                                                         >  1:3804536/1‑93 (MQ=255)
               ccGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCt                                                                         >  2:2572423/1‑93 (MQ=255)
               ccGGGCTGACGACACGCAACCGCCCTCCCGCAATATCTTTTCCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGctt                                                                         <  1:2327894/93‑3 (MQ=255)
                cGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTg                                                                        >  2:1070856/1‑93 (MQ=255)
                  ggCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGAc                                                                      >  1:466389/1‑93 (MQ=255)
                  ggCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGAc                                                                      <  1:3575892/93‑1 (MQ=255)
                   gCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACg                                                                     >  2:960115/1‑93 (MQ=255)
                      gacgacACGCAACCGCCCTTCCGAAATATCTTTTTCGACCATAGGATAGGGCAGCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGttt                                                                  <  2:2950336/93‑1 (MQ=255)
                       acgacACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTa                                                                 <  2:1639116/93‑1 (MQ=255)
                         gacACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTAtt                                                               >  2:3795881/1‑92 (MQ=255)
                         gacACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATc                                                               >  2:2520440/1‑93 (MQ=255)
                               cAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAAt                                                         <  2:2051218/93‑1 (MQ=255)
                                aaCCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATc                                                        >  2:451463/1‑93 (MQ=255)
                                aaCCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATc                                                        >  1:3759179/1‑93 (MQ=255)
                                 aCCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATct                                                       >  2:2837798/1‑92 (MQ=255)
                                 aCCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCg                                                       <  1:4384205/93‑1 (MQ=255)
                                   cGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTg                                                     >  2:2890735/1‑93 (MQ=255)
                                   cGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTg                                                     >  2:2890736/1‑93 (MQ=255)
                                    gCCCTTCCGCAATATCTTTTTCGATCTTATGATACGGCATCGTCGCCACTCCAACCCCCGCCAGTAATGCCTGACGTTTACCTTCAATTGTGc                                                    >  2:4349764/1‑93 (MQ=255)
                                    gCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGc                                                    >  1:4664380/1‑93 (MQ=255)
                                    gCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGc                                                    <  1:2022408/93‑1 (MQ=255)
                                    gCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGc                                                    >  1:78059/1‑93 (MQ=255)
                                       cTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCa                                                 >  2:290542/1‑93 (MQ=255)
                                        ttCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCAc                                                >  2:1791342/1‑93 (MQ=255)
                                            gCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGtt                                            >  2:2781603/1‑93 (MQ=255)
                                            gCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGtt                                            >  2:4469635/1‑93 (MQ=255)
                                             cAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTa                                           >  1:1450891/1‑93 (MQ=255)
                                              aaTATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTaa                                          <  2:4185111/93‑1 (MQ=255)
                                              aaTATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTaa                                          >  2:764824/1‑93 (MQ=255)
                                              aaTATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTaa                                          <  2:2585224/93‑1 (MQ=255)
                                               atatCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAg                                         >  1:4466802/1‑93 (MQ=255)
                                               atatCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAg                                         >  2:269806/1‑93 (MQ=255)
                                               atatCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGCATATCTTCCATCGTGCTCACCGTTaac                                         >  1:4466801/1‑92 (MQ=255)
                                                tatTTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAgc                                        <  2:2626062/93‑1 (MQ=255)
                                                tatCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAgc                                        >  2:2029523/1‑93 (MQ=255)
                                                  tCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAgcgc                                      >  1:3966811/1‑93 (MQ=255)
                                                  tCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAgcgc                                      >  1:1879983/1‑93 (MQ=255)
                                                  tCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAgcgc                                      <  1:1646942/93‑1 (MQ=255)
                                                  tCTTTTGCGACCATAGGATACGGCATCGTCGCCCCGCCCCCCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAgcgc                                      <  1:1646941/93‑1 (MQ=255)
                                                    tttttCGACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCgg                                    >  1:3810386/1‑93 (MQ=255)
                                                    tttttCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCgg                                    >  2:1139322/1‑93 (MQ=255)
                                                     ttttCGCCCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGc                                   <  1:3351609/93‑1 (MQ=255)
                                                       ttCGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTg                                 >  1:3274297/1‑93 (MQ=255)
                                                         cGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGtt                               <  1:3911777/93‑1 (MQ=255)
                                                         cGACCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGtt                               <  1:322549/93‑1 (MQ=255)
                                                           aCCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTCAGCGCGGCTGTTTg                             >  1:2724979/1‑93 (MQ=255)
                                                           aCCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTg                             >  1:2724978/1‑93 (MQ=255)
                                                           aCCATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTg                             >  1:3400273/1‑93 (MQ=255)
                                                             cATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTc                           >  1:3261640/1‑93 (MQ=255)
                                                             cATAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTc                           >  1:944571/1‑93 (MQ=255)
                                                              aTAGGATACGGCATCGTCGCCACGCCAATCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGATTGTcc                          >  1:845214/1‑93 (MQ=255)
                                                              aTAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTcc                          >  2:2126732/1‑93 (MQ=255)
                                                              aTAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTcc                          >  1:3673831/1‑93 (MQ=255)
                                                               tAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCca                         >  2:880846/1‑93 (MQ=255)
                                                               tAGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCca                         >  1:3780145/1‑93 (MQ=255)
                                                                aGGATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCcag                        <  1:83328/93‑1 (MQ=255)
                                                                 ggATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCcagc                       >  2:376853/1‑93 (MQ=255)
                                                                  gATACGGCATCGTCTCCACGCTAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCcagca                      <  2:3267561/93‑1 (MQ=255)
                                                                  gATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCcagca                      >  2:1988738/1‑93 (MQ=255)
                                                                  gATACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCcagca                      <  1:1198379/93‑1 (MQ=255)
                                                                    tACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCagcagc                    >  2:4378613/1‑93 (MQ=255)
                                                                    tACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCagcagc                    >  2:2181905/1‑93 (MQ=255)
                                                                    tACGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCagcagc                    <  1:4003158/93‑1 (MQ=255)
                                                                      cGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTg                  <  1:2929170/93‑1 (MQ=255)
                                                                      cGGCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTg                  >  2:162602/1‑93 (MQ=255)
                                                                       ggCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGt                 <  2:1567557/93‑1 (MQ=255)
                                                                       ggCATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGt                 >  2:4452452/1‑93 (MQ=255)
                                                                         cATCGTCGCCACGCCACCCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTAc               <  2:168728/93‑1 (MQ=255)
                                                                         cATCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTAc               <  1:916848/93‑1 (MQ=255)
                                                                         cATCGTCCCCACGCCACCCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTAc               <  1:752007/93‑1 (MQ=255)
                                                                          aTCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACg              >  1:2884451/1‑93 (MQ=255)
                                                                          aTCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACg              >  2:1835380/1‑93 (MQ=255)
                                                                          aTCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACg              >  2:3864803/1‑93 (MQ=255)
                                                                          aTCGTCGCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACg              >  2:3113183/1‑93 (MQ=255)
                                                                           tcgtcgCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACgg             >  1:651319/1‑93 (MQ=255)
                                                                           tcgtcgCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACgg             >  1:1964128/1‑93 (MQ=255)
                                                                           tcctcgCCACGCCAACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACgg             <  2:1409395/90‑1 (MQ=255)
                                                                               cgccacgccaCCCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCaa         <  1:4570244/93‑1 (MQ=255)
                                                                                gccacgccaACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAac        >  1:1072239/1‑93 (MQ=255)
                                                                                 ccacgccaACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAaca       >  2:1372696/1‑93 (MQ=255)
                                                                                  cacgccaACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAacac      >  1:1786669/1‑93 (MQ=255)
                                                                                    tgccaACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCg    <  2:2890722/92‑1 (MQ=255)
                                                                                     gccaAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCgg   >  2:4039049/1‑93 (MQ=255)
                                                                                     gccaACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCgg   >  2:4131137/1‑93 (MQ=255)
                                                                                      ccaACCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCggg  >  2:929654/1‑93 (MQ=255)
                                                                                          |                                                                                        
TCTCGCTGGTCGATTCCGGGCTGACGACACGCAACCGCCCTTCCGCAATATCTTTTTCGACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGG  >  NC_000913/3251050‑3251228

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑