Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,003,4800AG21.1% 145.3 / 18.7 71G289G (GGA→GGGxdhBputative xanthine dehydrogenase FAD‑binding subunit XdhB
Reads supporting (aligned to +/- strand):  ref base A (23/33);  new base G (15/0);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.80e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCAAACAA  >  NC_000913/3003350‑3003614
                                                                                                                                  |                                                                                                                                      
cTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCaa                                                                                                                               <  1:73085/140‑1 (MQ=255)
  aTCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATg                                                                                                                             >  1:211705/1‑140 (MQ=255)
  aTCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATg                                                                                                                             >  1:16920/1‑140 (MQ=255)
  aTCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATGCaaaa                                                                                                                             >  2:1110640/1‑138 (MQ=255)
     aGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAAt                                                                                                                          <  1:774823/140‑1 (MQ=255)
     aGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAAt                                                                                                                          <  2:1031987/140‑1 (MQ=255)
      gCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTTCAAAGaaaa                                                                                                                         >  2:1224549/1‑138 (MQ=255)
               tCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGaaa                                                                                                                >  1:1063423/1‑140 (MQ=255)
               tCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGaaa                                                                                                                <  1:1090627/140‑1 (MQ=255)
                   gCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAAt                                                                                                            <  1:970987/140‑1 (MQ=255)
                    cAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGAAACAAtt                                                                                                           >  1:706145/1‑140 (MQ=255)
                     aaGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTa                                                                                                          <  1:1192298/140‑1 (MQ=255)
                      aGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGAAAAAAATAc                                                                                                         >  1:1097919/1‑140 (MQ=255)
                      aGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAATGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTAc                                                                                                         >  1:1265271/1‑140 (MQ=255)
                       gATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTAcc                                                                                                        <  1:389207/140‑1 (MQ=255)
                        aTGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCa                                                                                                       <  2:100262/140‑1 (MQ=255)
                        aTGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCa                                                                                                       <  2:1265271/140‑1 (MQ=255)
                         tGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCTCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGaa                                                                                                                 <  2:1190574/129‑1 (MQ=255)
                         tGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGaa                                                                                                                 >  1:1190574/1‑129 (MQ=255)
                         tGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATGCAAAGAAATACAGCGAAAAAAATACCac                                                                                                      >  2:1147855/1‑139 (MQ=255)
                         tGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCAt                                                                                                      >  1:1062001/1‑140 (MQ=255)
                         tGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGGGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGAAACAAATACCAt                                                                                                      >  1:720461/1‑140 (MQ=255)
                          gTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATc                                                                                                     <  2:1165682/140‑1 (MQ=255)
                            cGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGGGTTTCGTCTGCATCTCATCCGGACGATGGCCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAAATGGAATTAAATACAGGGGAACAATTACCATCGa                                                                                                   >  1:540770/1‑140 (MQ=255)
                              ccccGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAAt                                                                                                 <  2:162563/140‑1 (MQ=255)
                              ccccGCGTTCTTCATGGCGGGCCAGTAAAGAGGTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATGCCATGAATTACAGCGCAACAATTACCATCGaaa                                                                                                 >  2:176632/1‑139 (MQ=255)
                               cccGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATg                                                                                                <  2:1460367/140‑1 (MQ=255)
                                ccGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGc                                                                                               <  1:1464097/140‑1 (MQ=255)
                                  gcgTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAACGAAACAATTAACAACGAATGCAc                                                                                             >  1:1005865/1‑140 (MQ=255)
                                        ttcATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGAAACAATTACCATCg                                                                                                    >  1:1256027/1‑127 (MQ=255)
                                        ttcATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCg                                                                                                    <  2:1256027/127‑1 (MQ=255)
                                              gCGGGCCAGTAAAGAGGTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATGCAATGAATGACAACGAAACAAATACCATCGAATGCCCCCTTTACGGGAt                                                                                 >  2:1082180/1‑140 (MQ=255)
                                                  gCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATGCAATGAATCACAGCGAAACAAATACCATCGAAAACACCAATAAAGGGATGGCt                                                                             >  2:1027542/1‑140 (MQ=255)
                                                  gCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCt                                                                             >  1:467899/1‑140 (MQ=255)
                                                          agagTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCt                                                                     <  2:1062001/140‑1 (MQ=255)
                                                           gagTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCtt                                                                    >  1:1129861/1‑140 (MQ=255)
                                                            agTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTc                                                                   >  2:294782/1‑140 (MQ=255)
                                                            agTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTc                                                                   >  1:1293226/1‑140 (MQ=255)
                                                              tttCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAc                                                                 >  2:1142391/1‑140 (MQ=255)
                                                              tttCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAc                                                                 <  1:207338/140‑1 (MQ=255)
                                                                tCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAcgc                                                               >  2:761088/1‑140 (MQ=255)
                                                                tCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTAAcgc                                                               <  2:1129861/140‑1 (MQ=255)
                                                                 cgccTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGCGGGAAAATTGCCACGAATCACAGGGAAAAAATTACCATCGAAAGCACCATTAACGGGATGCCTTTTCAGGTTCAcacc                                                              >  2:393898/4‑140 (MQ=255)
                                                                 cGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGGACCATTAACGGGATGCCTTTTCAGCTTCAcgcc                                                              >  1:1303970/1‑140 (MQ=255)
                                                                        aTCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAg                                                       <  2:706145/140‑1 (MQ=255)
                                                                         tctcATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAcg                                                                >  2:344866/1‑130 (MQ=255)
                                                                         tctcATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAcg                                                                <  1:344866/130‑1 (MQ=255)
                                                                             aTCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGACAcg                                                  >  1:1311013/1‑140 (MQ=255)
                                                                              tCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAATTGCAATGAATCAACGCGAAACAATTACCATCGAATGCACCATTAAAGGGATGTCTTTTCAGGTTGACGCCCCACCAGGGAccc                                                 >  2:372660/1‑138 (MQ=255)
                                                                              tCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCATCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAAGCAcgc                                                 >  2:860058/1‑140 (MQ=255)
                                                                              tCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCAcgc                                                 >  2:481427/1‑140 (MQ=255)
                                                                                 aGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCAcgccgc                                              >  1:414081/1‑140 (MQ=255)
                                                                                  gACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGct                                             <  2:1063423/140‑1 (MQ=255)
                                                                                     gATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCtctc                                          >  2:463598/1‑140 (MQ=255)
                                                                                      aTGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCg                                         <  2:36698/140‑1 (MQ=255)
                                                                                       tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCgg                                        >  2:1335564/1‑140 (MQ=255)
                                                                                       tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCgg                                        >  2:334831/1‑140 (MQ=255)
                                                                                         aCCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGaa                                      >  1:733172/1‑140 (MQ=255)
                                                                                              aaaaGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACt                                 >  1:1314163/1‑140 (MQ=255)
                                                                                               aaaGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTc                                <  1:1082180/140‑1 (MQ=255)
                                                                                               aaaGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTc                                <  2:414081/140‑1 (MQ=255)
                                                                                                 aGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCcg                              <  1:874429/140‑1 (MQ=255)
                                                                                                 aGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCcg                              <  1:463598/140‑1 (MQ=255)
                                                                                                 aGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCcg                              >  1:420143/1‑140 (MQ=255)
                                                                                                            aaGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGAc                   >  2:556/1‑140 (MQ=255)
                                                                                                                cGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCt               <  2:467899/140‑1 (MQ=255)
                                                                                                                cGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCt               <  1:1110640/140‑1 (MQ=255)
                                                                                                                  tCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTaa             <  1:294782/140‑1 (MQ=255)
                                                                                                                  tCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTaa             <  2:211705/140‑1 (MQ=255)
                                                                                                                        gggcggGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATCCACCATAAACGGGATCCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCa       <  2:1263619/139‑1 (MQ=255)
                                                                                                                          gcggGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCaaa     >  2:1176220/1‑140 (MQ=255)
                                                                                                                          gcggGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCaaa     <  1:1147855/140‑1 (MQ=255)
                                                                                                                            ggggggAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCAAACa   <  2:1097919/140‑1 (MQ=255)
                                                                                                                             gggggAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCAAACaa  >  1:57019/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                      
CTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAGTGTCAAACAA  >  NC_000913/3003350‑3003614

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