Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 111,615 T→C intergenic (+182/+34) mutT → / ← yacG 7,8‑dihydro‑8‑oxoguanine‑triphosphatase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,6150TC85.1% 700.0 / 115.1 303intergenic (+182/+34)mutT/yacG7,8‑dihydro‑8‑oxoguanine‑triphosphatase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/19);  new base C (140/117);  total (167/136)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 312 total aligned reads displayed.

CTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGA  >  NC_000913/111525‑111707
                                                                                          |                                                                                            
cTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCt                                                                                            >  2:956681/1‑93 (MQ=255)
cTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCt                                                                                            >  2:7007917/1‑93 (MQ=255)
 tCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCtt                                                                                           >  2:6751766/1‑93 (MQ=255)
   aTTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTAt                                                                                         >  1:5109834/1‑93 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        >  2:365313/1‑93 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        >  2:31207/1‑93 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        <  1:7062688/93‑1 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        >  2:5579830/1‑93 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        >  1:6498717/1‑93 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        >  1:5222687/1‑93 (MQ=255)
    ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATc                                                                                        <  1:1712535/93‑1 (MQ=255)
     tGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCa                                                                                       >  2:44361/1‑93 (MQ=255)
     tGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCa                                                                                       >  1:101692/1‑93 (MQ=255)
      gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAg                                                                                      >  1:4878832/1‑93 (MQ=255)
      gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAg                                                                                      <  1:5016215/93‑1 (MQ=255)
      gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAg                                                                                      <  1:6954240/93‑1 (MQ=255)
      gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAg                                                                                      >  2:793242/1‑93 (MQ=255)
       aTAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAgg                                                                                     <  1:1288119/93‑1 (MQ=255)
       aTAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAgg                                                                                     >  1:3320731/1‑93 (MQ=255)
          aaTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCt                                                                                  >  2:323668/1‑93 (MQ=255)
           aTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTa                                                                                 >  2:1570315/1‑93 (MQ=255)
           aTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTa                                                                                 <  1:5572912/93‑1 (MQ=255)
           aTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTa                                                                                 >  1:1505561/1‑93 (MQ=255)
            tttGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTAc                                                                                >  2:1430249/1‑93 (MQ=255)
              tGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGt                                                                              <  2:5888643/93‑1 (MQ=255)
                tttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGc                                                                            >  1:4626648/1‑93 (MQ=255)
                tttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGc                                                                            <  2:1322097/93‑1 (MQ=255)
                 ttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCg                                                                           <  1:322883/93‑1 (MQ=255)
                    tGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCa                                                                        >  2:2813193/1‑93 (MQ=255)
                    tGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCa                                                                        >  1:4837672/1‑93 (MQ=255)
                      tAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATc                                                                      <  1:965662/93‑1 (MQ=255)
                        ggTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAg                                                                    >  2:7319550/1‑93 (MQ=255)
                          tCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGAc                                                                  <  2:218992/93‑1 (MQ=255)
                           cGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACa                                                                 <  1:2127901/93‑1 (MQ=255)
                            ggATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACaa                                                                >  1:4947514/1‑93 (MQ=255)
                             gATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                               >  1:5312082/1‑93 (MQ=255)
                             gATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                               <  1:4546181/93‑1 (MQ=255)
                             gATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                               >  1:2756589/1‑93 (MQ=255)
                              aTAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAAt                                                              >  2:914749/1‑93 (MQ=255)
                               tAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATg                                                             >  1:2366433/1‑93 (MQ=255)
                                 aGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTc                                                           >  1:1783450/1‑93 (MQ=255)
                                   gCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTGCGTGCGGCATCAGACAAATGTCAc                                                         <  2:473059/93‑1 (MQ=255)
                                     gTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTg                                                       >  2:1182469/1‑93 (MQ=255)
                                      ttCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGc                                                      >  1:43798/1‑93 (MQ=255)
                                       tCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCt                                                     >  2:3287240/1‑93 (MQ=255)
                                        cACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCtt                                                    >  2:4559541/1‑93 (MQ=255)
                                         aCGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTGCGTGCGGCATCAGACAAATGTCACTGCttt                                                   <  2:2027205/93‑1 (MQ=255)
                                         aCGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                                   <  1:5866374/93‑1 (MQ=255)
                                         aCGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                                   >  2:508556/1‑93 (MQ=255)
                                          cgccgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTg                                                  >  2:6312160/1‑93 (MQ=255)
                                             cgcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGtt                                               >  1:1411140/1‑93 (MQ=255)
                                              gcATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttc                                              >  1:631138/1‑93 (MQ=255)
                                                aTCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttctt                                            >  2:6904595/1‑93 (MQ=255)
                                                 tCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                                           >  1:7256779/1‑93 (MQ=255)
                                                 tCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                                           <  2:3761217/93‑1 (MQ=255)
                                                 tCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                                           >  2:1158521/1‑93 (MQ=255)
                                                 tCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                                           >  2:7369/1‑93 (MQ=255)
                                                    gACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCgnn                                        >  1:3268356/1‑91 (MQ=255)
                                                    gACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCt                                        >  2:5774994/1‑93 (MQ=255)
                                                     aCATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTc                                       >  2:2903284/1‑93 (MQ=255)
                                                      cATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                                      <  1:6765102/93‑1 (MQ=255)
                                                       aTTTGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCa                                     >  1:5165422/1‑93 (MQ=255)
                                                        tttGCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAg                                    <  1:369992/93‑1 (MQ=255)
                                                           gCACAAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCa                                 <  2:4360724/93‑1 (MQ=255)
                                                             acaAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATc                               <  1:5889499/93‑1 (MQ=255)
                                                             acaAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATc                               >  1:1885585/1‑93 (MQ=255)
                                                             acaAGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATc                               <  1:3987129/93‑1 (MQ=255)
                                                              caAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCg                              <  2:1795847/93‑1 (MQ=255)
                                                                aGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCt                            >  2:1062451/1‑93 (MQ=255)
                                                                aGATGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCt                            >  1:5532724/1‑93 (MQ=255)
                                                                  aTGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCttt                          <  1:716146/93‑1 (MQ=255)
                                                                   tGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                         >  2:7218693/1‑93 (MQ=255)
                                                                   tGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                         >  1:7417626/1‑93 (MQ=255)
                                                                   tGCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                         >  1:4060805/1‑93 (MQ=255)
                                                                    gCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCg                        <  1:1198206/93‑1 (MQ=255)
                                                                    gCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCg                        >  1:1467510/1‑93 (MQ=255)
                                                                    gCCTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCg                        <  2:2681023/93‑1 (MQ=255)
                                                                     ccTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCgg                       <  2:6563539/93‑1 (MQ=255)
                                                                     ccTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCgg                       >  2:5986592/1‑93 (MQ=255)
                                                                      cTGATGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGa                      >  2:5268912/1‑93 (MQ=255)
                                                                       tGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGACACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaa                     >  1:771904/1‑93 (MQ=255)
                                                                         aTGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                   >  2:3924744/1‑93 (MQ=255)
                                                                          tGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGTTTTCGGAAAGa                  >  1:3123519/1‑93 (MQ=255)
                                                                          tGCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                  <  1:1998897/93‑1 (MQ=255)
                                                                           gCGACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGAt                 <  2:5305502/93‑1 (MQ=255)
                                                                             gACGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCg               >  2:3178043/1‑93 (MQ=255)
                                                                              aCGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc              <  2:3594935/93‑1 (MQ=255)
                                                                              aCGCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc              <  2:4142301/93‑1 (MQ=255)
                                                                                gCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa            <  1:6111110/93‑1 (MQ=255)
                                                                                gCTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa            <  2:5149500/93‑1 (MQ=255)
                                                                                 cTGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAc           <  1:6859636/93‑1 (MQ=255)
                                                                                  tGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAct          >  1:5776645/1‑93 (MQ=255)
                                                                                  tGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAct          <  2:1713880/93‑1 (MQ=255)
                                                                                  tGTCCGCGCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAct          <  1:6044163/93‑1 (MQ=255)
                                                                                      cgcgCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActgctg      >  1:3791536/1‑93 (MQ=255)
                                                                                      cgcgCCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActgctg      >  1:3033901/1‑93 (MQ=255)
                                                                                       gcgcCTTATCAGGCCTACGTGCGGCATCAGACAATTGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTgg     <  1:2409871/93‑1 (MQ=255)
                                                                                       gcgcCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTgg     <  2:4667908/93‑1 (MQ=255)
                                                                                       gcgcCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTgg     <  2:1940062/93‑1 (MQ=255)
                                                                                       gcgTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTgg     >  1:2687716/1‑93 (MQ=255)
                                                                                          |                                                                                            
CTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGA  >  NC_000913/111525‑111707

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