Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,925 +A coding (35/1431 nt) yaaJ ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,9221.A97.1% 212.5 / ‑0.7 68coding (38/1431 nt)yaaJputative transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (31/35);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTC  >  NC_000913/7791‑8054
                                                                                                                                     |                                                                                                                                     
cGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aa‑aagga                                                                                                                               <  2:305201/139‑5 (MQ=255)
  ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACg                                                                                                                              <  2:805360/139‑1 (MQ=255)
  ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACg                                                                                                                              >  1:351324/1‑139 (MQ=255)
   cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGc                                                                                                                             >  1:472319/1‑139 (MQ=255)
   cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGc                                                                                                                             >  2:156303/1‑139 (MQ=255)
           ggATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATg                                                                                                                     >  2:722343/1‑139 (MQ=255)
            gATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGa                                                                                                                    >  2:705582/1‑139 (MQ=255)
             aTGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGaa                                                                                                                   <  1:156303/139‑1 (MQ=255)
                aaTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGAGTAACTGCACAAATCCGGTGCGAAAAGTGAACCAACGACCTGCGCCGATGATCAGGTAAATCATTACCGTTCCGCA‑AA‑AGGACGCTGTTAaggaagga                                                                                                                >  1:373963/1‑139 (MQ=255)
                            aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCGACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    <  1:722343/139‑1 (MQ=255)
                            aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    <  2:30214/139‑1 (MQ=255)
                                 ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCt                                                                                                                            >  2:1008292/1‑110 (MQ=255)
                                 ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCt                                                                                                                            <  1:1008292/110‑1 (MQ=255)
                                  tttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                              >  2:487052/1‑139 (MQ=255)
                                   ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                             >  2:791249/1‑139 (MQ=255)
                                   ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCCGgcatgcat                                                                                             >  2:650430/1‑139 (MQ=255)
                                     gCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCatat                                                                                           <  2:205342/139‑1 (MQ=255)
                                       cAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                         <  1:791249/139‑1 (MQ=255)
                                                   tGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                             >  1:417437/1‑139 (MQ=255)
                                                      aGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTc                                                                          >  1:426805/1‑139 (MQ=255)
                                                        cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                        >  2:373347/1‑139 (MQ=255)
                                                          aaaCTGCATAAATCCGTTGCGCAAAGTGAACGAACAACCAGCGGCGAATAGCAGGTAAATCAATACCGATCC‑CC‑AAGAGGACGGAGCTAATGACGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGa                                                                      <  1:775593/139‑1 (MQ=255)
                                                          aaaCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGa                                                                      <  1:554090/139‑1 (MQ=255)
                                                            aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                         <  1:279085/118‑1 (MQ=255)
                                                            aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                         >  2:279085/1‑118 (MQ=255)
                                                                cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                >  1:177023/1‑139 (MQ=255)
                                                                  caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:426805/139‑1 (MQ=255)
                                                                  caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:472319/139‑1 (MQ=255)
                                                                     aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGc                                                                                                                             <  2:149154/73‑1 (MQ=255)
                                                                     aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGc                                                                                                                             >  1:149154/1‑73 (MQ=255)
                                                                       ccGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtgtg                                                         <  2:351324/139‑1 (MQ=255)
                                                                         ggTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATTCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTa                                                       <  1:373347/139‑1 (MQ=255)
                                                                         ggTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTa                                                       <  2:40529/139‑1 (MQ=255)
                                                                             cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          >  1:747794/1‑100 (MQ=255)
                                                                             cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          <  2:747794/100‑1 (MQ=255)
                                                                             cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACg                                                   <  1:433204/139‑1 (MQ=255)
                                                                             cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACg                                                   >  2:646297/1‑139 (MQ=255)
                                                                                    tGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCaa                                            <  2:177023/139‑1 (MQ=255)
                                                                                     gAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  2:965079/1‑120 (MQ=255)
                                                                                     gAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  1:965079/120‑1 (MQ=255)
                                                                                     gAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAGACGTTACCAAt                                           <  1:482690/139‑1 (MQ=255)
                                                                                     gAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                           <  1:457908/139‑1 (MQ=255)
                                                                                       aCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  1:387832/118‑1 (MQ=255)
                                                                                       aCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  2:387832/1‑118 (MQ=255)
                                                                                        ccAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGt                                        >  1:932227/1‑139 (MQ=255)
                                                                                          aacaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGttt                                      <  1:646297/139‑1 (MQ=255)
                                                                                           acaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTa                                     >  1:159027/1‑139 (MQ=255)
                                                                                            caacCTGCGCCGAAGAGCAAGTAAATCATTACCGATCC‑CCAAA‑AGGAAGCTGTTAATGAAGGAGAAAAAATGTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTGCATGTGTAAACGTTACCAAATGTTTaa                                    >  1:357698/1‑139 (MQ=255)
                                                                                              acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                  <  1:400922/139‑1 (MQ=255)
                                                                                              acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                  <  2:932227/139‑1 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACg                                                                                                                              <  2:6112/42‑1 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACg                                                                                                                              >  1:6112/1‑42 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAATGAGAAAAAATCTgg                                                                                                    <  1:221352/68‑1 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    <  1:540635/68‑1 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    >  2:221352/1‑68 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                    >  2:540635/1‑68 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtgtg                                                         <  2:369098/111‑1 (MQ=255)
                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCtgtg                                                         >  1:369098/1‑111 (MQ=255)
                                                                                                     ccGAAGAGCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtatat                           <  1:487052/139‑1 (MQ=255)
                                                                                                          gagCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtata                            <  2:20638/133‑1 (MQ=255)
                                                                                                          gagCAGGTAAATCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGtata                            >  1:20638/1‑133 (MQ=255)
                                                                                                                   aaTCATTACCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTAc             <  2:417437/139‑1 (MQ=255)
                                                                                                                          aCCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:630900/83‑1 (MQ=255)
                                                                                                                          aCCGATCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  1:630900/1‑83 (MQ=255)
                                                                                                                             atccc‑ca‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                             <  2:734247/61‑1 (MQ=255)
                                                                                                                             atccc‑ca‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc  >  2:485791/4‑139 (MQ=255)
                                                                                                                              aTCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCg                                                                             >  1:734247/1‑64 (MQ=255)
                                                                                                                              aTCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                   >  1:222903/1‑106 (MQ=255)
                                                                                                                              aTCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                   <  2:222903/106‑1 (MQ=255)
                                                                                                                               tCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                           >  2:181813/1‑97 (MQ=255)
                                                                                                                               tCC‑CCAAA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                           <  1:181813/97‑1 (MQ=255)
                                                                                                                               ccc‑ca‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                <  1:860145/74‑1 (MQ=255)
                                                                                                                               ccc‑ca‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                >  2:860145/2‑75 (MQ=255)
                                                                                                                               ccc‑ca‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTAc             <  1:876431/125‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑AA‑AGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTC  >  NC_000913/7791‑8054

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑