Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,420,727 G→A 100% G36D (GGC→GAC)  yfcG → GSH‑dependent disulfide bond oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,7270GA98.9% 302.2 / ‑4.5 95G36D (GGC→GAC) yfcGGSH‑dependent disulfide bond oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (48/44);  total (49/46)
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.88e-01

TCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAG  >  NC_000913/2420631‑2420823
                                                                                                |                                                                                                
tctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACgg                                                                                               <  2:2397252/100‑1 (MQ=255)
 ctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAACAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGt                                                                                              >  2:495747/1‑100 (MQ=255)
  tATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTc                                                                                             >  1:1580570/1‑100 (MQ=255)
   aTTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGCAGCAGAGCTGGATTATCCCTTGATTCAGGTAGACCTGGGGTAAGAAGGTc                                                                                             >  2:1485914/1‑99 (MQ=255)
   aTTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCa                                                                                            >  2:1471168/1‑100 (MQ=255)
     ttCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGt                                                                                          >  1:878915/1‑100 (MQ=255)
       cGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGttt                                                                                        <  2:1309970/100‑1 (MQ=255)
        gCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTc                                                                                       >  2:1279024/1‑100 (MQ=255)
        gCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTc                                                                                       <  2:883378/100‑1 (MQ=255)
           ccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTc                                                                                    <  2:1712838/100‑1 (MQ=255)
            cGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTcc                                                                                   >  1:297532/1‑100 (MQ=255)
             gACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCg                                                                                  >  1:1849265/1‑100 (MQ=255)
               cacCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGa                                                                                >  2:1603132/1‑100 (MQ=255)
                 cccAATGGTCACACAAATACGCTGTTTCTCGAAGAAGCAGAGCGGGATTATCGCTTGATTAAGGTAGACCAGGGGAAAGACGGTCAGTTTCGTCCGTAAt                                                                              >  2:188451/1‑100 (MQ=255)
                 cccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAt                                                                              >  2:2086673/1‑100 (MQ=255)
                 cccAATGGACACACAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATCAAGGTAGACCGGGAGAAAGACGGTCAGTTTCGTCCGGAAt                                                                              >  2:995926/1‑100 (MQ=255)
                   cAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAttt                                                                            >  2:1998692/1‑100 (MQ=255)
                   cAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAttt                                                                            >  1:2543677/1‑100 (MQ=255)
                   cAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAttt                                                                            >  1:2501805/1‑100 (MQ=255)
                      tGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTg                                                                         <  2:2098321/100‑1 (MQ=255)
                      tGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTg                                                                         <  1:2514161/100‑1 (MQ=255)
                        gTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcg                                                                       <  2:1252188/100‑1 (MQ=255)
                         tCATAAAATGACTCTGTTTCTCGGAGAAGCAGAGCGGGATTATCGGTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcg                                                                       <  2:1252587/99‑1 (MQ=255)
                         tCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcgc                                                                      <  2:1163787/100‑1 (MQ=255)
                          cacaAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAttttttttt                                                                      >  2:2546365/1‑95 (MQ=255)
                          cacaAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCa                                                                     >  1:1560549/1‑100 (MQ=255)
                          cacaAAAATACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCa                                                                     >  1:1345626/1‑100 (MQ=255)
                            caAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCAtt                                                                   >  1:1546215/1‑100 (MQ=255)
                              aaaTTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTc                                                                 <  1:769496/100‑1 (MQ=255)
                                 ttatgttGTTTCTCGTAGAAGCAGAGCTGGATTATCGCTTGATTTATGTAGACCTGGGGAAAGTCGGTCAGTTTCGTCCGGAATTTTTGCGTATTTCGcc                                                              <  2:216784/94‑1 (MQ=255)
                                 ttACGCTGTTGCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGAATTTCGcc                                                              <  1:700708/100‑1 (MQ=255)
                                  tACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGccc                                                             <  1:431467/100‑1 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                           >  1:1418143/1‑100 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                           >  2:1045138/1‑100 (MQ=255)
                                      cTGTTTCTGGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaaca                                                         <  2:2521257/100‑1 (MQ=255)
                                      cTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaaca                                                         <  1:1964091/100‑1 (MQ=255)
                                      cTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaaca                                                         >  2:57266/1‑100 (MQ=255)
                                       tGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaa                                                        <  2:766819/100‑1 (MQ=255)
                                         tttCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCACTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaaca                                                      <  2:1370689/100‑1 (MQ=255)
                                            ctcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACaaaa                                                   >  1:214219/1‑100 (MQ=255)
                                             tcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAt                                                  >  2:187765/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                                 >  2:1621127/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                                 >  1:1337280/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                                 >  1:650877/1‑100 (MQ=255)
                                               gcagaagCAGAGCAGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTc                                                >  1:1740088/3‑100 (MQ=255)
                                               gaagaaGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTc                                                >  2:464971/1‑100 (MQ=255)
                                                  gaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGTGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCgg                                             <  2:1049435/100‑1 (MQ=255)
                                                   atgCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGc                                            <  2:1557700/98‑1 (MQ=255)
                                                   aagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGc                                            <  2:293851/100‑1 (MQ=255)
                                                   aagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGc                                            <  1:1383565/100‑1 (MQ=255)
                                                     gCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCaa                                          >  1:207301/1‑100 (MQ=255)
                                                        gagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttg                                       <  1:913378/100‑1 (MQ=255)
                                                        gagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttg                                       >  2:1909473/1‑100 (MQ=255)
                                                          gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                     <  2:1928092/100‑1 (MQ=255)
                                                             ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGAt                                  <  2:836884/100‑1 (MQ=255)
                                                             ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATGTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGAt                                  <  2:100549/100‑1 (MQ=255)
                                                               aTTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCa                                >  1:1069457/1‑100 (MQ=255)
                                                                   tCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtct                            <  1:2071616/100‑1 (MQ=255)
                                                                    cGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtctc                           <  2:1024861/100‑1 (MQ=255)
                                                                    cGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtctc                           <  2:940383/100‑1 (MQ=255)
                                                                    cGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtctc                           >  1:1156979/1‑100 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                          >  1:1953440/1‑100 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                          <  1:683682/100‑1 (MQ=255)
                                                                         gATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGcc                      <  2:528257/100‑1 (MQ=255)
                                                                         gATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGcc                      <  1:1667050/100‑1 (MQ=255)
                                                                         gATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGcc                      >  1:2090017/1‑100 (MQ=255)
                                                                           ttAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGa                    <  1:2086687/100‑1 (MQ=255)
                                                                            tAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGAt                   >  2:1096686/1‑100 (MQ=255)
                                                                             aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATg                  <  2:800869/100‑1 (MQ=255)
                                                                             aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATg                  >  2:2480534/1‑100 (MQ=255)
                                                                              aGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATgg                 >  2:2425051/1‑100 (MQ=255)
                                                                                gTAGACCTGGGGAAAGACGGTGAGTTTCGTCCGGAATTTTTGCGCACGTCGCCCAACAACAAGATTCCAGCAATTGTTGATCATTCTCCTGCCGATgg                 <  2:500792/98‑1 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg              <  2:2087887/100‑1 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg              >  2:1542149/1‑100 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg              >  1:1410982/1‑100 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg              <  1:1932988/100‑1 (MQ=255)
                                                                                  aGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcggc             <  1:1017992/100‑1 (MQ=255)
                                                                                   gACCTGGGGAAAGACGGTCAGTTTCTTCCGGAATTTTTGCGCATTTCGCCCAAGCACAAAATTCCGGCAATTGTTGATCACTCTCCTGCCGATGgcggcg            <  2:606788/100‑1 (MQ=255)
                                                                                   gACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGgcggcg            >  2:901039/1‑100 (MQ=255)
                                                                                    aCCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGa           <  2:2539884/100‑1 (MQ=255)
                                                                                    aCCTGGGGAAAGACGGTCAGTTAAGTAAGGAATTTTTGAGCAAATCGACCAACAACAAAAATCCGGCCATTGTTGATCATTCTCCTGCCGATGGCGGCGa           >  1:1493775/1‑100 (MQ=255)
                                                                                     ccTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGaa          >  1:1381972/1‑100 (MQ=255)
                                                                                     ccTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGCGGCGAAc         >  2:2295087/1‑100 (MQ=255)
                                                                                     ccTGGAGAAAGCCGGCCAGTTTCGTCCGCAACTGTTGCGCATTGCGCCCACCAACAACATTCCGGCAATTGTTGATAATTCGCCTGCCGATGGCGGCGaa          <  2:811695/100‑1 (MQ=255)
                                                                                      cTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAc         <  1:901047/100‑1 (MQ=255)
                                                                                       tGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAcc        <  2:2543661/100‑1 (MQ=255)
                                                                                        ggggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCg       >  1:157918/1‑100 (MQ=255)
                                                                                         gggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGc      >  1:1954723/1‑100 (MQ=255)
                                                                                          ggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCt     <  1:768664/100‑1 (MQ=255)
                                                                                          ggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCt     >  1:1145265/1‑100 (MQ=255)
                                                                                           gAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTa    <  2:63506/100‑1 (MQ=255)
                                                                                            aaaGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTaa   <  1:1075189/100‑1 (MQ=255)
                                                                                             aaGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAg  <  1:1046693/100‑1 (MQ=255)
                                                                                             aaGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAg  >  1:2121641/1‑100 (MQ=255)
                                                                                             aaGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAg  >  2:2141464/1‑100 (MQ=255)
                                                                                                |                                                                                                
TCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAG  >  NC_000913/2420631‑2420823

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