Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 634,310 G→A 5.7% G242D (GGC→GAC)  ybdL → methionine‑oxo‑acid transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913634,3100GA 5.7% 224.0 / 2.2 70G242D (GGC→GAC) ybdLmethionine‑oxo‑acid transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTT  >  NC_000913/634163‑634457
                                                                                                                                                   |                                                                                                                                                   
ccGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCt                                                                                                                                                   <  1:105776/150‑1 (MQ=255)
 cGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTg                                                                                                                                                  >  1:239621/1‑150 (MQ=255)
       cGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGACTGGAAAGt                                                                                                                                            >  2:412440/1‑150 (MQ=255)
        gagaTTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTg                                                                                                                                           >  1:698194/1‑150 (MQ=255)
         agaTTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                                                                                                                                          <  2:341696/150‑1 (MQ=255)
         agaTTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTgg                                                                                                                                          <  1:408720/150‑1 (MQ=255)
          gaTTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTggg                                                                                                                                         <  2:1120513/150‑1 (MQ=255)
          gaTTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTggg                                                                                                                                         >  2:315854/1‑150 (MQ=255)
                tGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATtg                                                                                                                                   <  1:150496/150‑1 (MQ=255)
                    aTTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGtt                                                                                                                               >  1:212469/1‑150 (MQ=255)
                       aGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTgcc                                                                                                                            >  1:183658/1‑149 (MQ=255)
                       aGCGATGAAGTATACGATCAAATAAACTTTTCAAAACAGGGCCATGCCATTGTTCTGGCGCATACGAATATGCGTGAGCGGTAATTGGCGGTATCTAAATTTGGAAAGACATATAATATAACAAGCTGGAAAGTGGGTTATTGTGTTgcg                                                                                                                            >  2:1393716/1‑150 (MQ=255)
                        gCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTgcgc                                                                                                                           >  1:1053369/1‑150 (MQ=255)
                         cGATGAAGGCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCGGCGTGAGCGGGCAGGGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGcc                                                                                                                          <  1:238970/150‑1 (MQ=255)
                          gATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCa                                                                                                                         >  1:292821/1‑150 (MQ=255)
                             gAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGTAAATTTGGTTATTGTGTTGCGCCAgcg                                                                                                                      <  2:601899/150‑1 (MQ=255)
                              aaGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAgcgc                                                                                                                     <  2:1214506/150‑1 (MQ=255)
                              aaGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAgcgc                                                                                                                     >  2:1245697/1‑150 (MQ=255)
                                       aGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAgcg                                                                                                            >  1:1400991/1‑150 (MQ=255)
                                          acaTCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCg                                                                                                         <  2:407364/150‑1 (MQ=255)
                                           caTCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGACTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGa                                                                                                        >  2:491793/1‑150 (MQ=255)
                                            aTCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGaa                                                                                                       >  1:456825/1‑150 (MQ=255)
                                              cAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAAt                                                                                                     >  1:306072/1‑150 (MQ=255)
                                                 cTTTTCACCACAGGGCCATGCCAGTGTGCTGGCGAATCCCCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCg                                                                                                  >  2:1166530/1‑150 (MQ=255)
                                                 cTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCg                                                                                                  <  1:927934/150‑1 (MQ=255)
                                                 cTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCg                                                                                                  <  1:379605/150‑1 (MQ=255)
                                                 cTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCg                                                                                                  <  2:1181428/150‑1 (MQ=255)
                                                     tCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAg                                                                                              <  1:1387959/150‑1 (MQ=255)
                                                     tCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGACTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAg                                                                                              <  1:491793/150‑1 (MQ=255)
                                                       acaacaGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGt                                                                                            <  2:611838/150‑1 (MQ=255)
                                                        caacaGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCGTCAGCGCCGAAATTCGCAAGGTa                                                                                           <  2:426708/150‑1 (MQ=255)
                                                        caacaGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTa                                                                                           <  1:1393716/150‑1 (MQ=255)
                                                        caacaGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTa                                                                                           >  1:280008/1‑150 (MQ=255)
                                                           caGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACAt                                                                                        >  1:707076/1‑150 (MQ=255)
                                                            aGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATc                                                                                       <  1:519277/150‑1 (MQ=255)
                                                             gggCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCa                                                                                      >  2:1172773/1‑150 (MQ=255)
                                                              ggCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAg                                                                                     >  2:1115459/1‑150 (MQ=255)
                                                              ggCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGACTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAg                                                                                     <  1:412440/150‑1 (MQ=255)
                                                                 cATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTAt                                                                                  >  1:855043/1‑150 (MQ=255)
                                                                   tGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCt                                                                                >  1:1028179/1‑150 (MQ=255)
                                                                   tGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCt                                                                                >  2:931516/1‑150 (MQ=255)
                                                                   tGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCt                                                                                <  2:696999/150‑1 (MQ=255)
                                                                     ccAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGa                                                                              <  2:698194/150‑1 (MQ=255)
                                                                       aGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGAcc                                                                            <  2:863079/150‑1 (MQ=255)
                                                                          gtgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGCCCgtt                                                                         >  2:1039875/1‑150 (MQ=255)
                                                                                gcgcATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTg                                                                   <  2:456825/150‑1 (MQ=255)
                                                                                     tCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATAc                                                              <  2:855043/150‑1 (MQ=255)
                                                                                        gCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATAcccc                                                           <  2:707076/150‑1 (MQ=255)
                                                                                         cAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCg                                                          >  2:752260/1‑150 (MQ=255)
                                                                                             tGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGcc                                                                                                                    >  1:598595/1‑88 (MQ=255)
                                                                                             tGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGcc                                                                                                                    <  2:598595/88‑1 (MQ=255)
                                                                                                gTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGc                                                   <  2:212469/150‑1 (MQ=255)
                                                                                                   aGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTgg                                                >  1:605332/1‑150 (MQ=255)
                                                                                                           gTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCt                                        >  1:779650/1‑150 (MQ=255)
                                                                                                            tGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTg                                       <  2:779650/150‑1 (MQ=255)
                                                                                                                     cttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCTCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCGCAGCTGGCGCTTGCTGATATGCTAc                              <  2:1400991/150‑1 (MQ=255)
                                                                                                                     cttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTAc                              <  1:1172773/150‑1 (MQ=255)
                                                                                                                      ttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACg                             >  1:1350922/1‑150 (MQ=255)
                                                                                                                         aTTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGc                          >  1:1142271/1‑150 (MQ=255)
                                                                                                                             ggCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGaa                      >  2:1423560/1‑150 (MQ=255)
                                                                                                                              gCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAAc                     >  2:512959/1‑150 (MQ=255)
                                                                                                                                 aGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTg                  <  2:306072/150‑1 (MQ=255)
                                                                                                                                    ccTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGc               >  2:1319764/1‑150 (MQ=255)
                                                                                                                                      tATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTg                  >  2:497435/1‑145 (MQ=255)
                                                                                                                                      tATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTg                  <  1:497435/145‑1 (MQ=255)
                                                                                                                                       aTCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCAtt            >  2:551123/1‑150 (MQ=255)
                                                                                                                                           tatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCt        >  1:392830/1‑150 (MQ=255)
                                                                                                                                             tGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCTGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTg      >  2:503307/1‑150 (MQ=255)
                                                                                                                                               aCCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCg    >  1:1407095/1‑150 (MQ=255)
                                                                                                                                                 cGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGtt  >  1:802400/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTT  >  NC_000913/634163‑634457

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