Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,142,885 A→G 8.1% G35G (GGT→GGC)  hybC ← hydrogenase 2 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,142,8850AG 8.1% 149.0 / 2.8 65G35G (GGT→GGC) hybChydrogenase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base G (4/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTT  >  NC_000913/3142750‑3143026
                                                                                                                                       |                                                                                                                                             
cAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTCCATCCGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAg                                                                                                                                        >  1:512966/1‑143 (MQ=255)
 aGCGCACTTTCTGCCGCACGACCGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCGTCCATGCATCGCACCGATCCCGGTTTTTCCCGACGTCTTCAACACTCCGCCACTTGGTACCagcagc                                                                                                                                       >  2:825576/1‑143 (MQ=255)
  gcgcACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCGCTTCCATGCCGGGCCACATCGTACCGGACGcc                                                                                                                                      >  2:745914/1‑143 (MQ=255)
   cgcACTGTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATGCGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGccc                                                                                                                                     <  2:145977/143‑1 (MQ=255)
    ccACTTTCCGCCGCACGAACGGAAGACCGCGCGTGAGCGGCAGTACATACGCCACAGCTACGTTGCACAATCCTCCATGCATCGCGCGGATCGCGGGTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCa                                                                                                                                    <  1:786159/142‑1 (MQ=255)
     cACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCAt                                                                                                                                   <  1:726452/143‑1 (MQ=255)
      gcTTTCTGCCGCACGAGCCATAGACAGCGCTTTGGTTGTAGTATATACGCCACTGATGCGTTGCACAATCATCCTTGTATCGCGCGGATCGCGGTTTTTCACGCTCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATg                                                                                                                                  <  2:771203/142‑1 (MQ=255)
      aCTTTCTGCCGCACGAACGGAAGACAGCGCGTGCGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCAGGCATCGCGCGGATCGCGGTTTTTCAGGATCTCTTCCATCCCGAGCCACAGGGTGCCGGAAGCCCAtc                                                                                                                                  >  2:869110/1‑142 (MQ=255)
      aCTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGCTCTTCACGATCTCTTCCATGCAGCGCCACACGGCACCGGAAGcgcgtg                                                                                                                                  >  2:561769/1‑138 (MQ=255)
           cTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCCTGCATCGCGCGGATCGCGGTTTTTCCCGCTCTCTTCCATGCCGCGCCCCATGGTACCGGAAGCCCAtccctcc                                                                                                                             >  1:674242/1‑137 (MQ=255)
              ccGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTGCCGGAAGCCCATGCTTTCGaa                                                                                                                          >  1:405204/1‑143 (MQ=255)
              ccGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTGCCGGAAGCCCATGCTTTCGaa                                                                                                                          >  1:879921/1‑143 (MQ=255)
              ccGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGaa                                                                                                                          >  2:745614/1‑143 (MQ=255)
               cGCACGAACGGAAGCCAGCGCGTGGGTGGGAGGACATCCGCCACAGATACGGTGCACAATCAGCCATGCAGCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGaaa                                                                                                                         <  1:100134/143‑1 (MQ=255)
                   cGAGCGGAAGACAGCGCGTGAGGGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAacgac                                                                                                                     <  2:37198/143‑1 (MQ=255)
                    gAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAacgacg                                                                                                                    >  1:544617/1‑143 (MQ=255)
                    gAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGGTTTTCACGATCTCTTCCATGCCGCGCCGCATGGTACCGGTAGCCGATGCTTTCTAAtcggcg                                                                                                                    >  1:448089/1‑143 (MQ=255)
                      aCGGAAGACGGCGCGTGCGTGGTGGGACAACCCTTGTAGAGAGGTCGCACAATCATCCATGCCTCGCGCGGATCGCGGTCTTTCACGATCTCGTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGcc                                                                                                                  <  2:865011/143‑1 (MQ=255)
                         gAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCAACCATGCATCACGCGGATCGCGGTGTTTCACGATCTCTTCCATGCAGCGCCACCGGGTACCGGAAACACTTTCATTCCAAACGACGCCAtt                                                                                                               >  2:97485/1‑143 (MQ=255)
                            gACCGCGCGTGAGTGGTAGTACATACGCCACAGATACGTCGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTc                                                                                                            <  1:357436/143‑1 (MQ=255)
                             aCAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCg                                                                                                           >  1:773077/1‑143 (MQ=255)
                               aGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGAt                                                                                                         <  1:825414/143‑1 (MQ=255)
                                        gTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATc                                                                                                >  2:619125/1‑143 (MQ=255)
                                         tGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCg                                                                                               >  1:638569/1‑143 (MQ=255)
                                          ggTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCACTCGa                                                                                              >  2:367210/1‑143 (MQ=255)
                                            tagtaCATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATg                                                                                            <  2:542221/143‑1 (MQ=255)
                                            gagcaCATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTCCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATg                                                                                            <  1:20952/139‑1 (MQ=255)
                                                aCATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTa                                                                                        <  1:300740/143‑1 (MQ=255)
                                                  aTACGCCACAGATACGTTGCACAATCATCCATTCATCTCGCGCATCGTGGTTTTTCACGCTCCCTCCCACGCCGCACCACATGGTACCGTAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTaaa                                                                                      <  2:210268/143‑1 (MQ=255)
                                                       ccACGGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGTTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAACTGccc                                                                                 <  2:835129/143‑1 (MQ=255)
                                                       ccACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGccc                                                                                 >  1:667041/1‑143 (MQ=255)
                                                             aTACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAAt                                                                           <  2:366375/143‑1 (MQ=255)
                                                                 gTTGCACAATCATCCATGCATCGCGCGGATCGCGGTGTTTCACGATCTCTTCCATGCAGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGGACAAATGCCATCCAATACga                                                                       >  2:111758/1‑142 (MQ=255)
                                                                  ttGCACACTCATCCATGCATCGCGCGGATCGCGGTTTTCCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACggg                                                                      <  2:711359/143‑1 (MQ=255)
                                                                  ttGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACggg                                                                      <  2:232296/143‑1 (MQ=255)
                                                                        aaTCATCCATGCATCGCGCGGACCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCg                                                                <  1:760167/143‑1 (MQ=255)
                                                                          tcatcCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCgca                                                              >  1:217869/1‑141 (MQ=255)
                                                                           catcCCTGCATCGCGCGGATCGCGGTTTTTCACGCTCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGAt                                                             <  1:718284/143‑1 (MQ=255)
                                                                             tcCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCa                                                           >  2:101412/1‑143 (MQ=255)
                                                                               cgtGCATCGCGCGGATCGCGGTTTTTCACGATCTCTGCCATGCCGCGCCACATGGTACCGGAAGCCCAGGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAAt                                                         <  2:108611/141‑1 (MQ=255)
                                                                                   cATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCCTTTTCGATCTCGCAATCGATGCCTCATAGCCCCTCCACACGCGTTACCGCAACAATAGTa                                                     >  2:183436/1‑143 (MQ=255)
                                                                                         gcgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGCCGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATtctc                                               >  1:698303/1‑143 (MQ=255)
                                                                                         gcgGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCTGATCAAGAGTACTtctc                                               >  2:563956/1‑143 (MQ=255)
                                                                                           ggATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCACATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATAGATGCGGAAATGCCCCTCAATACGGGTTATCGGAGGAGTAGTAGttgagtg                                             >  2:444622/1‑138 (MQ=255)
                                                                                           ggATCGCGGTTTTTCACGATCTATTCCATGCCGCGCCACATGGTGCCGGAAGCACATGCTTTCCAAACGACGCTATTTTCGATTTCGCAATCGATGTGTTAGTGCAACGCAATGCGTGGTGCCGAACCAGTAGTAATtccgtg                                             >  2:100304/1‑139 (MQ=255)
                                                                                              tCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGCTCTATAGTGATACCCTGGCt                                          >  1:501506/1‑143 (MQ=255)
                                                                                                 cgGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCAt                                       <  2:773077/143‑1 (MQ=255)
                                                                                                 cgGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCAt                                       >  1:691848/1‑143 (MQ=255)
                                                                                                  ggTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTGCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATg                                      <  2:405204/143‑1 (MQ=255)
                                                                                                  ggTTTGTCACGATCTCTTCCATGCCGCGCCACATGGTGCCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGAATCAATAGTAATTCTCTGGCTCATg                                      <  2:879921/143‑1 (MQ=255)
                                                                                                      ttgcccGATCTCTTCGATGCCGCGGCACAGGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGAGTTCGCAAGCGATGCGTAAATGCCCCGCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCttt                                  <  2:556195/138‑1 (MQ=255)
                                                                                                      gttCACGATCTGTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGAGTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCttt                                  <  2:169396/142‑1 (MQ=255)
                                                                                                            gATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCACTACGGGTTACCGGATCAATAGTAATTCTCTGGCTCCTGCTTTGCTCGc                            >  1:829477/1‑143 (MQ=255)
                                                                                                                  ttCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATg                      >  2:119316/1‑143 (MQ=255)
                                                                                                                  ttCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAACTCTCTGGCTCATGCCTAGCTCGcccactg                      >  2:549384/1‑138 (MQ=255)
                                                                                                                   tCCATGCCGCGCCACGGGGTACCGGAAGCCCAGGCTTTCGAAACGCCGCCGTTTTCGATTTCGCGATCGATGCATAAATGCCCCTCAATACGGGTTACCGGATTAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGa                     <  1:486818/143‑1 (MQ=255)
                                                                                                                      ttGCCGCGCCACTTGTTACGGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGa                  <  2:117191/142‑1 (MQ=255)
                                                                                                                      atccctcgcgACATGGTACCTGAACCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCATTACGGGTTGCCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGa                  <  1:167542/137‑1 (MQ=255)
                                                                                                                      aTGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTCACCGGATCAATAGGAATTCTCTGAGTCAAGCTTTGCTGGCCCCATGACgg                  >  1:648019/1‑142 (MQ=255)
                                                                                                                             gcCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGAt           <  2:570779/143‑1 (MQ=255)
                                                                                                                              ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGTTTTGCTCGCCTCATGACGATTATGATc          <  1:572395/143‑1 (MQ=255)
                                                                                                                              ccACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATc          <  1:315338/143‑1 (MQ=255)
                                                                                                                                    ggTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttg    <  1:745614/143‑1 (MQ=255)
                                                                                                                                      tACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttgtt  <  2:206767/143‑1 (MQ=255)
                                                                                                                                      tACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCCCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCgttttgtt  <  1:619125/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                             
CAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCACAATCATCCATGCATCGCGCGGATCGCGGTTTTTCACGATCTCTTCCATGCCGCGCCACATGGTACCGGAAGCCCATGCTTTCGAAACGACGCCATTTTCGATTTCGCAATCGATGCGTAAATGCCCCTCAATACGGGTTACCGGATCAATAGTAATTCTCTGGCTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTT  >  NC_000913/3142750‑3143026

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