Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,082,696 G→C 11.3% Q22E (CAG→GAG)  fdoH ← formate dehydrogenase O subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,082,6960GC11.3% 150.7 / 13.5 63Q22E (CAG→GAG) fdoHformate dehydrogenase O subunit beta
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base C (3/4);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

GGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACG  >  NC_000913/4082555‑4082838
                                                                                                                                             |                                                                                                                                              
ggTTGTCGTACACCCCAATGTTATAGCCGACGGTATCGCGAATGTCGTTCCAATCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGa                                                                                                                                               <  1:680622/143‑1 (MQ=255)
  ttGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAgg                                                                                                                                             >  2:259190/1‑143 (MQ=255)
       gTACACGCAAATGGTATTGCCGAGGGCATCGGGAATGTCGTTCCACTCTGATCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCgg                                                                                                                                        <  2:789542/143‑1 (MQ=25)
        tACCCCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACTGGTGAAAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                       <  1:847976/143‑1 (MQ=37)
        tACCCCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                       <  1:707801/143‑1 (MQ=255)
        tACACCCCAATGTTACTGCCCCCGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACCGGCTTTTCGTCCCATCCAGGTGGTACCTTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggg                                                                                                                                       <  1:512167/143‑1 (MQ=255)
        tACACCCAAATGTTATTGCCGACGGTATCGCGAAGGTCGTTCCACGCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggg                                                                                                                                       <  2:335826/143‑1 (MQ=255)
        gacgccccAATGCTATTGCCGACGGTATCGCGAATGTCGTTCCCCTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCggg                                                                                                                                       <  2:556307/139‑1 (MQ=255)
         acacCCCAATGTTACTGCCGACGGTATCGCGAATGTCGTTCCACTCCGATCACGCCACCTGACAGGCTTTACAGCCGATACAGGGGGTACCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTCAGGCGCggga                                                                                                                                      <  1:347623/143‑2 (MQ=25)
             cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCa                                                                                                                                  >  1:357978/1‑143 (MQ=255)
             cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGACGATACAGGAGGTAACGTCGATGAGTTTCGGCACTTCTTACTAGAAGTCCCGCGCCTGAGGCGCCGGGGTCa                                                                                                                                  >  2:441339/1‑143 (MQ=37)
             cccAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCAATCTGAACACGCCACCAGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTGCGTGGAAGTACCGCGCCTGAGGCGCGGGGGTcg                                                                                                                                  >  2:503215/1‑142 (MQ=25)
               cAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGa                                                                                                                                >  1:474420/1‑143 (MQ=255)
                aaTGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAc                                                                                                                               >  2:569488/1‑143 (MQ=255)
                     gaTTGCCGACGGTATCGCGAATGTCGTGCCACTCTGAACACGCCGCCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTa                                                                                                                          <  1:384133/142‑1 (MQ=255)
                       tcgCCGACGGTATCGCGTAGGTCGTCGCCCGCCACCCACGCCTCCTCCCAGGCGTAACCGCCGATACAGCTCGTAACGTCGCCGAGTTTCCCCACCTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGt                                                                                                                        <  2:176764/141‑1 (MQ=255)
                             aCGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGCTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCACGGGCCACACCGCTAGTcgagaa                                                                                                                  >  2:672815/1‑139 (MQ=25)
                                     gcgAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGCGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGGCAGACCGTGAGTTGCGGAAggacgac                                                                                                          >  1:779364/1‑142 (MQ=37)
                                              gTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCtt                                                                                                 <  1:259190/143‑1 (MQ=255)
                                                 ccgctctGAACACGCCACCTGACAGGCTTTACAGCCGATACCGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgag                                                                                              <  1:657758/140‑1 (MQ=255)
                                                  catgctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                             <  1:659369/139‑1 (MQ=255)
                                                  cACTCCGAGCTCGCCACCTGACAGGTTTGACAGCCGATACAGGTGGTAGCGTGGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGGGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTgaga                                                                                             <  2:418809/143‑1 (MQ=255)
                                                    ctctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGAtt                                                                                           >  1:486920/1‑143 (MQ=255)
                                                     tctGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTg                                                                                          >  2:565619/1‑143 (MQ=255)
                                                       tGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGAt                                                                                        >  1:326458/1‑143 (MQ=255)
                                                         aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGCTTgctca                                                                                      >  2:837478/1‑139 (MQ=255)
                                                         aaCACGCCACCTGACAGGCTTTACAGCCGATACAGGAGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCAGAGGCGCGGGGGGCAGACCGGTAGTCGCGGAACGACGAACGAGAGCTTGAGATAGATaa                                                                                      >  2:557500/1‑143 (MQ=37)
                                                                cACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGAACGGTAGTCATGGGACGACCAATGATGTCTTGAGAATGATAAGCCatat                                                                               >  2:628930/1‑143 (MQ=37)
                                                                   cTGACAGGCTGTACAGCCGATACAGGTGGTAACGTCGATGAGTGTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATgtc                                                                            <  1:856562/143‑1 (MQ=255)
                                                                        aGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGtctc                                                                       >  1:481914/1‑143 (MQ=255)
                                                                          gCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGGCTCCg                                                                     >  2:801595/1‑143 (MQ=255)
                                                                           cTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGACTCCGt                                                                    >  1:69847/1‑143 (MQ=255)
                                                                            tttACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGAGGGGGGCAGACCGTTAGTCGCGCCACGACGAATGATATCGTGAGAATGGTAGCCCCTATGTCCTCTCCGtt                                                                   >  2:264043/1‑143 (MQ=255)
                                                                                   ccGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCtt                                                            >  2:246507/1‑143 (MQ=255)
                                                                                      aTACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTtctc                                                         >  2:605775/1‑143 (MQ=255)
                                                                                            gtggtACCCTCGATGAGGTGCGCCACTTCTTCCTGGAAGTCCCGCGCCGGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt                                                   <  2:232970/143‑1 (MQ=255)
                                                                                            gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt                                                   >  1:397741/1‑143 (MQ=255)
                                                                                            gtggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACAtt                                                   >  2:38606/1‑143 (MQ=255)
                                                                                              ggtAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTca                                                 <  1:521667/143‑1 (MQ=255)
                                                                                                tAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGCGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAGGCCATATGTCGTCTCCGTTACACCTTGTCCACAtcccca                                               >  2:512232/1‑138 (MQ=255)
                                                                                                 aaCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCTACCTTCACaa                                              >  2:121368/1‑143 (MQ=255)
                                                                                                 aaCGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaa                                              >  1:497564/1‑143 (MQ=255)
                                                                                                 aaCGGCGCTGAGTTTTGCCACTTCTTGCGGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaa                                              <  2:112378/143‑1 (MQ=255)
                                                                                                 aaCGGCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACaa                                              <  1:134418/143‑1 (MQ=255)
                                                                                                          gAGGTTCGCCACTTCTTCCTGGAAGGCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACt                                     <  1:300446/143‑1 (MQ=255)
                                                                                                             tttCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACCAGGCAGGAcctca                                  >  2:680775/1‑139 (MQ=255)
                                                                                                                  ccACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTcc                             >  2:373069/1‑143 (MQ=255)
                                                                                                                      ttcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCg                         >  1:774551/1‑143 (MQ=255)
                                                                                                                       gcttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGt                        <  1:514533/142‑1 (MQ=255)
                                                                                                                         ttcCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGACGTCt                      >  2:112835/1‑143 (MQ=255)
                                                                                                                           ccTGGAAGTCCCGCGCCTGAGGCGTGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACATGTATGGTCTTAAACTCCGGCGTCTGc                    <  2:110107/143‑1 (MQ=255)
                                                                                                                           ccTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGc                    <  1:672815/143‑1 (MQ=255)
                                                                                                                             tGGAAGTCGCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                  <  2:157420/143‑1 (MQ=255)
                                                                                                                             tGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgt                  <  2:326458/143‑1 (MQ=255)
                                                                                                                              ggAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATACGTTGAGATTGATAAGCCATAAGCCGTCTCCGTTACACCTTTTCCCCATTCAAAGGGAAGGACTTAAACTCCAGCGACAGCgtg                 >  2:564544/1‑143 (MQ=255)
                                                                                                                               aatgttccGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCgtgt                >  2:364227/7‑143 (MQ=255)
                                                                                                                                caGTCCCGCGCCGTCTGCGCGGGGGTCATACCGTTAGTCGCGGAACGCCGACCGATCTCCTGATGTTGATATGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGtt               <  2:535015/142‑1 (MQ=255)
                                                                                                                                 agccccGCGCCTGAGGCGCGTGGGTCAGCCCGTTAGTCCCGGATCGACGAACGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTGTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTc              <  2:309087/140‑1 (MQ=255)
                                                                                                                                   tCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGc            <  1:12966/143‑1 (MQ=255)
                                                                                                                                   cccggcgcCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCCTCTGCGGGTTCGc            >  2:646124/5‑143 (MQ=255)
                                                                                                                                   acccGCGCCTGAGGCGCTGGGGTCAGACCGTTAGTCGCCGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGc            <  1:880012/142‑1 (MQ=255)
                                                                                                                                    cccGCGCCTCAGGCACAAGGGTGATACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTACCACATTCACAAGGAAGGACATAAACTCAGCCGTCGCCGTGTTcaca           >  2:111245/1‑143 (MQ=255)
                                                                                                                                     ccGCGCCTCAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATACGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTACACTCCGGCGTCTGCGTGTTCGCAt          >  1:870679/1‑143 (MQ=255)
                                                                                                                                             gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTAACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACg  >  2:7934/1‑143 (MQ=255)
                                                                                                                                             gAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTATTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTCCGTGTTCGCATCACCAACg  >  2:70855/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GGTTGTCGTACACCCCAATGTTATTGCCGACGGTATCGCGAATGTCGTTCCACTCTGAACACGCCACCTGACAGGCTTTACAGCCGATACAGGTGGTAACGTCGATGAGTTTCGCCACTTCTTCCTGGAAGTCCCGCGCCTGAGGCGCGGGGGTCAGACCGTTAGTCGCGGAACGACGAATGATATCTTGAGATTGATAAGCCATATGTCGTCTCCGTTACACCTTTTCCACATTCACAAGGAAGGACTTAAACTCCGGCGTCTGCGTGTTCGCATCACCAACG  >  NC_000913/4082555‑4082838

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: