Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,632,992 A→T 5.4% W88R (TGG→AGG)  pgsA ← CDP‑diacylglycerol‑glycerol‑3‑phosphate 3‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,632,9920AT 5.4% 229.4 / 2.5 74W88R (TGG→AGG) pgsACDP‑diacylglycerol‑glycerol‑3‑phosphate 3‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base T (2/2);  total (37/37)
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.27e-01

CCAGCATCTGCGCCGCAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGCCAATCGGTGGCGGCGGCTACG  >  NC_002947/4632851‑4633140
                                                                                                                                                |                                                                                                                                                    
ccAGCATCTGCGGCGCAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTg                                                                                                                                               <  6:26001/149‑1 (MQ=255)
     aTCTGCGCCGCAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTg                                                                                                                                          >  6:75278/1‑149 (MQ=255)
      tCTGCGCCGCAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTgg                                                                                                                                         >  2:260639/1‑149 (MQ=255)
          cgccgcAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTggg                                                                                                                                        <  2:192526/146‑1 (MQ=255)
            ccgcAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGc                                                                                                                                   <  6:71284/149‑1 (MQ=255)
               cAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCAc                                                                                                                                 >  6:227253/1‑148 (MQ=255)
               cAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATACACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTCCAc                                                                                                                                 >  1:198647/1‑148 (MQ=255)
                      ccACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAgcagca                                                                                                                         >  3:202395/1‑149 (MQ=255)
                       cACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTTGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCAc                                                                                                                        >  3:17192/1‑149 (MQ=255)
                         cTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCa                                                                                                                      <  4:278050/149‑1 (MQ=255)
                            aCCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAggg                                                                                                                   >  1:203448/1‑149 (MQ=255)
                             ccGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAgg                                                                                                                    >  2:170181/1‑147 (MQ=255)
                             ccGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAgg                                                                                                                    <  1:170181/147‑1 (MQ=255)
                              cGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGc                                                                                                                  >  1:196735/1‑148 (MQ=255)
                              cGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGc                                                                                                                  >  5:17505/1‑148 (MQ=255)
                                       acacCGCCACATGCGCCGGCGCGCCTAGCTTGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGccacc                                                                                                        <  4:17192/149‑1 (MQ=255)
                                       acacCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTgg                                                                                                                                         <  6:276093/116‑1 (MQ=255)
                                       acacCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTgg                                                                                                                                         >  5:276093/1‑116 (MQ=255)
                                       acacCGCCACAATCGCCCGCGCGGCTATCTCGGCCATCCACTCGCGCCGCGCATCTATACCTACTACCAC‑AAGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGcca                                                                                                          <  1:260639/149‑1 (MQ=255)
                                         acCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTg                                                                                                                                               >  5:81507/1‑108 (MQ=255)
                                         acCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTg                                                                                                                                               <  6:81507/108‑1 (MQ=255)
                                           cGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGccaccgt                                                                                                      >  4:73782/1‑145 (MQ=255)
                                            gCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCCTCAg                                                                                                   <  3:73782/149‑1 (MQ=255)
                                             ccACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGATTCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGc                                                                                                  <  1:105817/149‑1 (MQ=255)
                                               acaTGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCTGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCtt                                                                                                >  6:264396/1‑149 (MQ=255)
                                               acaTGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCtt                                                                                                <  8:92754/149‑1 (MQ=255)
                                               acaTGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCACCGGCAGGGTCAGCCTGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCtt                                                                                                >  4:152542/1‑149 (MQ=255)
                                                  tGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTc                                                                                             <  6:17505/149‑1 (MQ=255)
                                                    cgcCCGCTCGCCTAGCTCGGCGATCCTCTCGCGCAGCGCA‑‑GACACT‑ACTTTCTCGCGGCCGATGATGACAGCCGCAGTCAGGGTCAGCCAGAAGTTGTCGTGGGTCTGCACCATCAGCACGAGGGCCACTGCCACCATCAGCTTGTCg                                                                                            <  2:198647/148‑1 (MQ=255)
                                                      cccGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGcc                                                                                         >  3:17157/1‑149 (MQ=255)
                                                              ccTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTc                                                                                 >  8:206575/1‑149 (MQ=255)
                                                               cTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTc                                                                                 >  1:161876/1‑148 (MQ=255)
                                                               cTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTc                                                                                 >  3:239272/1‑148 (MQ=255)
                                                               cTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTc                                                                                 >  5:289956/1‑148 (MQ=255)
                                                               cTAGCTCGGCCATCCACTCGCGCAGAGCA‑‑GACACC‑ACTATCTCGCGGACGATGATGACCGCAGCCGGCAGGGTCAGACAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGTCCACTGCCACCATCAGCTTGTCGGCCACCGGGTc                                                                                 >  1:166907/1‑148 (MQ=255)
                                                                  gCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCggg                                                                                   >  2:42633/1‑143 (MQ=255)
                                                                  gCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCggg                                                                                   <  1:42633/143‑1 (MQ=255)
                                                                  gCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGa                                                                               <  2:206429/147‑1 (MQ=255)
                                                                     cGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAg                                                                          <  2:196735/149‑1 (MQ=255)
                                                                      ggCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAgg                                                                         >  4:146977/1‑149 (MQ=255)
                                                                      ggCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAg                                                                          <  2:161876/148‑1 (MQ=255)
                                                                      ggCCATCCACTCGAGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCTCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAg                                                                          <  2:166907/148‑1 (MQ=255)
                                                                       gCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGgc                                                                        <  5:227253/149‑1 (MQ=255)
                                                                              aCTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAAc                                                                 >  5:52237/1‑149 (MQ=255)
                                                                                  gcgcagcgca‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCTTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTc                                                                                             <  7:86032/117‑1 (MQ=255)
                                                                                  gcgcagcgca‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTc                                                                                             >  8:86032/1‑117 (MQ=255)
                                                                                  gcgcagcgca‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGc                                                              <  7:90102/148‑1 (MQ=255)
                                                                                   cgcagcgca‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCTGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGt                                                            <  5:264396/149‑1 (MQ=255)
                                                                                   cgcagcgca‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCACCGGCAGGGTCAGCCTGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGt                                                            <  3:152542/149‑1 (MQ=255)
                                                                                    gcagcgca‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCg                                                             >  4:171618/1‑147 (MQ=255)
                                                                                    gcagcgca‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCg                                                             <  3:171618/147‑1 (MQ=255)
                                                                                    gcagcgca‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTg                                                           <  3:107848/149‑1 (MQ=255)
                                                                                         gca‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTctg                                                      <  4:202395/149‑1 (MQ=255)
                                                                                           a‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTctgct                                                    <  2:203448/149‑1 (MQ=255)
                                                                                           a‑‑gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTctgct                                                    <  7:206575/149‑1 (MQ=255)
                                                                                              gACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTctgctg                                                   <  5:127535/149‑1 (MQ=255)
                                                                                               acacC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTctgctg                                                   >  1:103081/1‑148 (MQ=255)
                                                                                                cacc‑acGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCa                                                 <  6:289956/149‑1 (MQ=255)
                                                                                                   c‑acGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGc                                                              >  8:150215/1‑133 (MQ=255)
                                                                                                   c‑acGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGc                                                              <  7:150215/133‑1 (MQ=255)
                                                                                                   c‑acGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTc                                              >  2:276381/1‑149 (MQ=255)
                                                                                                     acGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGc                                                              <  2:38960/132‑1 (MQ=255)
                                                                                                     acGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGc                                                              >  1:38960/1‑132 (MQ=255)
                                                                                                             gcgGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGcc                                      <  2:103081/148‑1 (MQ=255)
                                                                                                              cgGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAg                                    <  6:52237/149‑1 (MQ=255)
                                                                                                                  cGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAg                                >  8:105312/1‑149 (MQ=255)
                                                                                                                   gatgatGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCCAGGAAGGCGCCGAACGCCCTGCTATGCTGCAGTCGACGGGCCAGGTAGc                               >  7:214821/1‑149 (MQ=255)
                                                                                                                         gaCCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTc                           >  8:34679/1‑147 (MQ=255)
                                                                                                                         gaCCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTc                           >  5:10251/1‑147 (MQ=255)
                                                                                                                          aCCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCa                         <  7:34679/148‑1 (MQ=255)
                                                                                                                            cgccgccgGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGc                       >  7:40380/1‑148 (MQ=255)
                                                                                                                              ccgccgGCAGGGTCAGCCAGAAGTTGGCGTGGGGCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGCCaa                    <  5:75278/149‑1 (MQ=255)
                                                                                                                               cgccgGAAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGCCaa                    >  6:19850/1‑148 (MQ=255)
                                                                                                                                              ccAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGCCAATCGGTGGCGGCGGCTa    >  8:105984/1‑149 (MQ=255)
                                                                                                                                                aGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGCCAATCGGTGGCGGCGGCTACg  >  8:280751/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCAGCATCTGCGCCGCAGTCTTCCACTTACCGAGGTTGGACACCGCCACATGCGCCCGCGCGCCTAGCTCGGCCATCCACTCGCGCAGCGCA‑‑GACACC‑ACGATCTCGCGGCCGATGATGACCGCCGCCGGCAGGGTCAGCCAGAAGTTGGCGTGGGTCTGCACCAGCAGCACCAGGGCCACTGCCACCATCAGCTTGTCGGCCACCGGGTCGAGGAAGGCGCCGAACGGCGTGCTCTGCTGCAGTCGACGGGCCAGGTAGCCGTCCAGCCAATCGGTGGCGGCGGCTACG  >  NC_002947/4632851‑4633140

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