Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 145,794 T→C 7.2% V343A (GTG→GCG)  gltP → glutamate/aspartate‑proton DAACS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947145,7940TC 7.2% 198.8 / 4.7 69V343A (GTG→GCG) gltPglutamate/aspartate‑proton DAACS transporter
Reads supporting (aligned to +/- strand):  ref base T (35/29);  new base C (3/2);  total (38/31)
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.64e-01

CGAAGTCGATCTGCTCGTTCGTGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGACATGGCCCG  >  NC_002947/145653‑145941
                                                                                                                                             |                                                                                                                                                   
cGAAGTCGATCTGCTCGTTCGTGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGAcc                                                                                                                                              >  2:42492/1‑149 (MQ=255)
    gTCGATCTGCTTGTTCGTGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGa                                                                                                                                          <  2:335051/149‑1 (MQ=255)
             cTCGTTCGTGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCAc                                                                                                                                   <  5:241759/147‑1 (MQ=255)
              tcgttcgtGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCAAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACtt                                                                                                                                 >  3:102970/1‑148 (MQ=255)
                gttcgtGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGCCCCAGTACCAGAGCATCGCGGCGCTCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCAGACCCTGATGGTCACTTcc                                                                                                                               >  3:250811/1‑148 (MQ=255)
                    gtGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAg                                                                                                                           <  7:274087/148‑1 (MQ=255)
                      ggTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTa                                                                                                                        <  2:76951/149‑1 (MQ=255)
                           cgaccgGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATcgcg                                                                                                                   >  6:38337/1‑149 (MQ=255)
                            gacccGCTATTCGTTCCACCTAGACGGTTCGACCCAGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGAAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCgg                                                                                                                  <  4:42729/149‑1 (MQ=255)
                              ccgGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGgt                                                                                                                >  1:112515/1‑149 (MQ=255)
                               cgGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGgtg                                                                                                               >  3:201535/1‑149 (MQ=255)
                                ggCTACTCGTTCAACCTAGACGGTTCGACACTCTACCATAGCATCGCTGCGATCTTCATTGACCAGTTCTACTGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTTACCCTGATGGTCACTTCCAAAGGTATCTCGGgtgt                                                                                                              <  4:223942/149‑1 (MQ=255)
                                 gCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGtgtg                                                                                                             >  5:164206/1‑149 (MQ=255)
                                  cTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTc                                                                                                                                <  6:301652/129‑1 (MQ=255)
                                  cTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTc                                                                                                                                >  5:301652/1‑129 (MQ=255)
                                     ttCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCgg                                                                                                         <  2:181145/149‑1 (MQ=255)
                                         ttCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCTTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGa                                                                                                                                          <  1:269409/112‑1 (MQ=255)
                                         ttCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCAGTCGTGGAGCCAGCAGCTGCTGATGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCg                                                                                                      >  1:135358/1‑148 (MQ=255)
                                             aCCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTAtcg                                                                                                 >  5:274143/1‑149 (MQ=255)
                                              ccTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTAtcgt                                                                                                >  6:71690/1‑149 (MQ=255)
                                                   aCGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTg                                                                                           <  5:71690/149‑1 (MQ=255)
                                                     ggTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTgg                                                                                          >  6:11932/1‑148 (MQ=255)
                                                          gACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTgg                                                                                          >  7:163426/1‑143 (MQ=255)
                                                          gACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTgg                                                                                          <  8:163426/143‑1 (MQ=255)
                                                          gACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGAGGGTGTGCCGGGCGTAt                                                                                                   >  2:93400/1‑134 (MQ=255)
                                                          gACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTACTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTActgc                                                                                    >  4:238917/1‑149 (MQ=255)
                                                            cccTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGAGTGCCGGGCGTATCGTTCGTGGTActgctg                                                                                  >  4:132679/1‑149 (MQ=255)
                                                             ccTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTgg                                                                                 >  6:276137/1‑149 (MQ=255)
                                                               tGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTActgctg                                                                                  <  8:183341/146‑1 (MQ=255)
                                                               tGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTActgctg                                                                                  >  7:183341/1‑146 (MQ=255)
                                                                    cAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCtt                                                                          >  4:323948/1‑149 (MQ=255)
                                                                     agagCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGCGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCggg                                                                                                        <  5:231074/118‑1 (MQ=255)
                                                                     agagCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGCGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCggg                                                                                                        >  4:265991/1‑118 (MQ=255)
                                                                     agagCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGCGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCggg                                                                                                        >  6:231074/1‑118 (MQ=255)
                                                                     agagCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGCGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCggg                                                                                                        <  3:265991/118‑1 (MQ=255)
                                                                     agagCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGCGATGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCggg                                                                                                        >  4:261429/1‑118 (MQ=255)
                                                                         cATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGAAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCTGGCGTATCGTTCGTGGTACTGCTGGCCACCTTggg                                                                       >  8:311216/1‑147 (MQ=255)
                                                                           tCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTg                                                                         <  3:148188/143‑1 (MQ=255)
                                                                           tCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTg                                                                         >  4:148188/1‑143 (MQ=255)
                                                                               ggcgATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTg                                                                                                                                                    <  8:298722/64‑1 (MQ=255)
                                                                               ggcgATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTg                                                                                                                                                    >  7:298722/1‑64 (MQ=255)
                                                                               ggcgATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTgg                                                               <  3:132679/149‑1 (MQ=255)
                                                                                     cTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTc                                                         >  3:142075/1‑149 (MQ=255)
                                                                                     cTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTc                                                         >  1:52215/1‑149 (MQ=255)
                                                                                        cATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCGAGCAGCTGCTGCTGGTGCTGACCCTGATgg                                                                                                                                       <  8:88871/68‑1 (MQ=255)
                                                                                        cATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATgg                                                                                                                                       >  7:88871/1‑68 (MQ=255)
                                                                                         aTTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGc                                                      <  6:274143/148‑1 (MQ=255)
                                                                                           tGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGTCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGc                                                             <  1:248816/139‑1 (MQ=255)
                                                                                           tGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGc                                                             >  2:248816/1‑139 (MQ=255)
                                                                                                aGTTGTACGGTATCGACCTATCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAggc                                              >  8:53573/1‑149 (MQ=255)
                                                                                                  ttGTACGGTATCGACCTGTCGTGGAGCCTGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAggcct                                            >  1:192770/1‑149 (MQ=255)
                                                                                                  ttGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAggcct                                            <  2:112515/149‑1 (MQ=255)
                                                                                                     tACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCTCTGGAAggcctggc                                         <  4:142075/149‑1 (MQ=255)
                                                                                                     tACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAggcctggc                                         <  2:52215/149‑1 (MQ=255)
                                                                                                          tATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTc                                     >  6:126058/1‑148 (MQ=255)
                                                                                                            tCGACCTGTCGTGGAGCCAGCAGCTACTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATc                                  <  3:238917/149‑1 (MQ=255)
                                                                                                                ccTGTCGTGGAGCCAGGAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCgg                              <  5:276137/149‑1 (MQ=255)
                                                                                                                  tGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGc                            >  1:247029/1‑149 (MQ=255)
                                                                                                                    tCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGt                          >  8:64339/1‑149 (MQ=255)
                                                                                                                      gTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCATCGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTgg                        <  4:102970/149‑1 (MQ=255)
                                                                                                                      gTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTg                         >  8:292423/1‑148 (MQ=255)
                                                                                                                         gAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGaa                                                 <  3:120953/121‑1 (MQ=255)
                                                                                                                         gAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGaa                                                 >  4:120953/1‑121 (MQ=255)
                                                                                                                           accAGCAGCAGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCg                    >  6:170190/2‑148 (MQ=255)
                                                                                                                                      ctgctgGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGACa         <  3:323948/148‑1 (MQ=255)
                                                                                                                                         ctgGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGCCATgg      <  7:311216/148‑1 (MQ=255)
                                                                                                                                         ctgGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGACATGGc     <  2:192770/149‑1 (MQ=255)
                                                                                                                                            gTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGACATGGCCCg  <  5:11932/149‑1 (MQ=255)
                                                                                                                                            gTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGACATGGCCCg  >  7:119155/1‑149 (MQ=255)
                                                                                                                                             tGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAggcctgg                                          <  3:143068/108‑1 (MQ=255)
                                                                                                                                             tGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAggcctgg                                          >  4:143068/1‑108 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CGAAGTCGATCTGCTCGTTCGTGGTACCGACCGGCTATTCGTTCAACCTAGACGGTTCGACCCTGTACCAGAGCATCGCGGCGATCTTCATTGCCCAGTTGTACGGTATCGACCTGTCGTGGAGCCAGCAGCTGCTGCTGGTGCTGACCCTGATGGTCACTTCCAAAGGTATCGCGGGTGTGCCGGGCGTATCGTTCGTGGTACTGCTGGCCACCTTGGGCAGCGTGGGCATTCCGCTGGAAGGCCTGGCCTTCATCGCGGGCGTGGACCGCATCATGGACATGGCCCG  >  NC_002947/145653‑145941

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