Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,435,363 A→G V86A (GTG→GCG)  uspF ← nucleotide binding filament protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,435,3630AG98.4% 227.7 / ‑4.8 63V86A (GTG→GCG) uspFnucleotide binding filament protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (35/27);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

AGCTGCGGCGTTGGAACCGAGCAGATAAGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTAT  >  NC_000913/1435227‑1435499
                                                                                                                                        |                                                                                                                                        
aGCTGCGGCGTTGGAACCGAGCTGATAAGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGctc                                                                                                                                       <  2:58583/140‑1 (MQ=255)
aGCTGCGGCGTTGGAACCGAGCAGATAAGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGctc                                                                                                                                       <  1:597551/140‑1 (MQ=255)
aGCTGCGGCGTTGGAACCGAGCAGATAAGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGctc                                                                                                                                       <  2:820423/140‑1 (MQ=255)
                   aGCAGATAAGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTaaa                                                                                                                    >  2:456748/1‑140 (MQ=255)
                           aGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTaa                                                                                                            <  1:677680/140‑1 (MQ=255)
                              ggtgATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGa                                                                                                         <  1:845061/140‑1 (MQ=255)
                                  atatCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGAtctc                                                                                                     >  2:987159/1‑140 (MQ=255)
                                  atatCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGAtctc                                                                                                     >  1:899389/1‑140 (MQ=255)
                                  atatCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGAtctc                                                                                                     >  2:658573/1‑140 (MQ=255)
                                  atatCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGAtctc                                                                                                     >  2:76783/1‑140 (MQ=255)
                                       cGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCa                                                                                                >  1:1024890/1‑140 (MQ=255)
                                       cGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCa                                                                                                >  2:731446/1‑140 (MQ=255)
                                       cGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCa                                                                                                >  1:376150/1‑140 (MQ=255)
                                          tCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGtt                                                                                             <  1:405249/140‑1 (MQ=255)
                                                       atgatgATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTc                                                                                <  1:722802/140‑1 (MQ=255)
                                                        tgatgaTCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTcc                                                                               >  2:839994/1‑140 (MQ=255)
                                                        tgatgaTCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTcc                                                                               >  1:612577/1‑140 (MQ=255)
                                                          atgatCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGc                                                                             >  2:125842/1‑140 (MQ=255)
                                                            gatCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCtt                                                                           <  2:604278/140‑1 (MQ=255)
                                                             atCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCACttt                                                                          <  2:816544/140‑1 (MQ=255)
                                                              tCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTc                                                                         <  2:561400/140‑1 (MQ=255)
                                                               cATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCa                                                                        >  2:177829/1‑140 (MQ=255)
                                                               cATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCa                                                                        >  1:314142/1‑140 (MQ=255)
                                                                atatGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAg                                                                       >  1:942197/1‑140 (MQ=255)
                                                                atatGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAg                                                                       >  1:423207/1‑140 (MQ=255)
                                                                      gCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTcc                                                                               >  1:293674/1‑126 (MQ=255)
                                                                      gCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTcc                                                                               <  2:293674/126‑1 (MQ=255)
                                                                         gggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAt                                                              >  2:978893/1‑140 (MQ=255)
                                                                         gggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAt                                                              >  1:632555/1‑140 (MQ=255)
                                                                         gggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAt                                                              >  2:1006695/1‑140 (MQ=255)
                                                                         gggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAt                                                              >  2:935473/1‑140 (MQ=255)
                                                                         gggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAt                                                              >  1:459562/1‑140 (MQ=255)
                                                                          ggATCTTCTTCGCTAATTCAAAAAAGCAGTCTTTAGGCGAGTCTTCTTCAACATGGACATGCTCTCTGTCGGTTGGCAGTTTTAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAtt                                                             <  2:1018951/140‑1 (MQ=255)
                                                                          ggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAtt                                                             >  1:712427/1‑140 (MQ=255)
                                                                          ggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAtt                                                             <  2:376150/140‑1 (MQ=255)
                                                                          ggATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCAtt                                                             <  2:321086/140‑1 (MQ=255)
                                                                                 cttcGCCAATTCCAGAATGCGGTATTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTa                                                      <  1:76783/140‑1 (MQ=255)
                                                                                    cGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAAtt                                                   <  2:734330/140‑1 (MQ=255)
                                                                                    cGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAAtt                                                   >  1:474723/1‑140 (MQ=255)
                                                                                       cAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTg                                                <  1:987159/140‑1 (MQ=255)
                                                                                         aTTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCg                                              >  2:108204/1‑140 (MQ=255)
                                                                                           tCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGa                                            <  2:423207/140‑1 (MQ=255)
                                                                                            ccAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGaa                                           <  1:528065/140‑1 (MQ=255)
                                                                                                aaTGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACg                                       >  2:470746/1‑140 (MQ=255)
                                                                                                 aTGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGc                                      <  2:896500/140‑1 (MQ=255)
                                                                                                  tGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCt                                     <  1:456748/140‑1 (MQ=255)
                                                                                                    cGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTaa                                   >  1:622669/1‑140 (MQ=255)
                                                                                                       tCTTTGGGCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAAcc                                >  2:429108/1‑140 (MQ=255)
                                                                                                            gggCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCaga                            <  2:846713/139‑1 (MQ=255)
                                                                                                            gggCGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCaga                            >  1:846713/1‑139 (MQ=255)
                                                                                                               cGAGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAgg                        <  1:420533/140‑1 (MQ=255)
                                                                                                                 aGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCa                      >  1:90478/1‑140 (MQ=255)
                                                                                                                 aGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCa                      >  1:886975/1‑140 (MQ=255)
                                                                                                                 aGCCTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCa                      >  2:96198/1‑140 (MQ=255)
                                                                                                                   ccTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCAta                    >  1:413612/1‑140 (MQ=255)
                                                                                                                   ccTTCCTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCAta                    >  1:246133/1‑140 (MQ=255)
                                                                                                                       ccTCAACATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCAtagtag                <  2:314142/140‑1 (MQ=255)
                                                                                                                             cATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGt          <  1:658573/140‑1 (MQ=255)
                                                                                                                             cATGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGt          <  1:978893/140‑1 (MQ=255)
                                                                                                                               tGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGa        <  1:632521/140‑1 (MQ=255)
                                                                                                                               tGGACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGa        <  2:459562/140‑1 (MQ=255)
                                                                                                                                 gACATGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAg      >  1:688029/1‑140 (MQ=255)
                                                                                                                                     tGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTAt  <  1:125842/140‑1 (MQ=255)
                                                                                                                                     tGCGCTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTAt  >  2:282678/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
AGCTGCGGCGTTGGAACCGAGCAGATAAGTGGTGATATCCGGTCGATGGGAAGCAATGATGATCATATGAGCGGGGATCTTCTTCGCCAATTCCAGAATGCGGTCTTTGGGCGAGCCTTCCTCAACATGGACATGCACTCTGTCGGTTGGCAGTTTAAATTTTTTAATGATCTCTTCCAGTTGCGATTTGGCTTCCGCTTTCAGGTCATCCATTGCCGGTAATTCTGCGGAATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTAT  >  NC_000913/1435227‑1435499

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