Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,129,437 A→T 12.5% A40A (GCT→GCA gcvP‑I ← glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,129,4370AT12.5% 186.0 / 19.3 72A40A (GCT→GCAgcvP‑Iglycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base T (4/5);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

TGGGTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTGGC  >  NC_002947/1129288‑1129585
                                                                                                                                                     |                                                                                                                                                    
tgGGTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCa                                                                                                                                                      >  2:856245/1‑150 (MQ=255)
  ggTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGc                                                                                                                                                    <  2:243000/150‑1 (MQ=255)
  ggTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGc                                                                                                                                                    >  2:635854/1‑150 (MQ=255)
     gttgtAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGt                                                                                                                                                 >  2:1395901/1‑150 (MQ=255)
        gtagtaGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCAt                                                                                                                                              <  1:348129/150‑1 (MQ=255)
          agtagCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCTGCGGTCAtcc                                                                                                                                            <  1:1176646/150‑3 (MQ=255)
           gtagCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGc                                                                                                                                           <  2:107234/150‑1 (MQ=255)
               cccTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTcc                                                                                                                                       >  1:240489/1‑150 (MQ=255)
                ccTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCa                                                                                                                                      >  1:176148/1‑150 (MQ=255)
                          gTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGaa                                                                                                                            >  2:878988/1‑150 (MQ=255)
                           tAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAg                                                                                                                           >  1:1121325/1‑150 (MQ=255)
                               ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCca                                                                                                                       >  1:1322913/1‑150 (MQ=255)
                               ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCca                                                                                                                       <  2:1414112/150‑1 (MQ=255)
                               ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCca                                                                                                                       >  1:1057054/1‑150 (MQ=255)
                               ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCca                                                                                                                       >  1:1078642/1‑150 (MQ=255)
                                   gAACAGCTGGTTCTTGCCAGCAATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCTGCGGTCAtcc                                                                                                                                            >  1:209895/1‑123 (MQ=255)
                                   gAACAGCTGGTTCTTGCCAGCAATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCTGCGGTCAtcc                                                                                                                                            <  2:209895/125‑3 (MQ=255)
                                        gCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGa                                                                                                              >  1:780092/1‑150 (MQ=255)
                                         cTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAg                                                                                                             <  1:646217/150‑1 (MQ=255)
                                            gTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCAt                                                                                                          >  2:1270744/1‑150 (MQ=255)
                                             ttCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATg                                                                                                         <  2:1078642/150‑1 (MQ=255)
                                                ttGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGcc                                                                                                      >  2:564491/1‑150 (MQ=255)
                                                    cAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGc                                                                                                       >  2:524076/1‑145 (MQ=255)
                                                    cAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTATCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGc                                                                                                       <  1:524076/145‑1 (MQ=255)
                                                              cTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAgcggc                                                                                        >  2:869032/1‑150 (MQ=255)
                                                                          gagcgCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTAt                                                                            >  2:563247/1‑150 (MQ=255)
                                                                              gcgTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAg                                                                                                                           <  2:630341/99‑1 (MQ=255)
                                                                              gcgTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAg                                                                                                                           >  1:630341/1‑99 (MQ=255)
                                                                                    gCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCgg                                                                                  <  1:1146752/134‑1 (MQ=255)
                                                                                    gCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCgg                                                                                  >  2:1146752/1‑134 (MQ=255)
                                                                                      cTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGAt                                                                <  1:1395901/150‑1 (MQ=255)
                                                                                      cTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGAt                                                                >  2:262088/1‑150 (MQ=255)
                                                                                          ctctGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGCAGGGCCAGTCGCGATTCCTTATTTCTTTCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAAt                                                            <  1:1236956/150‑1 (MQ=255)
                                                                                             tGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCg                                                         <  1:262088/150‑1 (MQ=255)
                                                                                               gCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGtt                                                       >  2:1039668/1‑150 (MQ=255)
                                                                                                 cgtcgtGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGc                                                                   >  2:197599/1‑136 (MQ=255)
                                                                                                 cgtcgtGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGc                                                                   <  1:197599/136‑1 (MQ=255)
                                                                                                 cgtcgtGCGAGCCAAGCACGCTGGTACCTGTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTgg                                                     >  1:399576/1‑150 (MQ=255)
                                                                                                       gCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGc                                               <  1:869032/150‑1 (MQ=255)
                                                                                                       gCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGc                                               >  1:451251/1‑150 (MQ=255)
                                                                                                        cGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGcc                                              >  1:1145637/1‑150 (MQ=255)
                                                                                                        cGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGcc                                              >  1:1145568/1‑150 (MQ=255)
                                                                                                          aGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAGGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGa                                            <  1:878988/150‑1 (MQ=255)
                                                                                                          aGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGa                                            <  2:780092/150‑1 (MQ=255)
                                                                                                           gCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATg                                                               <  1:396830/130‑1 (MQ=255)
                                                                                                           gCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATg                                                               >  2:396830/1‑130 (MQ=255)
                                                                                                           gCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAg                                           <  2:1057054/150‑1 (MQ=255)
                                                                                                            ccAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTACAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAgg                                          >  1:928246/1‑150 (MQ=255)
                                                                                                             cAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGt                                         <  1:1039668/150‑1 (MQ=255)
                                                                                                                      tGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGa                                            <  2:814066/138‑1 (MQ=255)
                                                                                                                      tGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGa                                            >  1:814066/1‑138 (MQ=255)
                                                                                                                        gTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATg                              <  1:81926/150‑1 (MQ=255)
                                                                                                                           ccTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCg                           >  1:2882/1‑150 (MQ=255)
                                                                                                                            cTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCgg                          <  1:1146031/150‑1 (MQ=255)
                                                                                                                            cTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCgg                          <  2:2882/150‑1 (MQ=255)
                                                                                                                             tttGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGGAATAGTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCggg                         >  1:74188/1‑150 (MQ=255)
                                                                                                                                 aTGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGACATCGTAGTGGGTGCCGAGGTTGATGGTCATGGCGGAATTc                     >  2:690334/1‑150 (MQ=255)
                                                                                                                                    cTGTCGGGGACGACCGCAGCGGTCATGGCTTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTTATGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATCCCTc                  <  2:117906/150‑1 (MQ=255)
                                                                                                                                     tGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCa                 <  1:856245/150‑1 (MQ=255)
                                                                                                                                        gggggATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAGGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGc              <  1:81034/149‑1 (MQ=255)
                                                                                                                                          gccaTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCTTTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACgg                                                                     <  1:748520/90‑1 (MQ=255)
                                                                                                                                          gccaTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACgg                                                                     >  2:748520/4‑93 (MQ=255)
                                                                                                                                           ccaTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGc                                                                                                       <  2:1033559/56‑1 (MQ=39)
                                                                                                                                           ccaTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGc                                                                                                       >  1:1033559/3‑58 (MQ=39)
                                                                                                                                            gATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCa                                                                                                           >  1:692251/1‑53 (MQ=255)
                                                                                                                                            gATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCa                                                                                                           <  2:692251/53‑1 (MQ=255)
                                                                                                                                            caTGACCGCTGCGGTCATGGCGTCCCGGGCCACGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcgtcg          >  1:37260/2‑150 (MQ=255)
                                                                                                                                              tGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCCACACGCCGCTCTTGGTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcgtcgtc        <  2:37260/150‑1 (MQ=255)
                                                                                                                                                 ccGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAg                <  2:266585/139‑1 (MQ=255)
                                                                                                                                                 ccGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAg                >  1:266585/1‑139 (MQ=255)
                                                                                                                                                   gCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTgg   <  2:465708/150‑1 (MQ=255)
                                                                                                                                                   gCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTgg   >  1:465708/1‑150 (MQ=255)
                                                                                                                                                    cAGCGGACCTGCCGACCAGAGAGTCGAAGCCCAGCGCACTGAGCATGGCCTGTGCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCGGAGGGTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTGGc  >  2:1303224/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGGTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTGGC  >  NC_002947/1129288‑1129585

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