Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,148,110 C→T 6.3% S97S (AGC→AGT arcD‑II → arginine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,148,1100CT 6.3% 261.8 / 4.3 79S97S (AGC→AGTarcD‑IIarginine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base T (3/2);  total (41/38)
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.92e-01

CATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCGGTACA  >  NC_002947/1147975‑1148238
                                                                                                                                       |                                                                                                                                
cATGGTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAgc                                                                                                                                  <  2:281499/136‑1 (MQ=255)
       gCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCctggc                                                                                                                           >  2:137447/1‑136 (MQ=255)
             gCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTAGATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCa                                                                                                                     >  1:334025/1‑136 (MQ=255)
             gCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCa                                                                                                                     >  2:509980/1‑136 (MQ=255)
                 tCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGt                                                                                                                 >  1:380844/1‑136 (MQ=255)
                     gtTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGACGGc                                                                                                             >  1:46605/1‑136 (MQ=255)
                            aCCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCg                                                                                                      >  2:95041/1‑136 (MQ=255)
                                           aaGCCGGATCTGGTCTCGGGGGTATATGCCTTCGCCTAGGCCGGCTTTGGCGAGTACTTGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  2:457011/136‑1 (MQ=21)
                                           aaGCCGGATCTGGACTCGGGGGTATATGCCTACGCCAAGGCGGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTTGTGCTGCTGTTCAGCa                                                                                       <  1:397515/136‑1 (MQ=255)
                                           aaGCCGGATCTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:347638/136‑1 (MQ=255)
                                           aaGCCGGAACTGGTCTCGGGGGTATATGCCTACGCCTAGGGCGGGTTTGGCGAGTTCAAGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  2:511093/136‑1 (MQ=21)
                                           aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:103185/136‑1 (MQ=255)
                                           aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:24141/136‑1 (MQ=255)
                                           aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:95041/136‑1 (MQ=255)
                                           aaGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCAGCCTGGGGCTACTGGATCAGTGCCTGGCTGGGCAACGTCGGCTACTTCGTGCTGCTGTTCAGCa                                                                                       <  1:155145/136‑1 (MQ=255)
                                                    ctggcatCGGGGGGATATGCCTACGCCAATGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTTGATCATCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCt                                                                              <  2:441120/130‑1 (MQ=18)
                                                          tCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTTGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACt                                                                        >  1:474566/1‑136 (MQ=255)
                                                          ggtggggTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACt                                                                        >  1:129479/4‑136 (MQ=255)
                                                           cGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGGCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACtt                                                                       <  2:46605/136‑1 (MQ=255)
                                                                    atGCCTACGCCAAGGCCGGCTTTGGCGAGTTCATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTACTTCGTGCTGCTGTTCAGCACCCTGGGCTACTTCTTCCCgatctt                                                              >  2:363432/1‑131 (MQ=11)
                                                                    atGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGTGCCTGGCTGGGCAACGTCGGCTACTTCGTGCTGCTGTTCAGCACCCTGGGCTACTTCTTCCCgatctt                                                              >  1:426761/1‑131 (MQ=11)
                                                                     tGCCTACGCCTAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCAGCCTGGGGCTACTGGATCAGTGCCTGGCTGGGCaa                                                                                                                    <  2:435253/81‑1 (MQ=17)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  2:426125/110‑1 (MQ=255)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       >  2:25514/1‑110 (MQ=255)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       >  2:420808/1‑110 (MQ=255)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:25514/110‑1 (MQ=255)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:420808/110‑1 (MQ=255)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       >  1:426125/1‑110 (MQ=255)
                                                                     tGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCAGCCTGGGGCTACTGGATCAGTGCCTGGCTGGGCaa                                                                                                                    >  1:435253/1‑81 (MQ=255)
                                                                      gCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       >  1:425189/1‑109 (MQ=255)
                                                                      gCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  2:425189/109‑1 (MQ=255)
                                                                           cGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGaa                                                       >  1:518168/1‑136 (MQ=255)
                                                                              cAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGc                                                    >  1:266881/1‑136 (MQ=255)
                                                                              cAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGc                                                    >  1:515044/1‑136 (MQ=255)
                                                                              cAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGc                                                    <  1:73883/136‑1 (MQ=255)
                                                                               aaGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCTGCGAAGGCa                                                   <  1:404826/136‑1 (MQ=255)
                                                                               aaGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCa                                                   <  1:215426/136‑1 (MQ=255)
                                                                                   ccGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTTGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAacac                                               <  2:474566/136‑1 (MQ=255)
                                                                                     ggCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTAtt                                                                                                         <  2:386018/76‑1 (MQ=255)
                                                                                     ggCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTAtt                                                                                                         >  1:386018/1‑76 (MQ=255)
                                                                                     ggCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAacacac                                             <  1:509980/136‑1 (MQ=255)
                                                                                         ttGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAgg                                                     <  1:546790/124‑1 (MQ=255)
                                                                                         ttGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAgg                                                     >  2:546790/1‑124 (MQ=255)
                                                                                         ttGGCGAGTACATGGGCTTTTACTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACccaa                                         >  2:208409/1‑135 (MQ=255)
                                                                                               aGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGtt                                                                                            <  1:60590/79‑1 (MQ=255)
                                                                                               aGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGtt                                                                                            >  2:60590/1‑79 (MQ=255)
                                                                                               aGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACccatcgccat                                   <  2:495145/136‑1 (MQ=255)
                                                                                               aGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACccatcgccat                                   <  2:505981/136‑1 (MQ=255)
                                                                                               aGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACccatcgccat                                   <  2:206544/136‑1 (MQ=255)
                                                                                                    aTGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGtt                              >  2:445594/1‑136 (MQ=255)
                                                                                                       ggCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTtgtg                           <  2:515044/136‑1 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCaa                                                                                                                    >  1:443019/1‑45 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCaa                                                                                                                    >  1:437347/1‑45 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCaa                                                                                                                    <  2:437347/45‑1 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCaa                                                                                                                    <  2:443019/45‑1 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTtg                             <  2:365218/132‑1 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTtg                             >  1:365218/1‑132 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCa                         >  1:29240/1‑136 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCa                         >  1:133601/1‑136 (MQ=255)
                                                                                                         cTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCa                         >  1:327663/1‑136 (MQ=255)
                                                                                                             tCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCgc                     >  1:106756/1‑136 (MQ=255)
                                                                                                             tCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCgc                     >  2:408469/1‑136 (MQ=255)
                                                                                                              ccTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCgct                    >  1:376811/1‑136 (MQ=255)
                                                                                                                      tGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCtgt            <  1:236277/136‑1 (MQ=255)
                                                                                                                       ggggCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCtgtg           >  2:488950/1‑136 (MQ=255)
                                                                                                                       ggggCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCtgtg           <  2:129479/136‑1 (MQ=255)
                                                                                                                       ggggCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCtgtg           <  2:29240/136‑1 (MQ=255)
                                                                                                                          gCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTc        <  1:137447/136‑1 (MQ=255)
                                                                                                                           cTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGc                                                                                        <  1:374039/55‑1 (MQ=255)
                                                                                                                           cTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGc                                                                                        >  2:374039/1‑55 (MQ=255)
                                                                                                                           cTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       >  2:547516/1‑56 (MQ=255)
                                                                                                                           cTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCa                                                                                       <  1:547516/56‑1 (MQ=255)
                                                                                                                            tACTGGATCAGTGCCTGGCTGGGCAACGTCGGCTACTTCGTGCTGCTGTTCAGCACCCTGGGCTACTTCTTCCCGATCTTCGGTGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCgg      >  1:349116/1‑136 (MQ=14)
                                                                                                                            tACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTtg                             >  1:127593/1‑113 (MQ=255)
                                                                                                                            tACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTtg                             <  2:127593/113‑1 (MQ=255)
                                                                                                                            tACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCgg      >  2:483201/1‑136 (MQ=255)
                                                                                                                            tACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCgg      >  1:132123/1‑136 (MQ=255)
                                                                                                                            tACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCgg      >  1:473051/1‑136 (MQ=255)
                                                                                                                             aCTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCGGt     <  1:408469/136‑1 (MQ=255)
                                                                                                                                ggATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCGGTACa  >  1:453311/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                
CATGTTGGCCCTGGCCTTCGTGTTCCAGACCCTGGCCAACCGCAAGCCGGAACTGGACTCGGGGGTATATGCCTACGCCAAGGCCGGCTTTGGCGAGTACATGGGCTTTTCCTCGGCCTGGGGCTACTGGATCAGCGCCTGGCTGGGCAACGTCGGCTATTTCGTGCTGCTGTTCAGCACCCTGGGCTTCTACTTCCCGGTGTTCGGCGAAGGCAACACACCCATCGCCATCGGTTGTGCATCGCTGCTGCTGTGGTCGGTACA  >  NC_002947/1147975‑1148238

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