Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,781 A→G V147A (GTT→GCT)  yhiJ ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,7810AG87.8% 886.8 / 101.6 344V147A (GTT→GCT) yhiJhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/23);  new base G (126/176);  total (145/199)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 349 total aligned reads displayed.

GGGTTTTATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCAACCG  >  NC_000913/3629690‑3629870
                                                                                           |                                                                                         
 ggTTTTATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTgcagca                                                                                         >  1:512183/1‑93 (MQ=255)
   ttttATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGcatc                                                                                       >  1:4275039/1‑93 (MQ=255)
   ttttATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAAcatc                                                                                       >  2:711925/1‑93 (MQ=255)
     ttATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAAcatcat                                                                                     >  1:5806958/1‑93 (MQ=255)
      tATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCAtt                                                                                    <  1:2964531/93‑1 (MQ=255)
          aTTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCa                                                                                <  2:573117/93‑1 (MQ=255)
          aTTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCa                                                                                <  1:5984355/93‑1 (MQ=255)
           ttGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCAt                                                                               <  1:6394115/93‑1 (MQ=255)
             gTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTa                                                                             <  2:4725554/93‑1 (MQ=255)
             gTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTa                                                                             <  2:3315059/93‑1 (MQ=255)
             gTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTa                                                                             <  2:7310162/93‑1 (MQ=255)
             gTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTa                                                                             <  1:1766023/93‑1 (MQ=255)
                 cGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCAt                                                                         <  1:4308226/93‑1 (MQ=255)
                  gAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCAtt                                                                        <  1:7434786/93‑1 (MQ=255)
                    aTTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTa                                                                      <  1:5467560/93‑1 (MQ=255)
                    aTTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTa                                                                      >  2:3776304/1‑93 (MQ=255)
                    aTTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTa                                                                      >  2:233081/1‑93 (MQ=255)
                      ttAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACt                                                                    >  2:5035457/1‑93 (MQ=255)
                      ttAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACt                                                                    >  2:3134132/1‑93 (MQ=255)
                          tatCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCa                                                                <  1:5566942/93‑1 (MQ=255)
                            tCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCata                                                              >  1:1139924/1‑93 (MQ=255)
                            tCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCata                                                              <  1:3773297/93‑1 (MQ=255)
                            tCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCata                                                              <  2:1042406/93‑1 (MQ=255)
                            tCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCata                                                              <  1:4992939/93‑1 (MQ=255)
                            tCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCata                                                              >  1:84324/1‑93 (MQ=255)
                             cAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCatat                                                             <  2:6773301/93‑1 (MQ=255)
                              aGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCAtatt                                                            >  2:2225816/1‑93 (MQ=255)
                               gATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCAtatta                                                           <  1:1284715/93‑1 (MQ=255)
                                aTACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCAtattat                                                          >  1:5263292/1‑93 (MQ=255)
                                aTACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCAtattat                                                          <  2:5657121/93‑1 (MQ=255)
                                 tACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATattatt                                                         >  2:6277034/1‑93 (MQ=255)
                                 tACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATattatt                                                         <  1:6563274/93‑1 (MQ=255)
                                  aCATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATAttatta                                                        <  2:5305009/93‑1 (MQ=255)
                                    aTTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAg                                                      <  1:2537977/93‑1 (MQ=255)
                                     ttATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGt                                                     >  2:2092861/1‑93 (MQ=255)
                                      tataTTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGta                                                    <  1:1095046/93‑1 (MQ=255)
                                      tataTTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGta                                                    >  1:3186011/1‑93 (MQ=255)
                                       atatTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGtaa                                                   >  1:3894286/1‑93 (MQ=255)
                                        tatTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGtaaa                                                  >  1:7001775/1‑93 (MQ=255)
                                        tatTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGtaaa                                                  >  2:635215/1‑93 (MQ=255)
                                         atTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGtaaat                                                 >  2:1347783/1‑93 (MQ=255)
                                          ttGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGtaaata                                                <  1:3653503/93‑1 (MQ=255)
                                            gCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGtaaataaa                                              <  2:2618181/93‑1 (MQ=255)
                                            gCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGtaaataaa                                              <  2:4903555/93‑1 (MQ=255)
                                              gATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACa                                            >  2:685059/1‑93 (MQ=255)
                                              gATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACa                                            >  1:187886/1‑93 (MQ=255)
                                              gATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACa                                            >  1:4600619/1‑93 (MQ=255)
                                              gATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACa                                            >  1:2479877/1‑93 (MQ=255)
                                               atatTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACaa                                           >  1:6242840/1‑93 (MQ=255)
                                               atatTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACaa                                           >  1:5758762/1‑93 (MQ=255)
                                                tatTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACann                                          >  1:1643405/1‑91 (MQ=255)
                                                tatTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAAt                                          >  1:690357/1‑93 (MQ=255)
                                                tatTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAAt                                          <  1:1012034/93‑1 (MQ=255)
                                                 atTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATa                                         <  1:2043598/93‑1 (MQ=255)
                                                  ttttCCAGTTCCAGTTGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATAc                                        <  1:6280904/93‑1 (MQ=255)
                                                  ttttCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATAc                                        <  2:6928554/93‑1 (MQ=255)
                                                   tttCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACt                                       <  2:3914347/93‑1 (MQ=255)
                                                      ccagttccagATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCa                                    <  2:1460486/93‑1 (MQ=255)
                                                       nngttccagATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCAc                                   <  1:1330768/91‑1 (MQ=255)
                                                       cagttccagATGGGTTTGTGCGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCAc                                   >  2:642672/1‑93 (MQ=255)
                                                          ttccagATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCaa                                >  2:4075023/1‑93 (MQ=255)
                                                            ccagATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAAtt                              <  1:6820773/93‑1 (MQ=255)
                                                            ccagATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAAtt                              >  2:7033742/1‑93 (MQ=255)
                                                             cagATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTa                             >  1:138963/1‑93 (MQ=255)
                                                               gATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAat                           >  1:3981224/1‑93 (MQ=255)
                                                               gATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAat                           >  2:5188469/1‑93 (MQ=255)
                                                               gATGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAat                           >  1:3247997/1‑93 (MQ=255)
                                                                aTGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAata                          <  2:5894564/93‑1 (MQ=255)
                                                                aTGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAata                          <  1:6089219/93‑1 (MQ=255)
                                                                aTGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAata                          <  1:2202557/93‑1 (MQ=255)
                                                                aTGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAata                          <  2:4132785/93‑1 (MQ=255)
                                                                 tGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAatat                         >  1:6948012/1‑93 (MQ=255)
                                                                 tGGGTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAatat                         <  1:757644/93‑1 (MQ=255)
                                                                  gggTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATc                        <  2:2836784/93‑1 (MQ=255)
                                                                  gggTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATc                        <  2:7391038/93‑1 (MQ=255)
                                                                    gTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCgg                      <  2:3080633/93‑1 (MQ=255)
                                                                    gTTTGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCgg                      <  1:7272387/93‑1 (MQ=255)
                                                                     tttGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGt                     <  1:7398709/93‑1 (MQ=255)
                                                                     tttGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGt                     >  2:2943296/1‑93 (MQ=255)
                                                                     tttGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGt                     <  2:5992451/93‑1 (MQ=255)
                                                                      ttGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGta                    <  1:6491761/93‑1 (MQ=255)
                                                                       tGTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGtat                   <  1:2861860/93‑1 (MQ=255)
                                                                        gTACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGtatt                  <  1:4173399/93‑1 (MQ=255)
                                                                         tACGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGtatta                 <  2:1838600/93‑1 (MQ=255)
                                                                          aCGACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGtattat                <  2:4548678/93‑1 (MQ=255)
                                                                           cGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATc               >  2:75500/1‑93 (MQ=255)
                                                                            gACTCAGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCa              >  2:1636005/1‑93 (MQ=255)
                                                                             aCTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCaa             <  2:6033370/93‑1 (MQ=255)
                                                                               tCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAaga           >  1:7353228/1‑93 (MQ=255)
                                                                               tCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAaga           >  2:1802941/1‑93 (MQ=255)
                                                                                cAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAagag          >  1:3537376/1‑93 (MQ=255)
                                                                                 aGTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAagaga         >  2:799120/1‑93 (MQ=255)
                                                                                  gTAACTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAagagag        >  1:2359032/1‑93 (MQ=255)
                                                                                  gTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAagagag        <  1:6648815/93‑1 (MQ=255)
                                                                                     aCTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCaa     >  2:7174320/1‑93 (MQ=255)
                                                                                     aCTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCaa     <  2:4174392/93‑1 (MQ=255)
                                                                                      cTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCAAc    <  1:7120687/93‑1 (MQ=255)
                                                                                      cTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCAAc    >  2:2406435/1‑93 (MQ=255)
                                                                                      cTGCAGCATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCAAc    >  2:3606615/1‑93 (MQ=255)
                                                                                        gCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCAACCg  <  1:2671296/93‑1 (MQ=255)
                                                                                           |                                                                                         
GGGTTTTATCATTGTTTCGAATTTAATATCAGATACATTATATTGCGATATTTTCCAGTTCCAGATGGGTTTGTACGACTCAGTAACTGCAACATCATTAGCATTATCATTTACTGTCATATTATTAAGTAAATAAACAATACTGCACCAATTAATATCGGTATTATCAAGAGAGCAACCG  >  NC_000913/3629690‑3629870

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