Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,654 A→G Y98H (TAC→CAC) 
Y98H (TAC→CAC) 
fdhF ←
b4560 ←
selenopolypeptide subunit of formate dehydrogenase H
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,6540AG95.5% 517.3 / 9.5 179Y98H (TAC→CAC) 
Y98H (TAC→CAC) 
fdhF
b4560
selenopolypeptide subunit of formate dehydrogenase H
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base G (69/102);  total (73/106)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 185 total aligned reads displayed.

ACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGA  >  NC_000913/4296566‑4296742
                                                                                        |                                                                                        
aCGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTtc                                                                                      >  1:3866765/1‑93 (MQ=255)
 cGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTtct                                                                                     <  2:2078167/93‑1 (MQ=255)
  ggCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTtctc                                                                                    <  1:2705295/93‑1 (MQ=255)
  ggCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTtctc                                                                                    <  1:3447706/93‑1 (MQ=255)
  ggCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTtctc                                                                                    >  2:2402285/1‑93 (MQ=255)
      cgcgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTg                                                                                >  2:1147882/1‑93 (MQ=255)
      cgcgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTg                                                                                <  2:748560/93‑1 (MQ=255)
       gcgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGa                                                                               <  1:1625609/93‑1 (MQ=255)
        cgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGAt                                                                              <  1:114897/93‑1 (MQ=255)
        cgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGAt                                                                              <  1:1172855/93‑1 (MQ=255)
        cgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCATGCTTCTCTTTGAt                                                                              <  2:781188/93‑1 (MQ=255)
        cgcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGCAGTCTGGATGGCATCCGGACCGTGGTTCTCTTGGAt                                                                              >  1:1369097/1‑93 (MQ=255)
         gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATg                                                                             <  1:1825653/93‑1 (MQ=255)
         gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATg                                                                             >  1:3002458/1‑93 (MQ=255)
         gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATg                                                                             <  2:1931937/93‑1 (MQ=255)
         gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATg                                                                             <  2:1468351/93‑1 (MQ=255)
         gcAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATg                                                                             <  1:2620793/93‑1 (MQ=255)
           aaaTTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGgc                                                                           >  1:3954282/1‑93 (MQ=255)
           aaaTTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGgc                                                                           <  1:921830/93‑1 (MQ=255)
              tttttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCt                                                                        >  1:987625/1‑93 (MQ=255)
               ttttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTc                                                                       >  2:3049832/1‑93 (MQ=255)
                tttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCa                                                                      >  1:4489065/1‑93 (MQ=255)
                tttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCa                                                                      <  1:3765560/93‑1 (MQ=255)
                tttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCa                                                                      <  1:523996/93‑1 (MQ=255)
                tttGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCa                                                                      <  1:3632410/93‑1 (MQ=255)
                   tcATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGgc                                                                   <  1:2204048/92‑1 (MQ=255)
                    cATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGgcg                                                                  >  1:311164/1‑93 (MQ=255)
                       tACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTc                                                               <  1:3261647/93‑1 (MQ=255)
                       tACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTc                                                               >  2:4019147/1‑93 (MQ=255)
                       tACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTc                                                               <  2:2546624/93‑1 (MQ=255)
                       tACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTc                                                               >  1:3140810/1‑93 (MQ=255)
                         cATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACAGTGCTTCTCTTTGATGGCGCTCAGGCGCTCgg                                                             >  2:482970/1‑93 (MQ=255)
                          aTAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGc                                                            >  1:4485533/1‑93 (MQ=255)
                          aTAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGc                                                            <  2:3702639/93‑1 (MQ=255)
                          aTAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGc                                                            >  2:837293/1‑93 (MQ=255)
                           tAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCa                                                           >  1:3387118/1‑93 (MQ=255)
                           tAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCGCTTTGATGGCGCTCAGGCGCTCGGCa                                                           >  2:3796780/1‑93 (MQ=255)
                           tAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCa                                                           <  2:3340530/93‑1 (MQ=255)
                            aGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCaa                                                          >  2:3115040/1‑93 (MQ=255)
                            aGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCaa                                                          <  2:4486844/93‑1 (MQ=255)
                             gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAAc                                                         >  2:1110902/1‑93 (MQ=255)
                             gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAAc                                                         >  1:4667705/1‑93 (MQ=255)
                             gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAAc                                                         >  2:1648566/1‑93 (MQ=255)
                             gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAAc                                                         <  2:498377/93‑1 (MQ=255)
                             gTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAAc                                                         >  2:3030002/1‑93 (MQ=255)
                                ggTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTa                                                      >  1:3787091/1‑93 (MQ=255)
                                 gTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTaa                                                     <  2:3283367/93‑1 (MQ=255)
                                 gTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTaa                                                     >  1:2924406/1‑93 (MQ=255)
                                  tttCGTTCCCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAt                                                    <  2:1469372/93‑1 (MQ=255)
                                  tttCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAt                                                    <  2:1871550/93‑1 (MQ=255)
                                   ttCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAAtt                                                   <  1:2529237/93‑1 (MQ=255)
                                    tCGTTCCCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTc                                                  <  2:3018205/93‑1 (MQ=255)
                                    tCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTc                                                  <  2:853609/93‑1 (MQ=255)
                                    tCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTc                                                  <  1:3678390/93‑1 (MQ=255)
                                     cGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGCCCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCa                                                 <  1:4436125/93‑1 (MQ=255)
                                      gTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAg                                                >  1:2347663/1‑93 (MQ=255)
                                         aCCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCTGGCGATCGGCAACGTAATTCAGTGc                                             >  1:3806888/1‑93 (MQ=255)
                                          cccGTACCCCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGcc                                            >  2:4580358/1‑93 (MQ=255)
                                            cGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTc                                          >  1:1208254/1‑93 (MQ=255)
                                             gTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCa                                         <  1:16866/93‑1 (MQ=255)
                                             gTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCa                                         <  1:2421370/93‑1 (MQ=255)
                                               aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATc                                       >  1:960357/1‑93 (MQ=255)
                                               aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATc                                       <  2:2298483/93‑1 (MQ=255)
                                               aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATc                                       >  1:330611/1‑93 (MQ=255)
                                               aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATc                                       <  2:3891048/93‑1 (MQ=255)
                                               aCCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATc                                       <  1:1536242/93‑1 (MQ=255)
                                                  aCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCa                                    <  2:3906946/93‑1 (MQ=255)
                                                  aCGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGCAATTCAGTGTCTCATCCCa                                    >  2:1280822/1‑93 (MQ=255)
                                                    gcgAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAgg                                  >  1:2091789/1‑93 (MQ=255)
                                                       aggagCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGaaa                               >  2:3350442/1‑93 (MQ=255)
                                                       aggagCCGGCCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGaaa                               >  1:3369501/1‑93 (MQ=255)
                                                         gagCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACa                             <  1:1415032/93‑1 (MQ=255)
                                                         gagCCGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACa                             <  1:2820790/93‑1 (MQ=255)
                                                             cGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGtt                         <  2:4112608/93‑1 (MQ=255)
                                                             cGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGtt                         >  1:4310969/1‑93 (MQ=255)
                                                             cGGTCGTCTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGtt                         <  1:2208486/93‑1 (MQ=255)
                                                               gtcgtcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCg                       <  2:1794657/93‑1 (MQ=255)
                                                                tcgtcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                      <  2:1355091/93‑1 (MQ=255)
                                                                tcgtcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                      <  1:4148264/93‑1 (MQ=255)
                                                                tcgtcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                      >  2:3386084/1‑93 (MQ=255)
                                                                tcgtcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                      <  2:358756/93‑1 (MQ=255)
                                                                tcgtcTGGATGGCATCCGGACCGTGCTGCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                      <  1:668489/93‑1 (MQ=255)
                                                                tcgtcTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGa                      >  1:1925922/1‑93 (MQ=255)
                                                                 cgtcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAg                     <  2:2478618/93‑1 (MQ=255)
                                                                   tcTGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGtt                   <  2:2006989/93‑1 (MQ=255)
                                                                     tGGATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTg                 >  1:1034467/1‑93 (MQ=255)
                                                                       gATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTgcc               <  1:2853270/93‑1 (MQ=255)
                                                                       gATGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTgcc               <  2:1471427/93‑1 (MQ=255)
                                                                        aTGGCATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTgccg              <  1:2001011/93‑1 (MQ=255)
                                                                            cATCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCAc          <  2:583549/93‑1 (MQ=255)
                                                                              tCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGc        <  2:4558265/93‑1 (MQ=255)
                                                                              tCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGc        <  1:40395/93‑1 (MQ=255)
                                                                              tCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGc        <  1:620531/93‑1 (MQ=255)
                                                                              tCCGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGc        <  2:2226869/93‑1 (MQ=255)
                                                                               ccGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCt       <  1:399302/93‑1 (MQ=255)
                                                                                cGGACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTg      <  1:4543170/93‑1 (MQ=255)
                                                                                 ggACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTgg     >  1:4196911/1‑93 (MQ=255)
                                                                                 ggACCGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTgg     >  1:856679/1‑93 (MQ=255)
                                                                                    ccGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGa  <  1:4324561/93‑1 (MQ=255)
                                                                                    ccGTGCTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGa  >  2:2717170/1‑93 (MQ=255)
                                                                                        |                                                                                        
ACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGA  >  NC_000913/4296566‑4296742

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