Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,900,160 C→A 9.3% P175Q (CCG→CAG)  PP_3443 → glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,900,1600CA 9.3% 198.5 / 6.8 64P175Q (CCG→CAG) PP_3443glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base A (3/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

TCGGCAAGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTC  >  NC_002947/3900013‑3900302
                                                                                                                                                   |                                                                                                                                              
cccgCAATAACCTGAATGATTCCGAAACTTAATGCGCCCGCACCTGCGACTACATCGTCGACACTATCAACGCCGTCAAGGACATGGACCGCCTCTCCAGCCGCTTCCACCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGc                                                                                                                                              <  2:551591/147‑1 (MQ=255)
    cAAGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTggg                                                                                                                                          <  1:859778/150‑1 (MQ=255)
        aaCCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTg                                                                                                                                      <  1:939232/150‑1 (MQ=255)
        aaCCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTg                                                                                                                                      <  1:663252/150‑1 (MQ=255)
           cTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCAGCGGGGCGCGGcg                                                                                                                                   <  1:97479/150‑2 (MQ=255)
                 gATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCACGCCCAGCGGGGCGCGGcggatctg                                                                                                                             <  2:927484/150‑8 (MQ=255)
                        aaaaGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGc                                                                                                                      >  2:1055329/1‑150 (MQ=255)
                            ggAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTa                                                                                                                  <  1:895533/150‑1 (MQ=255)
                             gAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTAc                                                                                                                 >  1:809359/1‑150 (MQ=255)
                               aTTCGACCGCACCCCCGACTAGATCGTCAACACCATCAACTCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGTCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACaa                                                                                                               <  2:630743/150‑1 (MQ=255)
                                  cGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTa                                                                                                            <  2:248523/150‑1 (MQ=255)
                                   gACCGCAGCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTAc                                                                                                           <  2:102075/150‑1 (MQ=255)
                                      cGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCg                                                                                                        <  2:997636/150‑1 (MQ=255)
                                               gACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACgag                                                                                               >  1:1129477/1‑150 (MQ=255)
                                                       cgtcgACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTcac                                                                                       <  2:891148/150‑1 (MQ=255)
                                                        gtcgACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTg                                                                                                     >  1:339925/1‑135 (MQ=255)
                                                        gtcgACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTg                                                                                                     <  2:339925/135‑1 (MQ=255)
                                                            acacCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGc                                                                                  >  1:52043/1‑150 (MQ=255)
                                                              acCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTc                                                                                >  2:1258478/1‑150 (MQ=255)
                                                                    aaCGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCg                                                                          >  1:223501/1‑150 (MQ=255)
                                                                     aCGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCgg                                                                         >  1:1068253/1‑150 (MQ=255)
                                                                      cGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGgc                                                                        <  1:1055329/150‑1 (MQ=255)
                                                                        cccTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGgcgc                                                                      >  1:467040/1‑150 (MQ=255)
                                                                         ccTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCt                                                                     >  2:645441/1‑150 (MQ=255)
                                                                           tCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGa                                                                   <  1:1200158/150‑1 (MQ=255)
                                                                               ggACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCAt                                                               >  2:925634/1‑150 (MQ=255)
                                                                               ggACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCAt                                                               <  2:1129477/150‑1 (MQ=255)
                                                                                gACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATg                                                              <  1:925634/150‑1 (MQ=255)
                                                                                  ccTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATggg                                                            <  2:938888/150‑1 (MQ=255)
                                                                                   cTGGACCGCCGCTCCAGCCGCTTCGAGCTGGCGCAAGCCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGc                                                           >  1:1316841/1‑150 (MQ=255)
                                                                                   cTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGc                                                           >  1:1331824/1‑150 (MQ=255)
                                                                                         cgccgcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACg                                                     <  2:1174125/150‑1 (MQ=255)
                                                                                             gcTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGATGATAATGGGCAACAgggggg                                                 >  2:1371058/1‑150 (MQ=255)
                                                                                               tCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGgtggtg                                               >  1:267543/1‑150 (MQ=255)
                                                                                                ccAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGgtgt                                              >  1:463590/1‑150 (MQ=255)
                                                                                                 cAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGgtgt                                              >  1:1225984/1‑149 (MQ=255)
                                                                                                 cAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGgtgt                                              <  2:1225984/149‑1 (MQ=255)
                                                                                                 cAGCCGCTGCGAGCGGGAGCAAGGCACCGTGGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGtt                                             >  2:415554/1‑150 (MQ=255)
                                                                                                   gCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCa                                           >  1:803007/1‑150 (MQ=255)
                                                                                                    ccGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCaa                                          <  2:1331824/150‑1 (MQ=255)
                                                                                                       cTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGcc                                       >  2:376918/1‑150 (MQ=255)
                                                                                                        ttCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACTAATTCCCGCTTAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCg                                      >  2:115046/1‑150 (MQ=255)
                                                                                                         agcggCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATTATATGTAACACGGTGGTGTTCAAGCCgg                                     <  1:903307/146‑1 (MQ=255)
                                                                                                              cTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGc                                    >  1:422734/1‑146 (MQ=255)
                                                                                                              cTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGc                                    <  2:422734/146‑1 (MQ=255)
                                                                                                               tGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACaa                                                                                                               >  1:620348/1‑70 (MQ=255)
                                                                                                               tGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACaa                                                                                                               <  2:620348/70‑1 (MQ=255)
                                                                                                                ggAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGtt                              >  2:607703/1‑150 (MQ=255)
                                                                                                                       ggCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCaa                                          >  1:1073275/1‑131 (MQ=255)
                                                                                                                       ggCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCaa                                          <  2:1073275/131‑1 (MQ=255)
                                                                                                                       ggCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTa                       <  1:1258478/150‑1 (MQ=255)
                                                                                                                              tCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGa                >  1:435635/1‑150 (MQ=255)
                                                                                                                                ggCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATc              <  1:1371058/150‑1 (MQ=255)
                                                                                                                                 gCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATcc             <  1:645441/150‑1 (MQ=255)
                                                                                                                                    aGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGAc          >  2:643121/1‑150 (MQ=255)
                                                                                                                                    aGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGAc          <  2:809359/150‑1 (MQ=255)
                                                                                                                                     acagggccacGCCCAGCGGGGCGCGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCgg                                     >  2:1270035/6‑122 (MQ=255)
                                                                                                                                     acagggccacGCCCAGCGGGGCGCGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCgg                                     <  1:1270035/117‑1 (MQ=255)
                                                                                                                                      aTCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGc                                                                                  <  1:922540/76‑1 (MQ=255)
                                                                                                                                      aTCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTCAACGAGACCTTCACCACGc                                                                                  >  2:922540/1‑76 (MQ=255)
                                                                                                                                       agggccacgcCCAGCGGGGCGCGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCt       >  2:153098/8‑150 (MQ=255)
                                                                                                                                            cgcgCCTCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACACCTTCCCGCTGAACGAGATCTTCACCACGCTCATCCCGGGGCTGATCATGGGCAACACGGTAGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGACCCGACCCTTGCTc  <  2:492246/150‑1 (MQ=255)
                                                                                                                                            cgcgCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTc  <  2:463590/150‑1 (MQ=255)
                                                                                                                                            caggcCCAGCGGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTc  >  1:1423797/4‑150 (MQ=255)
                                                                                                                                                 cccGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATg                                                              >  1:1151690/1‑85 (MQ=255)
                                                                                                                                                 cccGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATg                                                              <  2:1151690/85‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TCGGCAAGAACCTGAAGGATTCGGAAAAGGAATTCGACCGCACCTGCGACTACATCGTCGACACCATCAACGCCCTCAAGGACCTGGACCGCCGCTCCAGCCGCTTCGAGCTGGAGCAAGGCACCCTCGGCCAGATCCGCCGCGCCCCGCTGGGCGTGGCCCTGTGCATGGGGCCCTACAACTACCCGCTGAACGAGACCTTCACCACGCTCATCCCGGCGCTGATCATGGGCAACACGGTGGTGTTCAAGCCGGCCAAGTTCGGCGTACTGCTGATCCGACCCTTGCTC  >  NC_002947/3900013‑3900302

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: